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SECTION  I 


INTRODUCTION 

Project  7090,  Operations/Intelligence  Techniques  Experimenta- 
tion, hM  as  its  primary  objective  the  development  of  man/machlne  in- 
terface, graphics  display,  and  data  processing  techniques  for  effec- 
tive applications  of  intelligence  data  in  support  of  operational  mis- 
sions. The  problems  of  processing,  interpreting  and  applying  intel- 
ligence data  - data  from  multiple  reporting  systems  whose  output 
differ  in  content  and  format  - are  numerous.  Further  complications 
arise  in  trying  to  use  this  data  to  develop  a unified  consistent  real- 
time picture;  the  separate  systems  usually  report  in  different  time 
spaces  and  on  Independent  non-contlnuous  attributes  of  the  force- 
elements/  forces  involved.  A significant  portion  of  this  available 
data,  however,  is  reported  in  formatted  messages,  contains  primarily 
positional  information  and  has  time  and  positional  errors  which  are 
either  tolerable  or  can  easily  be  corrected. 

Fundamental  to  any  C^  system  handling  positional  geographic- 
based  data,  is  the  ability  to  display  this  data  over  a map  background. 

Existing  systems  possess  this  capability  to  varying  degrees;  however, 
none  of  them  has  the  total  set  of  needed  capabilities.  Systems  with 
an  air  situation  display  function  normally  offer  fast  response  to 
operator  requests  but  provide  only  minimal  geographic  background; 
alternatively  ground-oriented  systems  provide  the  desired  geographic 
data  (sometimes  too  much)  but  suffer  from  slow  response  to  "simple" 
change  of  context  coimnands  (e.g.,  translate  to  a new  area  of  the  map, 
zoom  in  on  the  current  center  point). 

A basic  map  display  system  must  not  only  provide  for  rapid  dis- 
play of  operations/intelligence  data  in  the  context  of  major  physical 
and  political  boundaries,  but  must  also  provide  for  the  display  of 
multiple  classes  of  feature  data  - e.g.,  lines  of  communication 
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(rivers,  railroads,  roads),  bases,  cities,  terrain.  In  addition,  it 
naist  automatically  display  this  background  information  in  a manner 
consistent  with  the  scale  of  the  area  being  displayed.  When  a large 
area  is  being  pictured  the  amount  of  feature  data  and  the  degree  of 
detail  in  the  basic  geography  should  be  minimized  - e.g.,  only  major 
rivers  and  cities  are  shown,  while  small  Islands  2tnd  Inlets  are  ig- 
nored. If  a small  region  is  being  shown,  then  all  available  boundary 
data  as  well  as  all  cities.  Islands,  and  Inlets  in  the  region  should 
be  displayed.  Finally,  the  map  display  system  must  provide  rapid 
response  to  the  "simple"  change  of  context  commands. 

For  FY76,  therefore,  a primary  objective  of  Project  7090  has 
been  the  design  of  a map  display  facility  which  we  call  a Geographical 
Data  Display  Environment  (GDDE) . A major  requirement  for  our  system 
was  a digitized  geographical  data  base.  A wide  variety  of  map  features 
were  desirable  in  this  data  base,  including  topography,  rivers,  roads, 
railroads,  cities,  and  military  bases  in  addition  to  the  usual  geo- 
political boundaries.  However,  only  the  geopolitical  boundaries — 
coastlines  and  national  boundaries — ^were  available  in  digitized  form 
at  the  outset.  Thus  our  development  of  a geographical  data  base  had 
two  parts:  the  immediate  utilization  of  the  digitized  coastlines  and 
boundaries  available  in  a data  base  called  World  Data  Bank  I (WDBI) 
and  the  development  of  a digitizing  capability  for  the  generation  of 
the  remaining  map  features. 

The  utilization  of  WDBI  was  constrained  in  several  ways.  First 
of  all,  special-purpose  data  structures  and  data  manipulation  tectr- 
nlques  were  likely  for  the  GDDE,  because  of  the  goal  of  providing  ra- 
pid response  to  "simple"  change-of-context  commands.  This  fact  implied 
that  the  WDBI  data  format  would  not  be  useable  and  that  we  would  have 
to  reformat  WDBI  data  in  GDDE  form,  A second  structure  involves  the 
amount  of  detail  that  must  be  left  in  the  maps.  Two  factors  Influenced 
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this  situation:  the  raster-scan  displays  we  used  and  the  special  em- 
phasis in  the  GDDE  development  on  the  relation  of  scale  and  detail. 

This  restriction  requires  very  finely  tuned  tools  for  controlling  the 
amount  of  detail  present  in  a derivative  map.  The  last  constraint  is 
that  we  must  be  able  to  provide  special  handling  for  various  anomalies 
that  arise,  either  from  the  data  Itself  or  from  our  manipulation  of 
it.  Altogether,  we  require  a special  set  of  tools  that  can  modify 
WDBI  to  create  custom  map  data  bases  for  Europe,  our  initial  area  of 
interest,  while  retaining  the  ability  to  focus  on  other  regions  of 
the  world  or  even  to  use  different  Input  data  bases. 

The  set  of  tools  we  required  was  Implemented  as  a package  of 
FORTRAN  programs.  This  report  documents  that  package.  Section  II 
provides  a general  description  of  our  process  of  generating  maps  for 
use  in  an  interactive  computer  graphics  environment.  Section  III  is 
a user's  guide  for  the  package,  and  Section  IV  is  a programmer's  guide. 
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SECTION  II 


THE  MAP  GENERATION  PROCESS 


INTRODUCTION 

The  process  of  generating  maps  for  use  In  an  Interactive  com- 
puter graphics  environment  has  five  components.  These  components  and 
their  inter-relationships  are  shovm  in  Figure  1.  A discussion  of  each 
component  is  given  in  this  section. 

DATA  CREATION 

The  data  creation  process  consists  of  obtaining  a set  of  geo- 
graphic coordinates  (latitude  and  longitude)  for  the  region  to  be 
mapped.  To  date  we  have  not  had  to  digitize  data  ourselves.  The  data 
creation  process  consisted  of  making  some  initial  modifications  on  an 
already  existing  data  base.  World  Data  Bank  I. 

The  data  contained  in  World  Data  Bank  I (WDBI)  was  digitized  in 
1973  by  the  Federal  Systems  Division  of  IBM.  There  are  over  100,000 
points  in  the  data  base,  and,  as  suggested  by  its  name,  the  range  of 
the  points  is  the  entire  globe.  The  data  is  organized  in  three 
files:  coastline  and  Islands,  boundaries  and  lakes,  and  an  index  to 
the  first  two  files.  The  record  format  is  given  in  Appendix  A,  Table 
I. 

In  this  format  World  Data  Bank  I was  unacceptable  for  use  in  our 
raster  graphics  lab.  Two  modifications  were  made  to  the  data  before 
the  tape  could  be  read  in  our  lab.  The  first  consisted  of  blocking 
the  data  at  6400  bytes  instead  of  80.  This  was  essential  to  work  with 
the  data  due  to  an  eight  tenths  inch  magnetic  tape  inter-record  gap.  The 
second  change  was  an  EBCDIC  to  ASCII  conversion.  Both  of  these  steps 
were  performed  on  an  IBM  370/158. 
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i EDITING 

DATA  REDUCTION  .. ►(CONTENT  EVALUATION) 

I 

IMAGE  CONSTRUCTION 


Figure  1.  Map  Generation  Process 
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The  data  could  then  be  handled  In  our  minicomputer  facility. 

The  contents  of  the  three  files  were  printed  so  that  an  estimate  of 
the  number  and  location  of  points  In  the  European  region  we  wished  to 
map  could  be  made.  On  examining  the  Index  contained  In  the  third 
file  It  was  found  to  be  Inadequate.  A thorough  sort  was  not  completed 
on  this  file,  and  the  Index  Itself  Is  neither  complete  nor  totally 
accurate.  WDBI  line  segments  are  referenced  by  a map  code  of  the 
area  (continent,  country,  region)  to  which  they  belong.  In  many  cases, 
however,  the  Index  would  be  needed  to  find  the  location  of  a line  seg- 
ment given  the  line  segment  number.  In  this  respect  the  Index  was  In- 
efficient. To  solve  this  problem  the  Index  was  Inverted  and  sorted 
by  line  segment  numbers.  (Its  format  Is  given  In  Appendix  A,  Table 
IVO 

For  all  subsequent  work  with  VfDBl  It  was  necessary  to  keep  only 
the  latitude  and  longitude  In  radians  for  each  data  point.  Reduced 
tapes  for  both  coastline  and  boundary  data  were  made.  For  each  data 
point  a 2-word  line  segment  number  and  the  latitude  and  longitude  In 
radians  were  retained  (see  Appendix  A,  Table  II).  The  data  was  con- 
verted from  ASCII  to  binary  and  remained  blocked  at  6400  bytes,  (l.e., 
400  16-byte  records  per  block) . The  Input  tapes  for  the  projection 
component  of  the  map  generation  package  are  the  reduced  coastline  and 
boundary  data  files. 

PROJECTION 

The  projection  process  consists  of  transforming  a set  of  geograph- 
ic coordinates  Into  Cartesian  coordinates.  Since  this  Is  a transfor- 
mation from  the  globe  to  the  plane,  that  Is,  from  three-space  to  two- 
space,  some  stretching  and/or  shrinking  must  occur  In  the  region  which 
Is  being  mapped.  The  stretching  and  shrinking  can  result  In  distor- 
tion of  the  shape,  area,  distance  and  direction  measurements  or  any 
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combination  of  these  features  In  the  mapped  region.  The  type  and  ex- 
tent of  distortion  Is  related  to  several  characteristics  of  the  region 
(e.g..  Its  size,  shape  and  location  on  the  globe)  as  veil  as  the  type 
of  projection  which  Is  used. 

For  very  small  regions,  that  is,  regions  extending  over  only  a 
few  degrees  of  latitude  and  longitude,  the  amount  of  stretching  and 
shrinking  Is  also  small  regardless  of  the  projection  which  Is  used. 

It  Is  possible  In  such  cases  to  simply  plot  geographic  coordinates. 

On  the  other  hand  for  a region  extending  over  many  degrees  of  lati- 
tude and  longitude,  projection  type  has  a definite  iiq>act  on  the 
resultant  representation  of  the  region. 

The  European  region  which  Is  mapped  for  the  project  Is  of  the 

O DO 

latter  type.  It  extends  from  12  W longitude  to  39  E and  30  N to 

O 

68  N latitudes.  A secant  conic  projection,  also  called  a conic 
with  two  standard  parallels,  was  chosed  for  this  region.  The  decision 
to  use  this  projection  was  based  on  several  considerations  including 
the  following: 

• simple  representation  of  parallels  and  meridians  for  ease  in 
judging  direction  between  points; 

• no  scale  distortion  along  meridians  so  the  region  to  be 
mapped  has  no  longitudinal  restriction  • 

• small  percentage  of  stretching  and/or  shrinking  along  the 
parallels  relative  to  other  projections  for  the  European 
region  (a  comparison  of  scale  error  among  several  projections 
Is  given  In  Appendix  B); 

• ability  to  control  where  minimum  and  maximum  distortion  oc- 
curs within  the  region  through  choice  of  standard  parallels; 

• simplicity  of  calculations  (details  on  the  mathematics  of 
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I the  projection  are  given  In  the  progranmers  guide  under  SUB- 

SET) . 

Figure  2 gives  a breakdown  of  the  projection  process.  There 
are  actually  three  steps  to  be  performed.  The  first  Is  to  make  a 
first  cut  on  the  world  data  base  to  obtain  a region  slightly  larger 
than  that  to  be  mapped.  Points  within  this  first  cut  are  projected 
Into  the  plane  via  a secant  cone.  The  final  step  scissors  the  quad- 
rangle so  that  It  forms  a rectangle.  A scope  photograph  of  the  pro- 
jected European  data  Is  sho\m  In  Figure  3.  A geographic  coordinate 
grid  has  been  overlayed  to  aid  In  determining  direction.  Because  a 
conic  projection  was  used  the  meridians  appear  as  radial  lines  from 
the  apex  of  the  cone  and  at  their  true  angle,  and  the  parallels  are 
concentric  arcs  with  the  apex  of  the  cone  as  their  center. 

Although  the  secant  conic  Is  the  only  projection  In  the  package 
to  date,  the  projection  component  can  be  modified  to  map  data  by  any 
projection.  Modifications  would  be  necessary  to  only  one  subroutine. 


DATA  REDUCTION 

The  third  component  In  the  map  generation  process  Is  data  reduc- 
tion. This  step  Is  essential  for  two  reasons.  The  first  Is  a limit 
on  the  amount  of  storage  available  for  map  data.  The  equipment  In 
our  graphics  lab  Includes  two  minicomputers,  an  Interdata  Model  4 and 
an  Interdata  Model  70.  It  Is  a result  of  the  minicomputer  environ- 
ment and  the  nature  of  the  application  (l.e.,  map  data  Is  essential 
only  as  a background  for  military  situation  overlays)  that  this  re- 
striction Is  Imposed. 

The  second  reason  for  reducing  the  number  of  points  In  the  data 
base  Is  the  limited  resolution  of  the  display  medium.  Figures  4-6 
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show  three  maps  of  Norway,  Sweden  and  Denmark  with  varying  data  base 
sizes  and  display  scales.  From  these  display  photographs  It  Is  clear 
that  In  order  to  produce  the  best  representation  of  a region  It  Is 
necessary  to  relate  the  size  of  the  data  base  to  the  scale  of  the  dis- 
play. If  the  scale  of  display  is  too  large  for  the  size  of  the  data 
base,  the  map  is  very  angular  as  shown  in  Figure  4.  In  addition  to 
its  angularity  many  fjords  along  the  Norwegian  coast  which  would  re- 
solve on  the  display  screen  are  not  present  at  all.  Equally  unaccept- 
able is  the  case  In  which  the  scale  of  display  Is  too  small  for  the 
size  of  the  data  base  (Figure  5) . The  result  in  this  situation  Is  a 
map  in  which  coastline  and  boundary  line  segments  are  not  distinct. 
Islands  shrink  to  a single  point  and/or  do  not  resolve  with  nearby 
land  masses.  If  the  size  of  the  data  base  is  appropriate  for  the 
scale  of  the  display  the  map  does  not  have  problems  with  angularity 
or  resolution  as  shown  In  Figure  6. 

In  order  to  produce  data  bases  of  the  proper  size  both  for  dis- 
play at  a given  scale  and  to  be  handled  by  our  system,  two  approaches 
are  taken.  The  first  is  to  remove  islands  which  either  shrink  to  a 
single  point  or  are  too  close  to  nearby  land  masses  to  resolve  at  the 
display  scale  unless  the  Island  has  significance  with  respect  to 
map  applications.  In  order  to  Implement  this  approach  there  Is  a 
routine  In  the  package  which  determines  both  the  size  of  an  Island  in 
a data  base  and  its  distance  from  surrounding  land  (see  documentation 
for  ISLAND). 

The  second  approach  consists  of  removing  points  from  coastline 
and  boundary  line  segments.  Points  should  be  removed  If  they  do  not 
resolve  at  the  display  scale  and  consequently  are  unnecessarily  being 
stored,  or  if  they  create  bunching  which  results  in  a less  than  op- 
timal representation  of  the  region. 

After  the  projection  process  the  Cartesian  coordinates  of  the 
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Figure  5.  Map  of  Scandinavia  With  Too  Much  Data 


data  polnta  are  stored  In  chains.  A chain  Is  defined  as  a group  of 
adjacent  points  forming  part  of  a coastline  or  boundary  line  segment. 
(The  format  of  the  data  Is  given  In  Appendix  A,  Table  111.)  Data 
reduction  by  point  removal  Is  accomplished  by  first  assigning  a rank 
to  each  point  In  the  data  base.  Then  for  a given  display  scale  a 
limit  Is  chosen  and  only  those  points  with  a rank  equal  to  or  greater 
than  that  limit  will  be  displayed.  One  of  two  methods  Is  used  to  as- 
sign  a rank  to  a point.  In  either  case  the  data  base  Is  processed 
on  a chaln-by-chaln  basis.  The  rank  of  a point  Is  a function  of 
either  the  deviation  of  the  point  from  the  general  trend  of  the  chain 
to  which  It  belongs  or  the  product  of  the  deviation  and  the  length  of 
the  trend  line.  The  first  method  Is  used  to  produce  a data  base  whose 
size  has  been  determined  solely  by  the  resolution  of  the  display  me- 
dium at  a given  scale.  The  second  method  was  added  to  handle  cases 
In  which  the  storage  restriction  further  limited  the  size  of  data 
base.  Under  the  first  method  points  with  equal  deviation  are  assigned 
equal  rank.  For  example,  In  Figure  7 points  a and  b would  each  be 
assigned  a rank  based  solely  on  the  distance  d.  However,  If  only  one 
of  the  two  points  could  be  retained,  point  b should  be  deleted  be- 
cause Its  deletion  results  In  a smaller  loss  of  area  than  would  a's 
deletion.  By  using  method  two  the  deviation  of  a point  Is  set  equal 
to  the  area  swept  out  by  the  point  and  Its  associated  trend  line  (for 
more  Information  on  the  point  ranking  algorithm  refer  to  the  program- 
mer ' s guide  under  DETAIL) . 

EDITING 

The  fourth  component  in  the  map  generation  process  Is  editing 
and  there  are  several  routines  in  the  package  dedicated  to  'cleaning 
up'  a data  base. 

Two  types  of  errors  were  found  In  the  raw  WDBI  data.  As  expected 

2 

The  method  which  was  originally  Implemented  for  rank  assignment  Is  in 
use  at  the  Harvard  Lah  for  Computer  Graphics  and  Spatial  Analysis  and 
la  described  In  "POLYVRT  User's  Manual",  Version  1.1. 
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In  creating  a large  geographic  data  base  digitizing  errors  were  made. 

A typical  example  of  this  type  of  error  is  unmatched  coordinates  of 
end  points  from  line  segments  of  adjacent  countries.  In  such  cases 
the  endpoints  must  be  changed  to  coincide  either  by  altering  the  co- 
ordinates of  an  existing  point  or  adding  or  deleting  a point.  The 
other  type  of  error  is  the  creation  of  artificial  endpoints,  that  is, 
forming  two  or  more  line  segments  from  an  Internal  boundary  or  coast- 
line of  a country.  In  order  to  create  a world  data  base  with  enough 
detail  so  as  to  be  able  to  map  small  regions  reasonably  well,  the 
maps  from  which  the  data  is  digitized  must  be  detailed  (l.e.,  large 
scale).  For  any  extensive  region,  an  example  Is  our  European  region, 
the  data  must  be  digitized  from  several  sheets.  Some  coastline  and 
boundaries  must  inevitably  be  split.  Because  of  our  storage  restric- 
tion and  the  overhead  associated  with  storing  a line  segment,  broken 
chains  should  be  joined. 

In  addition  to  correcting  errors  in  World  Data  Bank  I data, 
clean-up  routines  are  needed  to  produce  data  bases  tailored  for  a 
particular  display  scale.  For  example,  entire  chains  may  need  to  be 
deleted  which  represent  Islands  too  small  to  appear.  In  data  bases 
whose  size  has  been  determined  by  the  amount  of  available  storage 
rather  than  the  resolution  of  the  display  certain  display  anomalies 
infrequently  arise.  Figure  8 shows  two  examples  of  the  problem.  In 
(a)  deletion  of  a point  results  in  the  creation  of  a loop  within  the 
line  segment,  and  in  (b)  deletion  of  an  entire  inlet  results  in  the 
intersection  of  that  line  segment  with  a neighboring  one.  Due  to  the 
complexity  of  adding  an  algorithm  to  the  point  ranking  routine  to  pre- 
vent this  type  of  error  occurring,  and  due  to  the  infrequency  of  its 
occurrence,  corrections  on  the  rank  of  points  are  made  during  the  edit- 
ing process. 

IMAGE  CONSTRUCTION 

The  final  component  in  the  map  generation  process  is  image 
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construction.  In  order  to  be  displayed  the  chains  of  a data  base 
must  be  collected  as  an  image.  The  Image  can  be  displayed  via  our 
graphics  handler.  Pallet.  It  is  during  this  process  that  a minimum 
rank  for  points  is  specified  and  only  points  of  appropriate  rank  are 
placed  in  a line  of  the  image.  Each  line  consists  of  a set  of  coor- 
dinates corresponding  to  points  which  can  be  connected.  The  lines 
are  created  by  calls  to  Pallet  subroutines. 

SUMMARY 

The  package  of  programs  that  generates  geographic  data  bases  cov- 
ers the  four  processes  of  data  projection,  reduction,  editing  and 
image  construction.  Up  to  this  point  we  have  not  had  to  contend  with 
digitizing  data  because  this  process  has  already  been  completed  for 
coastline  and  internal  boundary  data.  Data  creation  has  only  implied 
making  the  initial  modifications  to  V/DBI  which  were  discussed  earlier. 
It  should  be  noted  that  data  creation  and  projection  are  only  performed 
once  for  a given  region.  This  is  in  contrast  with  the  last  three  com- 
ponents of  which  many  iterations  are  generally  needed.  For  example 
digitizing  errors  are  not  apparent  until  a map  is  displayed  and  their 
removal  is  not  assured  until  the  map  is  again  displayed. 

In  addition  to  the  routines  required  for  these  processes,  two 
others,ENDPT  and  PRINTM,are  included  to  provide  a listing  of  the  line 
segments  and  a printer  map  of  the  data  base.  During  some  processes 
these  routines  provide  sufficient  information  about  the  data  base  to 
decide  if  the  process  was  completed  properly.  In  such  cases  the  image 
construction  step  can  be  bypassed  which  significantly  speeds  up  the 
map  generation  process. 


3 

Pallet  is  the  graphics  display  language  designed  for  our  lab.  It 
consists  of  a collection  of  FORTRAN  callable  subroutines  which  allow 
the  user  to  display  data  and  interact  with  the  display  terminal  to 
change  the  display. 
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SECTION  III 


USER’S  MANUAL 


INTRODUCTION 

This  section  provides  the  instructions  for  using  programs  in 
the  map  generation  package.  There  are  eleven  routines  in  the  package 
covering  the  areas  of  data  projection,  creation  and  reduction,  image 
construction,  data  base  content  evaluation  and  editing.  A breakdown 
of  the  routines  in  the  package  by  these  areas  is  given  in  Figure  9. 

For  each  routine  a brief  description  is  given  as  well  as  the 

4 

user  inputs  and  program  output.  Appendix  A is  referenced  for  the 
input/output  data  specifications.  A sample  run  is  also  Included. 


4 

Logical  unit  to  external  device  assignments  are  listed  for  each 
program  for  each  device  Invoked  by  the  program.  For  devices  not 
used  by  the  program  it  is  asataaed  that  the  standard  initial 
assignments  have  been  made. 
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PROCESS 
DATA  CREATION 

PROJECTION 

REDUCTION 

IMAGE  CONSTRUCTION 
CONTENT  EVALUATION 

EDITING 


Figure  9. 


PROGRAM(S) 

BOX 

SUBSET 

DETAIL 

EXCLUD 

ISLAND 

CHAP 

ENDPT 

PRINTM 

JOIN 
. EDITDB 
MERGE 


Breakdown  of  Map  Generation  Routines 
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DATA  CREATION 


BOX 

Box  creates  a chain  of  points  which  provide  a rectangular 
outline  for  the  map.  The  coordinates  of  the  comers  of  the  rectangle 
may  be  specified  by  the  user,  or  calculated  In  the  program  to  be  the 
minimum  and  maximum  x-  and  y-  coordinates  of  points  in  the  data  base. 

Inatmctlons 
To  execute  BOX: 

• make  logical  unit  assignments 

Input  tape  on  drive  95 

output  tape  on  drive  85 
printer 

Carousel  terminal 

# load  BOX  from  loader 
LG  2E00 
LO  61  BOX 
END 

9 start  execution 
ST  5000 

Input 

Data;  projected  map  data  tape,  format  as  listed  In  Appendix  A, 

Table  III. 

User:  In  addition  to  specifying  the  logical  units  for  input  and  out- 

put tapes,  the  user  must  enter  the  coordinates  of  the  lower 
left  and  upper  right  hand  corners  of  the  rectangular  outline. 
If  the  user  would  like  these  coordinates  to  be  calculated  by 
•*  the  mlnlmta  and  maximum  x-  and  y-  coordinates  of  points 


AS  0195 
AS  0285 
AS  0313 
AS  0510 
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in  the  data  base,  correaponding  coordinates  of  the  tvo  points 
should  be  the  same.  A rank  within  the  range  specified  In 
DETAIL  for  the  points  created  by  BOX  must  also  be  entered. 

All  data  is  entered  from  the  Carousel. 


Output 

Data:  formatted  the  same  as  the  input  tape:  the  final  chain  on  the 

tape  is  the  box  chain  and  is  identified  as  such  by  a three  in 
the  third  field  of  the  chain  ID. 

Printer :at  the  conclusion  of  the  run  the  total  nusiber  of  chains  which 
have  been  processed  and  the  minimum  and  maximum  x-  and  y-  co- 
ordinates of  points  in  the  data  bMe  are  printed. 


Sample  Run 

iiT  0000 

CUTER  2 BOX  COOUDIUATE3  2E5.3 

-. 393-. 263 
0. 2730. 410 

ENTER  INPUT  DEVICE  HUMBER (UN) 
01 

ENTER  OUTPUT  DEVICE  NUMBER (UN) 
02 

ENTER  DEVICE  LEVEL  FOR  BOX 
11 


LOJ 
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DATA  PROJECTION 


SUBSET 

Subset  combines  three  processes:  projection  of  data  points  from 
the  globe  to  the  plane,  subsetting  the  data  base,  retaining  only  those 
points  within  a user-specified  area  and  reformatting  the  chains  which 
are  retained  for  ease  In  handling  in  subsequent  routines.  In  order 
to  decrease  processing  time  the  data  Is  first  subjected  to  a geographic 
coordinate  test.  Only  those  points  Inside  a specified  geographic 
range  are  projected  and  the  coordinates  thus  determined  are  checked 
against  a Cartesian  coordinate  range.  Those  points  Inside  this  range 
are  retained.  ^ 

Instructions 
To  execute  SUBSET: 

• make  logical  unit  assignments 

AS  0695  Input  tape  on  drive  95 

AS  0485  output  tape  on  drive  85 

AS  0313  printer 

AS  0104  card  reader 

• load  SUBSET  from  loader 
LG  2E00 

LO  61  SUBSET 
END 

• start  execution 
ST  5000 

Input 

Data:  WDBI  data,  format  aa  given  In  Appendix  A,  Table  I 

User:  The  user  must  specify  the  following  which  Is  read  from  four 

cards : 

^The  planar  region  describes  a rectangle  and  Is  referred  to  as  such  in 
the  sample  run. 
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• a single  digit  Integer  Identifying  the  type  of  data:^ 

1 ■ coastline 

2 ” boundary 

•two  three-digit  Integers  specifying  the  offset  into  the 
the  file  for  the  region  to  be  mapped,  and  the  number  of 
blocks  to  be  processed. 

• four  floating  point  numbers,  format  F6.2,  which  specify  the 
latitude  (North  to  South)  and  longitude  (West  to  East) 
limits  for  the  Initial  test 

• four  floating  point  numbers,  format  F8.5,  which  specify 
the  planar  limits;  min,  max  x-coordlnates,  followed  by 
min,  max  y-coordlnates 

Output 

Data:  The  output  data  from  SUBSET  Is  formatted  as  specified  In 

Appendix  A,  Table  III.  For  each  chain  of  points  there  Is  a 
header  for  Identification,  followed  by  the  Cartesian  coordi- 
nates of  the  elements  of  the  chain,  followed  by  an  end  of 
chain  marker. 

Printer :The  Input  data  which  was  read  from  cards  Is  printed  for  veri- 
fication. As  each  line  segment  Is  processed  diagnostics  are 
printed  which  include: 

• the  number  of  crossings  of  the  geographic  coordinage  range 
and  where  these  crossings  occur; 

• if  a line  segment  leaves  the  Cartesian  coordinate  range, 
where  it  leaves  and  If  It  re-enters,  where  it  re-enters; 


The  type  of  data  may  be  coastline,  boundary  box,  other.  For  appro- 
priate Integer  code  see  Appendix  A,  Table  III. 
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• the  number  of  generated  chains  from  a line  segment.  (A 
chain  consists  of  a series  of  points  inside  the  planar  re 
glon.  Thus,  from  a single  line  segment  in  WDBI  several 
chains  could  be  generated  if  the  line  segment  leaves  and 
re-enters  the  region.) 

At  the  conclusion  of  the  run  the  following  summary  statistics 
are  printed: 

• total  nuiri>er  of  line  segments  Inside  the  specified  region 

• number  of  line  segments  processed  and  totally  outside  the 
region 

• total  number  of  data  points  inside  the  region 

• number  of  generated  chains 

Ssmple  Run 

Card  Reader  Input 
2 

HH75 

071.00028.00-25.00040.00 

-0 . 39  80000 . 27800-0 . 2630000 .41000 

Printer  Output 

•UIMRFN  OF  blocks  SKtPPEO  s M 
number  of  blocks  to  BF  PERO  s /« 

I AT  PANC.F  71.PB  TO  P0.0B  I ONB  PAN6T  •25,00  TO  4P.00 
K RAMCF  IS  -M.19BP0  TO  0,27000  Y PAN6F  IS  -0,20300 
to' 0,41000 
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tINt  SEGME.NT  nUMBEH  ■300t3(«0 


LCAVING  RECTANCLL  AT  0.16460 
CROSSING  LAT-LONG  BOUNOANV  AT 
CROSSING  LAT«tONC  0OUNDAMV  AT 
iNfCRlNG  rectangle  AT  0.16S71 
LLAVING  rectangle  *T  0,1747R 
CNTtRlNG  rectangle  AT  0.193S4 


0.41000  FROM  0.16476 
0.124049671  01 
0.123S5607I  01 
0.41000  FROM  0.16867 
0.41000  FROM  0.17453 
0.41000  FROM  0.19334 


0.41086 

0.40147053E  00 
0.4926344SE  00 
0.40991 
0.41080 
0.40974 


TOTAL  M'MHFM  OF  C«0S3INC»  OF  LAT-L0M6  BOUNDS  • 2 

rniMBET^  OF  OFMF.HATtO  CHAINS  ■ 3 

NUMBER  OF  POINTS  INSIDE  hecTANRLE  ■SAP 

LAT  RANGE  IS  63.B1P  TO  71,P7ft 

LONG  range  is  7,b48  TO  il,Pt« 


LINE  SEGHfcNT  NiihHER  *3^01340 


TOTAL  NUNBFR  OF  CROSSINGS  OF  LAT^LONG  BOUNDS  ■ R 

number  of  generated  chains  ■ 1 

NUMBER  OF  POINTS  INSIDE  RECTANGLE  ■l!2 
LAT  RANGE  IS  B3.0PR  TO  05.937 
LONG  RANGE  IS  18.234  TU  24.P66 


F I ^Al. 
TOT  AL 

TOTAL 

TOTAL 

total 


SlMHARY 

1.3  INSIDE  AT  ALL  • 

U3  PROCESSED  BUT  TOTALLY  OUTSIDE  RECT, 
number  of  inside  pis  ■ 1579 
NMMBI R or  GENERATE^  CHAINS  ■ 55 


a 23S 
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DATA  REDUCTION 


DETAIL 

In  order  to  reduce  the  size  of  the  map  data  base  while  retaining 
an  adequate  number  of  points  for  an  accurate  representation  of  the  area, 
points  In  the  subsetted  data  base  are  assigned  a rank.  Then  only 
those  points  possessing  a rank  greater  than  or  equal  to  a specified 
value  are  displayed. 

DETAIL  measures  the  Importance  of  a point  In  delineating  a fea- 
ture and  assigns  a rank  based  on  that  measure.  The  Importance  of  a 
point  Is  measured  by  Its  deviation  from  the  general  trend  of  the  chain 
to  which  It  belongs.  The  calculation  of  the  deviation  of  a point  from 
Its  associated  trend  line  Is  the  first  task  performed  by  DETAIL.  Once 
the  deviation  of  a point  Is  known  a rank  Is  assigned  by  checking  this 
value  against  a series  of  user  specified  bandwldths  which  define  equiv- 
alence classes.  A point  becomes  a member  of  that  class  whose  band- 
width Is  the  largest  exceeded  by  the  deviation  of  the  point. 

Instructions 
To  execute  DETAIL: 

• make  logical  unit  assignments 

AS  0195  Input  tape  on  drive  95 

AS  0685  output  tape  on  drive  85 

AS  0313  printer 

AS  0510  Carousel  terminal 

• load  DETAIL  from  loader 
LG  2E00 

LO  61  SUBSET 
END 

• start  execution 
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ST  5000 


Input 

Data:  a projected  data  tape,  format  as  given  In  Appendix  A,  Table 

III 

User:  from  the  Carousel  the  user  is  asked  to  enter  the  following: 

• a single  digit  Integer  specifying  whether  or  not  the  user 
wants  diagnostics  printed  for  each  data  point: 

1 » yes 

0 ■ no 

• a single  digit  Integer  specifying  the  metric  to  be  used 
for  computing  the  deviation  of  points  from  their  associated 
trend  lines; 

1 * distance  of  the  point  from  the  trend  line 

2 > product  of  distance  of  the  point  from  the  trend  line 

and  the  trend  line  length . ^ 

• a two  digit  Integer  specifying  the  number  of  bandwldths 

• a single  digit  Integer  Indicating  the  bandwidth  specifica- 
tion method; 


bandwidth  (1) 

- BW(I) 

bandwidth  factor(l) 

- BWF(I) 

band  multiplier 

- BM 

minimum  Cartesian  coordinate  range 

- MCR 

bandwidth  Increment 

- INC 

minimum  bandwidth 

- MINBW 

Bandwldths  are  computed  by  one  of  three  methods: 

(1)  geometrical 

BWF(l)  - MINBW 
BWF(l)  - BWF(I-1)-BM 

^Originally  only  metric  1 was  Implemented.  The  'produce  metric'  was 
added  as  an  effort  towards  a more  even  point  distribution.  For  a 
description  of  the  two  metrics,  see  Section  II,  Data  Reduction. 
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BW(I)  - BWF(I)-MCR 

(2)  linear 
BWF(l)  - MINBW 

BWF(l)  - BWF(I-l)  + INC 
BW(I)  - BWF(I)*MCR 

(3)  arbitrary 

BWF(I)  Is  specified  by  the  user  for  all  I, 

BWF(I) > BW(I-l) 

BW(I)  - BWF(I)*MCR 

• for  geometric  bandwldths 

a floating  point  number,  format  E14. 7, specifying  mini- 
mum bemdvldth , MINBW 

a floating  point  number,  format  F5.3,  specifying  a band 
multiplier , BM 

a floating  point  number,  format  E14.7,  specifying  the 
minimum  Cartesian  coordinate  range  of  the  map , MCR 

• for  linear  bandwldths 

a floating  point  number,  format  E14.7,  specifying  mini- 
mum beuidwldth,  MINBW 

a floating  point  number,  format  E14.7,  specifying  band- 
width Increment , INC 

a floating  point  number,  format  E14.7,  specifying  the 
minimum  Cartesian  coordinate  range  of  the  map , mcr 

• for  arbitrary  bandwldths 

N floating  point  numbers,  each  having  format  E14.7 
specifying  bandwidth  factors  (N  ” number  of  bandwldths 
previously  specified  by  the  user),  BWF  (I),  I > 1,N 

Output 

Data:  The  third  and  fourth  fields  of  each  data  point  contain  the 

deviation  of  the  point  from  Its  trend  line  and  Its  rank  re- 
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spectively. 


Printer :Dlagnostlcs  are  printed  for  each  chain  which  has  been  pro- 
cessed including: 

• the  original  WDBI  line  segment  number  corresponding  to  the 
chain 

• Integer  code  for  type  of  chain  (coastline,  boundary,  etc.) 

• number  of  points  in  the  chain 

• point  distribution  by  bandwidth 

Carousel: A summary  is  given  for  each  bandwidth/equivalence  class 
including: 

• the  bandwidth  factor  and  bandwidth 

• the  number  of  points  assigned  to  that  class 

• the  number  of  points  assigned  to  that  class  and  those  of 
higher  priority 

• the  greatest  number  of  points  in  a single  chain  assigned 
to  that  class  and  those  of  higher  priority 

In  addition  the  total  number  of  chains,  total  number  of  points, 
the  number  of  points  in  the  longest  chain  and  the  nund>er  of  points 
in  shortest  chain  are  printed. 

Sample  Run 
Input : 

.sT  5000 

UCTAIL  UiJt'IliCK 
CNTCK  1 FOK  DEBUG,  ELSE  0 
0 

L.ri'EK  CODE  FOK  METKIC 

LNT.;;K  1 FOR  DEVIATION  FROM  TREND  LINE 

ENTER  2 FOR  PRODUCT  OF  DEVIATION  AND  TREND  LINE  LENGTH 
1 

ENTER  NO.  OF  DETAIL  LEVEES  (NN,  UP  TO  20) 

10 

ENVER  BANDWIDTH  SPEC  METHOD 
ENTER  1 FOR  GEOMETRIC 
ENTER  2 FOR  LINEAR 
ENTER  3 FOR  ARBITRARY 
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1 

ENTER  MIN  BANDWIDTH  (El  4. 7 ) 

0.0010000 

ENTER  BAND  MULTIPLIER  (N.NNN,OVER  1.0) 
1.  500 

ENTER  MIN  COORD  RANGE  (El  4. 8) 
0.68000000 


Output:  Point  Distribution 


S«q.  WDBIL.S.N. 

type 

Chain 

Chain 

#Pca. 

9 

1 

2 

3 

4 

5 

# 

bh 

% 

# 

27 

# 

62 

Ma  12 

6 

6 

6 

1 

0 

1 

b7 

1 

7M 

64 

**•  16 

6 

11 

2 

6 

0 

1 

Chms 

66 

1 

3t 

66 

H9  37 

M 

24 

2 

8 

0 

3 

6« 

^«»23S« 

1 

3? 

67 

Na  41 

6 

26 

3 

1 

4 

1 

r.H'ja 

76 

3el)23Q)l 

\ 

33 

66 

M*  17 

6 

13 

6 

6 

0 

Chh« 

71 

1IDI292II 

1 

36 

61 

NaJ28 

1 

23! 

16 

13 

16 

18 

CH«|* 

72 

3D*293(i 

1 

37 

62 

N*268 

6 

126 

21 

16 

16 

7 

rMM» 

73 

3itn294« 

1 

36 

93 

Hm  16 

6 

11 

6 

1 

0 

CHn* 

74 

3tlR2SS* 

1 

36 

94 

N8 

6 

6 

7 

6 

6 

6 

CHk« 

78 

3H0292I' 

1 

46 

98 

**•  12 

6 

7 

6 

1 

0 

1 

CNN* 

76 

3W4298H 

1 

41 

96 

N8  14 

6 

6 

1 

6 

• k 

0 

fMN* 

77 

3*82828 

1 

46 

161 

Nal97 

6 

128 

17 

IS 

9 

CN»'* 

76 

3882688 

1 

47 

162 

M*  42 

M 

22 

6 

6 

a 

3 

CHW* 

76 

3827718 

1 

46 

IMS 

Hm 

9 

6 

8 

6 

1 

6 

0 

fHN* 

66 

3882788 

1 

83 

166 

Hm 

7 

6 

4 

1 

6 

6 

0 

Summarization 


UL 

BWF 

BW 

NP 

CNP 

CLC 

CSC 

0 

0.00000 

0.00000 

2 

7550 

476 

4 

1 

0.00100 

0.00068 

48  19 

7543 

476 

4 

2 

0.00150 

0. 0 0102 

625 

2729 

201 

2 

3 

0.00225 

0.00153 

515 

2104 

142 

2 

4 

0. 00338 

0. 00230 

412 

1589 

104 

2 

5 

0.00506 

0. 00344 

319 

1177 

64 

2 

6 

0.00759 

0.00516 

203 

858 

42 

2 

7 

0.01139 

0. 0 0775 

166 

650 

28 

2 

8 

0.01709 

0.01162 

80 

434 

18 

2 

9 

0.02563 

0.01743 

54 

404 

12 

2 

10 

0.03344 

0.02614 

350 

350 

9 

2 

137  ClIAINE  PROCEEDED 
7550  POINTS  PROCESSED 
476  POINTS  IN  LONGEST  CHAIN 
4 POINTS  IN  SHORTEST  CHAIN 

DONE 
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EXCLUD 

EXCLUD  deletes  entire  chains  of  points  from  the  data  base. 

Instructions 
To  execute  EXCLUD 

• make  logical  unit  assignments 

AS  0195  input  tape  on  drive  95 

AS  0685  output  tape  on  drive  85 

AS  0510  Carousel  terminal 

• load  EXCLUD  from  loader 

LG  2E00 
LO  61  EXCLUD 
END 

• start  execution 

ST  8000 

Input 

Data:  projected  data  base,  format  as  given  in  Appendix  A,  Ta>^le  III 

User:  The  user  must  supply  a list  of  chain  numbers  corresponaing  to 

the  chains  to  be  deleted.  The  chains  must  be  listed  In  order 
of  their  appearance  on  the  input  tape. 

Output 

Data:  specific  r*'  .n'?  have  been  removed  from  the  data  base. 

Sample  Run 
oT  a 00  0 

l»Tc:k  input  device  numuer(nn) 

01 

I.NTEK  OUi’PUT  DEVICE  WUMBER(.'IN) 

06 
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L1JT2R  CHAINS  TO  HE  ELI  MI  NAT  KD  ( HNN  ) , 
CUTER  ZERO  AS  LAST  CHAIN 
057 
082 
161 
185 
003 
000 
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• ISLAND 


This  routine  is  an  aid  In  determining  those  Islands  which  are 
too  small  to  be  displayed  at  a given  scale.  An  Island  Is  deleted  If 
It  appears  as  a single  point  or  If  It  merges  ^.*lth  another  land  mass 
when  displayed.  ISLAND  computes  the  horizontal  and  vertical  extent 
of  all  chains  In  the  data  base  whose  endpoints  coincide.^  The  mea- 
surements are  made  In  terms  of  the  Cartesian  coordinates  of  the  data. 
To  test  the  proximity  of  an  Island  to  nearby  landmasses,  a second 
pass  through  the  data  Is  required.  The  coordinates  of  each  point  In 
the  vicinity  of  an  Island  are  compared  to  the  minimum  and  maximum 
X-  and  y-  coordinates  of  the  Island  calculated  In  the  first  pass 

through  the  data.  If  the  distance  between  these  points  Is  less 

g 

than  a minimum  value,  the  coordinates  of  the  points  and  their  re- 
spective chain  numbers  are  printed. 

Instructions 
To  execute  ISLAND 

• make  logical  unit  assignments 

AS  0683  Input  tape  on  drive  85 

AS  0313  printer 

• load  ISLAND  from  loader 

LG  2E00 
LO  61  ISLAND 
END 

• start  execution 

ST  5000 

Q 

Some  editing  Is  necessary  before  running  ISLAND  as  not  all  islands  in 
WDBI  have  coincident  endpoints.  ISLAND  will  not  perform  any  compu- 
tations for  such  chains. 

^This  value  was  chosen  In  relation  to  the  resolution  of  the  display 
medium.  With  256  x 240  addressable  positions  the  oilnimum  horizontal 
distance  was  4/256  and  the  minimum  vertical  distance  was  4/240. 
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Input 

Data:  projected  data,  format  as  given  In  Appendix  A,  Table  III 


Output 

Printer: From  the  first  pass,  for  each  chain  the  coordinates  of  the 

first  and  last  points  as  well  as  the  WDBI  line  segment  number 
and  chain  number  are  listed.  For  each  Island  the  minimum  and 
maximum  x-  and  y-  coordinates  are  listed.  From  the  second 
pass  the  coordinates  and  chain  nund>ers  of  points  less  than  the 
minimum  allowable  distance  apart  are  printed. 


Sample  Run 
First  Pass 


WDBLSN 

Chain  ty^ 

Chain  # Cub,  Chain  1 

x-atart 

y-atart 

x>and 

y*end 

1 

Its  171 

3,9783 

3, J748 

3,9788 

M.336I 

sttpstilfl 

1 

117  179 

3,1837 

3.1198 

3,9979 

8,2477 

1 

its  173 

3,1893 

3,1994 

3,1899 

8,1934 

Mm  V 

13 

p«i79e 

4T  V • 3,1833 

MAR  V 

rn 

PalSPP 

*T  V • 3,1918 

MfN  f 

13 

4T  » • 3,1737 

MBR  V 

IS 

3,1933 

AT  1 a •«10S4 

eriT* 

V • 

3,3139  nCLTt  V ■ 3,3144 

SHIAHS* 

1 

119  174 

3,1789 

3,1971 

P.17h5 

8,1971 

MIN  R 

m 

3,1731 

4t  T ■ 3,1934 

MAI  1 

ts 

3,1339 

»T  » • 3,1979 

Mm  y 

13 

3,1934 

4T  » • 3,1783 

MAR  V 

rs 

3,9334 

41  * • 3,1783 

nciT* 

1 • 

3,3133  0CLr4  V ■ 3,3373 

! 

193  173 

3,9997 

8,9749 

3,9783 

3,9733 

1 

191  178 

3,9733 

3,9988 

3,9994 

3,2793 

• M 

3 

1 177 

-3,9733 

•ff,400P 

Mm  I 

IS 

•3,4333 

4T  r • •3,9733 

MAR  R 

tn 

3,9933 

4T  V 3 .3,9733 

M|N  R 

IS 

•3,9733 

4T  X • .3,4333 

MAR  Y 

ts 

3,4133 

4T  X ■ 3,9833 

OUT* 

1 • 

B«6PPP  0C( 7A  y 0 P«60Pf 
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Second  Pass 


tfitsi  X 

Y Other  lecKl 

^ 0IA4B  chain# 

X 

y 

Delta  X 

Delta  y 

63 

M.H87S 

6.3961 

1 

6.6912 

6.3413 

6.6892 

63 

6.3901 

1 

6.6893 

6,3468 

6.6872 

63 

«>.M87e 

6.3961 

i 

6.6886 

6,3394 

6.6699 

145 

P,?746 

.6.1429 

4? 

6.2241 

•6.1378 

6.6669 

146 

8.7163 

.6.1389 

4? 

6.2241 

•6.1378 

6,6977 

145 

6.2169 

.6.1374 

4? 

6.2741 

-6.1378 

6,6664 

145 

6.7227 

.6.1439 

4? 

6.2241 

-6.1378 

6.6698 

;7 

6.6179 

6.2966 

in 

6.6166 

6.2497 

8.6671 

77 

6.6128 

6.2926 

in 

6.6166 

6,2497 

8.6638 

77 

6.61 79 

'«.?966 

in 

6. 6112 

6,2966 

6.6882 

77 

6.6126 

6,2928 

in 

6.6117 

6.2966 

8.6822 

77 

6.61 79 

6.2966 

in 

6.6114 

6,2913 

8,8699 

77 

6.6176 

6.7928 

in 

6.6114 

6.7913 

6.6619 

77 

6.6179 

6.2968 

in 

6.6121 

6.7926 

6,6648 

77 

6.6126 

6.7928 

in 

6.6121 

6.7926 

6.8668 

77 

6,6179 

6.2968 

in 

6.6137 

6.2924 

6.6849 

77 

6.6126 

6.2928 

in 

6.613? 

6.2924 

6.6864 

77 

6.6179 

6.2968 

in 

6.6142 

6.2927 

8.6641 

77 

6.6176 

6.2928 

50 

6.614? 

6.7927 

6.6888 

77 

6.6179 

6.2968 

50 

6.6148 

6,2927 

6,8641 

77 

6.6126 

6.2978 

50 

6.6148 

6,2927 

6.6666 
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IMAGE  CONSTRUCTION 


CHAP 


CHAP  creates  Images  of  the  data  points  via  subroutine  calls  to 
Pallet  and  stores  these  Images  in  a drum  file.  The  contents  of  this 
file  Is  properly  formatted  for  display. 

Instructions 
To  execute  CHAP: 

• make  logical  unit  assignments 

AS  0195  mount  tape  URG016 

AS  0685  input  tape  on  drive  85 

AS  0510  Carousel  terminal 

• load  COREDP 

RW  DE 
BI  DCOO 
LO  DE 

• load  URG016  from  tape  drive  95;  enter  these  commands  from 
Carousel  (COREDP  responses  not  given) 

ST  DCOO 

LO 

01 

0080,  FFFE 

• load  CMAP 

RW  60 
LO  60 

• start  execution 

ST  2E00 

Input 

Data:  projected  data, format  as  given  In  Appendix  A,  Table  III 

User:  the  user  will  be  asked  to  supply  an  eight  character  name  for 

the  Image  and  a two  digit  Integer  specifying  the  minimum  rcmk 
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Data 


(l.e.,  bandwidth  supplied  by  DETAIL)  which  points  must  have 
in  order  to  be  retained  in  the  map. 

Output 

: drum  file  containing  an  image  of  the  data  points  of  rank 

greater  than  or  equal  to  that  specified  by  the  user. 

Sample  Run 

ST  2L00 
ENTER  NAME 
MAPTHREE 

ENTER  DETAIL  LEVEL(NN) 

11 

.^TOP 

EOJ 
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CONTENT  EVALUATION 


ENDPT 


ENDPT  provides  a list  of  the  chains  In  the  map  data  base  and  the 
X-  and  y-  coordinates  of  the  first  and  last  points  In  each  chain. 


Instructions 
To  execute  ENDPT: 

make  logical  unit  asslgniaents 

AS  0685  Input  tape  on  device  85 

AS  0313  printer 

load  ENDPT  from  loader 
LG  2E00 
LO  61  ENDPT 
END 

start  execution 
ST  5000 

Input 

Data:  projected  data  base,  format  as  given  In  Appendix  A,  Table  III 


Output 


Prlnter:For  each  chain  the  coordinates  of  the  first  and  last  points 

are  listed  as  well  as  the  WDBI  line  segment  number  and  an  In- 


teger Indicating  the  source  of  the  chain  (l.e.,  coastline, 
boundary) 


Sample  Run 

TOTti  sro, 

14 

«e 
14 
fn 

t4 

3J 
14 
37 

3« 

43 


74T44T 

VilTAKT 

UND 

vmn 

4,1441 

4,3379 

4.4348 

• ^11114 

4,4184 

i«g|li99 

4,4484 

4.8984 

4.4448 

9, 4344 

4,4383 

4,4114 

‘•4,4146 

4.4311 

4.4384 

•4.8313 

PgP.MH 

4.3317 

4,4314 

4.4344 

4,1343 

4,3816 

4,4338 

4.4784 

4.3784 

•PgP40H 

8,3494 

4,3744 

4,3889 

4,3383 

1 1 9^ 

4,3437 

nin  to 

3441374 

3444M4 

J««4|4ll 

34*4 1 •» 

3444434 


smii 


CF 


04ICIN  srn. 
s 

14 

14 

14 

34 

34 

3« 

33 

3« 

37 

J* 

43 


3444444 

3444SI4 

34449S4 
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PRINTM 


PRINTM  prints  a crude  map  from  the  points  In  the  data  base.  Due 
to  the  low  resolution  of  the  printer  the  map  produced  by  PRINTM  Is  of 
limited  use.  It  does  provide  a check  on  the  extension  of  regions 
mapped  with  the  data,  and  most  major  problems  with  hose  (e.g., 

missing  chains)  will  be  made  obvious. 

Instructions 

To  execute  PRINTM: 

• make  logical  unit  assignments 

AS  0685  Input  tapes  on  device  85 

AS  0313  printer 

AS  0510  Carousel  terminal 

• load  PRINTM  from  loader 

LG  2E00 
LO  61  PRINTM 
END 

• start  execution 

ST  5000 

Input 

Data:  projected  data,  format  as  given  In  Appendix  A,  Table  III 

User:  The  user  Is  asked  for  a two  digit  integer  which  specifies  the 

minimum  rank  of  points  to  be  printed. 

Output 

Prlnter:The  output  from  PRINTM  Is  a 60  x 100  matrix  map  of  the  data 
base.  The  horizontal  and  vertical  extensions  of  the  map  are 
reflected  In  the  rows  and  columns  of  the  matrix.  An  element 
of  the  matrix  contains  an  X If  the  coordinates  of  any  point 
in  the  data  base  are  within  the  range  of  that  particular  cell. 
Otherwise  the  element  contains  a blank. 
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Carousel I At  the  conclusion  of  the  run  s PAUSE  Is  encountered.  The 
user  may  change  Input  tapes  and  continue  the  run  by  typing  CO, 
and  produce  the  map  with  an  overlay  at  the  same  level  of  detail. 
For  example.  If  coastline  and  boundary  data  were  on  separate 
tapes,  PRINTM  could  produce  a map  of  the  coastlines  first  and 
then  coastlines  with  a boundary  overlay.  The  coastlines  would 
be  Identified  by  the  character  'X',  and  boundaries  by  the  char- 
acter 'O'. 


Sample  Run 


ST  5000 

ENTER  DETAIL  LEVEL(NH) 
04 

PAUSE  99 
PAUSE 
RW  « 

ST  5000 

ENTER  DETAIL  LEVEL(RR) 
02 

PAUSE  99 

PAUSE  > 
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f 


yvvvvint 

y vvyy  vxyxvx 


yyy 

X 


xxxx 

XX 

xxxx 

XX  XX 


XXXXXXX  XXX 
XXXXX  X XXX 

X 

X 

XX 

X X 
XX  X 
X 

X X 
XX 
y 

XX 
X 


XXX 

X 


XX 

XX 

XX 

xxxx 

XXX  XXX 


xxxxxx 
XXX 
X X 
XX  X 
XX  X 
XX  X 
XX  XX 
XX  X X X 

vy  xxyy  y 

xy  xyy  x 

xyx  X 

XXX  XX 

XX  XXXXX  X 

X xxxx  X XX 

XXX  X X 

XX  X XXX 

XXXXX  XXX 
X XK 

XXX  XXX 
X xxxxxx  XX  X 
XX  XXXXXXX  X X 
X xxxxxx  XX 

XX  xxxxxxxxxxxxx 

XX  xxxxxx  XX  X 

xxxx  XXV  XX  V XX  XX 

XX  XXX  XX  XX 

VX  XXXXX  XX  XX  XX 

X r X r X X XXXXX  X 

X XX  X XX  XXX 

XX  XX  XX  X 

X X 


XXXXX 

xxxx 

X X 

X XX 

X 

X 

X 


XXX 

X 


X X 


XX 

XX 


XXX  X XX 
XXXXX  X X 
XX  XX  X 
XXX  X 

XXXXXXXX  XX 
XXXXXXXXX  X X 

XXX  xxxxy  X X 

V XXXXX  XX 
X XX  X 
XXX  X 
XXX 

XXVI  X 
XX 

XXX  X X XV 
XX 

XX  XX  X 
XX  XX 
XX 

xxx 
yy 

X 

xy 

XV 
VX 


X 
XX 
XX 

XX  X XX 

y XX  yyyvxxxvx  x 
XX  xxxx  X 

XX  XX  xxxxxx 
xxx  xxx 

XX  X IX  yy 

X XX 
XX  X 

X 

yxx 

xxxx  X yy 
X 

XX 
X 


y yy 
yyy 

X X XV 
XXXXXXX 
XX 


XX 


XX 

xxxxxx 

X 

X XX 
X X 


XVX 

xxx  xxx 
X XX  xxx 
XX 

XX  XX 
XX  X 


XX 


V 

XX 


yy 

V 


V 

yy 


X X 


y 

VX 


xxx  XX 

X 


xxxx  XX  X 


XX 


X X 
XX 
X X 
XX 

XX  X 
X 

X XX 
X 


X 

XX  XX 
XX 

XX  X 

X X xxx 
X X 
X XXXXX 
XX  XXXXX 
XX  xxx 

X xxxxxx 
X 

XX  XX 
X 


X X 
X 


XX 


XV 

XX 

XX 


XX 


X 

XX  X 
X 


X xxxx 
X 


XX 

XX 


XX 


X X 


xxx 

XX 


XX 

X 


X 

X X 
X 

XXV  XV 

X XX  XX 

xxxx  X y 
y vv 
viyy 


XV 


XX  X 
X 
X 

XX 


xxx 

X 


X 

xxx 

XX 

X 

y 

X X 


Figure  10.  Sample  Output:  Detail  Level  4 
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KIVW 


XMIX  VX 
XXXKX  X 


XXX 

X 


X xxxx 

XX 

X X XXXX 

xxxxx  xx 


XX 

XX 


XX 

X 


xxxxx  X X 
XXX  XX  X 
XXX  XX 

XXXXXXXX  XX 

XX  xxxxxxxxx  xxxx 

XX  XXX  XXXXXXX  X 
XX  XX  xxxxx  XX 
xxxx  X XX  X XX 
XXX  XXX  XX  XXX  X XX 

XXX  xxxx 

xxxx  X XX 

X XX  XX 

xxxxx 


XX 

XX 

XXX 

X 


XXXXXK 
XXX 
XX  X 
XX  X 
XX  X 
XX  X 
XX  XX 
XX  X XXX 
xxxxx  X 

XXX  X 

XXX  X 

xxxxx  XX 

xxxxx  X 

xxxx  X XX 

XXX  xxxx 

XX  xxxxx 

XXXXXX  XXX 
XX  XX 

XXX  X XX 
X XXXXXX  XXX  X 


X 

X 

xxxx 

XXX  X 
XX  X 

XX  X 
XX  X 
XXX  XXX 
XX  XX 
XX  XX 

XX  X 

XX  XX  xxxxxxxxxxx 
XXXXXXX  XX 

xxxxx  XXXXXX 

xxxx  XXX  X 

XX  X xxxxx 

X XX  X 

XX  X X 

XX  XX 

XXX 

XXXXXX  X X XX 


XXXXX  X 

xxxx 

X XXX 

X XX 

XX 

X 

XXX 

X 

xxxxx 
xxxx 
X xxxx 
XXXXXXXX 
XXX  XXX 


XXXXXXXX  XX  XX 


XXXXXX  XX 
xxxx  XX  X XX 
XXX  XX  XX 

xxxxx  XX  XX 


X 

XX 


X 

XX 

X 


XX  XX  X XX  XX 
XX  XX  XX  XXX  xxxxx 

XX  XXX  X X 

XXX  X 

XXX  XX 

XXX  XX 


XX  XX 

X X XXXXXXX  XX  xxxx 

XXXXXX  XXXXXX  XX  X 

XX  X XXX  xxxx  X 

xxxxx  XXX  X XX  XX 
XX  xxxx  X X 

XXX  XXXXXX  X X 

xxxx  XXX  XXXXXXX  xxxxx 


XX 

xxxx 

X 


X K 


XX 

XX 


XXX 

X 

XX 

XX 


X 

XX 

XX 

XX 


X XXX 
XXX 
XXXX 


X 

XX 

XX 


XX 

IX 

XX 

X XX  X 
IKX 

I XX 


XXX  XX  XXX 
X XX  XXX 
XX 

XX  XXX 

X XXX  X 

X XXXXXX 
IX 


XXX  XXX 

X xxxxx 
XX  xxxxx 
XX  XXX 

XI  XXXXXX 

xx  X X 

XXX  XI  XX  XX 

XX  II  X XX  X 


XXX  XX 
III 


xxxx  XX  I 


XXX 
XI  X 
XX 

X XX 
XX 


X I 
XXX 


IX 

XX 

IX 


xxxxx 

X 

XX 

X 

XXX 


XXX  XX 

XIX 
XI 

X 

XX 
X X 


X xxxx  xxxxx 
XXX  XX  X 
X X 
XIX 

XIX  xxxxx 
KXX  X X 
X XXX 

XX 


xxxxx 

XX  » 
XX 
I 

XXX  IX  X 
IX  I XXX 
X IX  X III 
XX  XX  XX 
xxxx  XX 
xxxxx  XX 

XXXXXX  X 
I XX 

xxxx 


Figure  11.  Soipl*  Output:  Detail  Level  2 
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EDITING 


JOIN 


JOIN  connects  chains  In  which  artificial  endpoints  were  created 
In  the  digitizing  process,  and  those  which  have  coincident  endpoints 
(e.g.,  coastline  of  neighboring  countries)  In  order  to  minimize  chain 
storage  overhead.  A maximum  of  20  chains  may  be  connected  In  a sin- 
gle run. 


Instructions 


To  execute  JOIN: 


make  logical  unit  assignments: 


AS  0195 
AS  0685 
AS  0212 
AS  0412 
AS  0510 


Input  tape  on  device  95 
output  tape  on  device  85 

drum  files  2 and  4 for  scratch 
files 

Carousel  terminal 


• load  JOIN  from  loader 


LG  2E00 
LO  61  JOIN 
END 

• rewind  drum  files 
RW  02 
RW  04 


• start  execution 
ST  8000 


Input 

Data:  projected  data,  format  as  given  in  Appendix  A,  Table  III 

User:  The  user  must  specify  the  Input,  output  and  scratch  files. 

If  several  runs  are  to  be  made  two  drum  files  and  one  mag 
tape  may  be  used  Instead  of  two  mag  tapes  and  one  drum  file 
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to  speed  up  the  process.  In  addition  the  user  will  Input  a 
series  of  four  digit  Integers  specifying  the  chains  to  be 
linked.  If  the  points  of  a chain  are  stored  In  reverse  order 
to  that  needed  for  linkage.  Input  a minus  sign  as  the  first 
digit  of  the  chain  number. 


To  Join  line  segment  0001  (A  B)  and  0002  (C  -►  B) , Input: 

0001 

-002 

0000 

Output 

Data:  User-specified  chains  have  been  linked.  The  new  linked  chain 

assumes  the  chain  number  of  the  first  chain  In  the  series. 

It  Is  placed  at  the  end  of  the  file.  The  third 
field  In  the  chain  ID  of  a Joined  chain  will  contain 
a four  (see  Appendix  A,  Table  III).  The  first  two  fields  of 
the  chain  ID,  whlth  hold  the  WDBI  line  segment  number 
of  the  chain,  contain  zeroes  because  a chain  produced  by  JOIN 
may  be  composed  of  points  from  several  WDBI  line  segments. 

Sample  Run 

Jr  8 0 0 0 

-NTER  INPUT  DEVICE  NUMBER  (UiJ) 

04 

ENTER  OUTPUT  DEVICE  NUMBER (UN) 

06 

ENTER  BCRATCri  DEVICE  NUMDEU(NU) 

02 
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LiLITUK  CilAI.jo  VJ  Dij  LI  i-iKED  ( lU-JUU  ) 

L'LGATIVL  FOA  OPPOSI'l'L  ORDER 

Ltli'ER  ZERO  .vS  END  OF  CHAINS  TO  BE  LINKED 
SELL  CHAIN  TO  EXIT 
1 1 
-23 
161 
0000 
..  rop 

LJJ 
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EDITDB 


EDITDB  functions  as  a point  editor.  With  It  the  coordinates  and/ 
or  rank  of  a point  can  be  changed.  Points  can  be  added  to  or  deleted 
from  a chain. 

Instructions 
To  execute  EDITDB: 

make  logical  unit  assignments 

AS  0195  Input  tape  on  device  95 

AS  0685  output  tape  on  device  85 

AS  0412  scratch  file  Is  drum  file  4 

AS  0510  Carousel  terminal 

load  EDITDB  from  loader 
LG  2E00 
LO  61  EDITDB 
END 

rewind  scratch  file 
KW  4 

start  execution 
ST  5000 

Input 

Data:  projected  data,  format  as  specified  In  Appendix  A,  Table  III 

User:  The  user  must  specify  the  Input,  output  and  scratch  files 

and  a series  of  four  digit  Integers  specifying  chains  to 
which  edited  points  belong.  Editing  Is  performed  on  a chain 
by  chain  basis.  Any  number  of  points  may  be  edited  In  a 
single  chain.  For  each  point  In  a chain  the  user  specifies 
• the  position  of  the  point  In  Its  chain  (Note:  Position  of 
the  point  must  reflect  previous  deletions  and  additions 
to  the  chain.  If  a point  Is  deleted,  positions  of  follow- 
ing points  are  decreased  by  one;  slmularly  If  a point  Is 
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added,  positions  of  following  points  are  Increased  by  one; 

• an  Integer  code  for  Che  type  of  editing  to  be  performed; 

# the  new  geographic  coordinates  and/or  rank  as  appropriate. 

Output 

Data:  User-specified  points  are  edited  and  the  chains  to  which  they 

belong  are  written  following  the  last  chain  In  the  data  base. 

Sample  Run 

<’  o 0 0 

i-'t.i'Kii  I i.Vi  '■  0 !•' V I v.’f’  M'.iLI' ( iJiJ ) 

•jf, 

i-.'.'iVK  OU ''I’J'i’  Of'  yU'i'  .iU.'.UKK  UJ.'I) 

K.i'i'i  i;  ok.'K/'  i'i'u  i>i  I ii.. ) 

■:  4 

.w  lit'  I * 1 4’ I’ 1) ; . 1 44 i’l ) 

of  KO  AJ)  I'nil  o*-'  Cltr>r:JS  ■'■’v;  tir  KDIT't'O 

0 0 0 4! 

00  ?i 
■:  0 0 •) 

I OiT  : 

1 -C  l4,‘.4.V;i  I L I 4J 

Jr  ill  ''Ail  IrViL 

"'-L.-'i'i. 

4 — ill)  »^.)  I iJ^l* 

' 1 -1 1 1 'Lr'i' I- 

I'llAI.  .JlTiUKH  2 

r .*ll*i<  r*.4  4“4  4tt  - .4  vlf  I^OIl.l  *0  Lyi*  l'ili**l'0(.4.Ji4) 

-50 

V.)  > 1 i JiJ  i'!LiMi>r'i<  =•  15'^  K iJ''’ r K 'i'/'3K  l4U.‘Uii^K  (-4) 

2 

r.Ji|-.<  Ll'Vr’:. 

0 1 

I Ui<'*'llr'rt  t*  D I i i 4^  S/'i-it'  t'll/' 1 .4  ? 
r .J  i'l- i<  0«i4O 

1 

I ..i  i H I'vJii  OK  »*oii4'p  I'o  lii'  Knr''i'()(.'uJ;i) 

151 

.■'UJ  i'iOi.  iJLMUKlJ  »=  151  KiJ'i’KK  '-VoK  .JUnUKI?  (J ) 

2 

I i4 ',r'\/KT-  (.4.4) 

0 1 
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FUKTUKK  I'UiriiJG  IlJ  SAi-ir'  C11AI;'J? 

KUTFK  l = jfKS,  0“iJO 

0 

OllAIiJ  iiUiiBFH  = ?1 

I'WTFK  P05ITI0U  Of  ro  UF  f DITFD  ( NNU ) 

f.7  4 

.-'OSiriOy  yUi'liiFR  = *>74  fi«ITF.(  TASK  WUHBFH  (N) 

5 

K U'r<Tlif;<  FDI'i’IlJG  ;.S  SAilF  CHAIiJ? 
fL4Tfl<  1-YfS,  0=NO 

0 

STOP 

FOJ 
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MERGE 

MERGE  combines  two  data  tapes  which  were  created  independently 
by  the  projection  routine,  SUBSET,  (e.g.,  coastline  and  boundary 
data) . 

Instructions 
To  execute  MERGE: 

# make  logical  unit  assignments 

AS  0195  output  tape  on  drive  95 

AS  0685  input  tape  on  drive  85 

^ Carousel  terminal 

® load  MERGE  from  loader 
LG  2E00 
LO  61  MERGE 
END 

• start  execution 

ST  5000 

Input 

Data:  projected  data,  format  as  specified  in  Appendix  A,  Table  III 

User:  The  user  specifies  the  logical  units  for  input  and  output. 

The  output  tape  is  initially  blank.  The  input  tapes  are 
copied  to  the  output  tape  one  after  the  other;  successive 
input  tapes  are  mounted  on  tape  drive  85,  £uid  the  logical 
unit  number  is  re-entered.  The  user  can  merge  any  number  of 
tapes  using  this  procedure. 

Output 

Data:  The  merged  tape  has  the  same  format  as  the  individual  input 

tapes  (see  Appendix  A,  Table  III).  However,  the  sixth  field 
of  the  chain  ID  contains  the  cumulative  chain  number.  For 
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example,  if  the  boundary  tape  was  copied  after  the  coastline 
tape,  the  new  chain  mmd>er  of  the  first  boundary  chain  would 
be  1 -t-  (the  number  of  coastline  chains).  On  the  Input  bound- 
ary tape  this  field  for  the  first  chain  would  contain  a one. 

Sample  Run 

0 r 5 0 0 0 

ENTLK  lUPUT  UbVlCt:  UUMBb.UUd) 

06 

bN'i'bK  OUTPUT  UbVICb  UUMDEMCjN) 

01 

LNTUK  INPUT  OeVICb  UUMDliR(UW)  - IF  FIMISlIIiU  CUTrJR  0 0 

00 

TOTAL  NUMBER  OF  CHAINS  172 
EOJ 
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SECTION  IV 


PROGRAMMER'S  GUIDE 


INTRODUCTION 

Additional  documentation  Is  provided  here  on  each  of  the  routines 
for  which  operating  instructions  are  given  In  Section  III.  The 
function  of  each  routine  and  a description  of  the  procedure  Implemented 
to  perform  that  function  are  stated.  This  Is  followed  by  a list  of 
the  common  blocks  used  and  subroutines  called  by  the  program.  A 
high-level  flow  chart,  program  listing  and  load  module  map  complete 
the  documentation.  Similarly  for  each  subroutine  called  by  the 
main  program,  the  function,  procedure  description,  common  block  and 
subroutine  list  and  program  listing  are  given. 


61 


PFECEDiro  paoe^blank-not 


PIIMCD 


BOX 


Purpose 

BOX  creates  a chain  of  points  which  form  a rectangular  outline 
for  a map. 

Procedure  Description 

The  routine  copies  all  of  the  data  to  the  output  tape  until  the 
end-of-data  mark  Is  encountered.  The  data  Is  read  Into  core  and 
written  to  tape  In  6400-byte  blocks.  In  core  It  Is  examined  four 
words  at  a time  for  the  end  of  the  data.  The  variable  IHEAD 
Identifies  the  section  of  a chain  being  examined (1 . e . , header,  set 
of  coordinates  or  end  of  chain  - see  Appendix  A,  Table  III) . When 
the  end  of  the  data  has  been  reached  the  five  point  box  chain  Is 
placed  In  the  buffer  and  written  to  tape.  A record  Is  kept  of  the 
minimum  and  maximum  x-  any  y-  coordinates  while  the  data  points 
are  being  copied  to  the  output  file  so  that  these  may  be  used  as 
the  box  coordinates  should  the  user  so  specify. 

Conmon  Blocks 

None 

Subroutines 

None 

The  end  of  the  data  Is  signalled  by  two  end-of-chaln  marks.  An 
end-of-chaln  Is  a four  word  entry  consisting  of  a floating  point 
1000.0  followed  by  three  blanks. 
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TBiirr  f 

TBIIFF  (I«3)aB 
Ial«4 

rrfi.LT.jfiOB)  Ro  m jh? 
aBTTE(TOllT)  buff 
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SUBSET 


Purpose 

SUBSET  performs  three  steps  in  the  map  generation  process; 
subsetting,  projecting  and  reformatting  VTOBI  data. 

Procedure  Description 

Each  data  point  In  V/DBI  Is  first  subjected  to  a gross  test  on  Its 
latitude  and  longitude.  The  geographic  (lat.,  long.)  coordinates  for 
this  test  are  input  by  the  user.  Only  those  points  inside  the 
specified  range  are  projected.  Once  a point  has  been  projected,  Its 
Cartesian  coordinates  are  tested  to  determine  whether  or  not  the 
point  falls  Inside  the  rectangle.  If  It  does.  Its  coordinates  are 
written  on  tape.  If  a rectangular  boundary  Is  crossed  in  drawing 
a line  from  the  point  previously  examined  to  the  point  currently 
being  processed,  the  coordinates  of  the  Intersection  of  this  line 
and  the  boundary  are  determined  and  written  on  tape.  Whenever 
a line  segment  extends  beyond  the  Cartesian  coordinate  range  of 
the  display  the  segment  Is  truncated  and  an  end-of-chaln  Is  written 
on  tape.  The  end-of-chain  marker  Is  also  written  after  the  last 
point  of  a line  segment  inside  the  Cartesian  coordinate  range  has 
been  processed. 

Conmon  Blocks 


Block  Name 

Contents 

Description  of  Contents 

/BLKl/ 

XYPTS  (1600) 

output  buffer 

/BLK2/ 

WLAT,  BLAT,  WLONG, 

ELONG 

geographic  coordinate  range 

/BLK3/ 

XMIN,  XMAX,  YMIN, 

YMAX 

Cartesian  coordinate  range 

/BLK4/ 

USTDP,  RUSTDP, 
CTRLME,  RGLOBE,  PI, 

PIl,  CONST. 

constants  for  projection 
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/BLK5/  IRITE  flag  signalling  that  an 

end-of-chaln  has  just  been 
written  In  output  buffer 

Subroutines 

PROJEC:  Projection  of  point  from  the  globe  Into  the  plane  via 

a secant  cone. 

CALL  PROJEC  (RLAT,  RLONG,  X,  Y) 

RLAT:  latitude  of  point  In  radians 

RLONG:  longitude  of  point  In  radians 

X:  projected  x-coordlnate , returned  by  routine 

Y:  projected  y-coordlnate,  returned  by  routine 

TESTl:  tests  the  geographic  coordinates  of  a point. 

CALL  TESTl  (PREVPT,  PRESET,  RLAT,  RLONG) 

PREVPT:  flag  signalling  location  of  last  point  processed 

with  respect  to  geographic  coordinate  range, 
returned  by  routine 

PRESET:  flag  signalling  location  of  point  currently  being 

processed  with  respect  to  geographic  coordinate 
range,  returned  by  routine 
RLAT:  latitude  of  point  In  radians 

RLONG:  longitude  of  point  In  radians 

TEST2:  tests  x and  y coordinates  of  projected  point. 

CALL  TEST2  (OLDFL,  NEWL,  X,  Y,  BOUNDX,  BOUNDY) 

OLDFL:  flag  signalling  location  of  last  point  processed 

with  respect  to  Cartesian  coordinate  range, 
returned  by  routine 

NEWFL:  flag  signalling  location  of  point  currently  being 

processed  with  respect  to  Cartesian  coordinate 
range,  returned  by  routine 
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X:  x-coordlnate 

Y:  y-coordlnate 

BOUNDX;  x-coordlnate  of  boundary  point  If  one  has  been 
generated  by  the  routine 

BOUNDY:  y-coordlnate  of  boundary  point  If  one  has  been 

generated  by  the  routine 

SAVELS:  writes  chain  Identification  header  In  output  buffer. 

CALL  SAVELS  (LSNl,  LSN2,  SOURCE,  NSUM,  POINT) 

LSNl:  first  three  digits  of  WDBl  line  segment  number 

LSN2:  last  four  digits  of  vn)Bl  line  segment  number 

SOURCE:  Integer  specifying  the  type  of  data  In  the  chain® 

NSUM:  new  chain  number 

POINT:  pointer  In  XYPTS,  the  output  buffer 

SAVEXY:  writes  x-  and  y-  coordinates  of  a point  In  output  buffer. 

CALL  SAVEXY  (X,  Y,  POINT)  or  CALL  SAVEXY  (BOUNDX,  BOUNDY,  POINT) 
X:  x-coordlnate  of  a projected  VfDBl  data  point 
Y:  y-coordlnate  of  a projected  WDBl  data  point 
BOUNDX:  x-coordlnate  of  a generated  boundary  point 

BOUNDY:  y-coordlnate  of  a generated  boundary  point 

POINT:  pointer  In  XYPTS,  the  output  buffer 

SAVEND:  writes  end-of-chaln  In  output  buffer 

CALL  SAVEND  (DUMMYX,  POINT) 

DUMMYX:  end-of-chaln  mark  • 1000.0 

POINT:  pointer  In  XYPTS,  the  output  buffer 
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01  ME  NS  TOM  RLOCK(  l6ffi«)»iniOCK(l«ef<),IXVPTS(l60P) 

CnHNON/BLKl /XYPTSCl ®00)  /BLK2/IILAT,HL*T,WL0NG,f LONG 

common  /RLK3/XMTN,XM4X, vmin, Vm*k 

COMMON  /RLK4/lJST0P,«USTr>P,CTRLMF,RGL08r,Bl,Pll, CONST 

COMMON  /BLK9/IR1TC 

rauiVALFNCf  fBL0CKn)#I8L0CKf  J)) 

EOUTVALFHCF  (XYPTSf  l)»IXYPT8(n) 

TNTECCP  CMAlN,GRTST»NEWFL.0L0FL,P0rNT,PRF3PT,PRrVPT 

integer  T0TCHA,T0T90m,sumin 

iNTtGCP  SOURCE 

DFAI  LARGL*»LAR6L0 

ft-TTE  «R 

I V 1 »(' 

I 

BOUNPX«R, 

BOUNOY«P, 

LSINbH 

ISOUTbO 

TOTSUMbO 

TOTCMAbr 

OUMMYXBIPIRCI, 

NSIIMBi) 

POINTbJ 
PI»3. 14159 
PIIbPI/IRP. 

USTOP»57,R0 

RSTOPbAI 

RGL0BE»1, 

CTRUMf  B7.47 
CT»LMC»CTRLMr*Pl 1 

rH)STOP»?,*Pl*RGLOBr  * t (II9TOP-R8TOPY7360,)b(  (COSfllSTOPBPIl  ) ) / 
CfCOSfflSTfiPBPIl  )-C0S(U3T0P*Pr  I ) ) ) 
fNTRLPB(MSTDP*R8TOP)/2, 

rnNStBRUST0P42.*P!BRGL0RF*(l'9T0P-CNTRl P)735P. 

C 

C REAO  from  CAROS  » • OF  BLOCKS  TO  RF  S«TPPEO»  • OF  BLOCKS  TO  RF  RFAO» 

CLaT-long  Range  (ofgrFFS)»  X^y  range 
C 

REA0(1,5RR)  SOURCE 

REAOn.ARP)  N,M 

RF  A0(  1 »41R)  IILAT.RI  ATjMLONC.FLOnG 
RFAO  (1,4?R)  XMIN,XMAX, YMIN,YMAX 
WRITE  f3,R3R)N 
WRIU(3,535>  M 

WRITE  n,S4R>  i/LAT,B|.AT,wlPNC,FLONG 
WRITE  (3»55R)  XM^n, XMAX , YMTN, YMAX 
Ul ATbOLATbPII 
Rl ATbRLATbPII 
WLONGbwi onG*PI1 
f IONGbEIONGbPII 
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LJ  O LJ  CVCV  U LJ  U 


SKIP  M BLOCKS 

IMN.ro. B)  GO  TO  ?6 

00  ?S  I«1,N 
PF*0  (6)  BLOCK 
CONTINUF 
on  25B 

r(E*0(6)  BLOCK 
on  ?3M  lat»l597,4 

J«I*? 

1 «T^3 

RL*T«flt.OCKrJ) 
r«L0NG»BLOCK(L) 

IFfLSNl.ME.IBLOCKd))  CP  TP  75 
TML8N2.FO.IBI  OCK(K)I  GP  TO  iPO 

PRPCEO'IRF  for  first  POINT  IN  LlNf  SFGMfNT 

75  ir(LSNl.EO.B)  GO  TO  100 
C 

C give  SFG»*LNT  Skmhary  foB  previous  line  segment 

c 

IMINPIT.FO.0)  LS0IiT»L80UT*1 
IMINBiT.EO.n  LSIN«LSIN*1 

TnTSl)H«T0TSIJN4SIIMIH 

IF(INBlTt.EO.B)  CH*In»CM*IN-1 
T0TCHA«T0TCH*4CH*IN 

wPTTFf3,430)  GRTST  ^ 

NRITF  (3,440)  chain 
WRITE  (3,450)  SUHIN 
Sm4lLA«ShalL*/PI  1 
I ARGLA«L*RGL*/PI1 
SMALLO»SHAl.L0/PI1 
LARGLO»L*RGLO/PI1 
wRrTF(3,K60)  3MAtL*»LARGL4 
WRITE  (3,570)  SHALL0,L*RGL0 

IF((IBITE,CO,0),*Nn,(NSLiH,GT.0))  CALI  SAVrNP(PUMHVK, POINT) 
100  LSN1»IBL0CK(I)  0 

LSN2»IBL0CK(K) 

SmaLLAaBLOCK(J) 

I ARGLAiBl nCK ( J) 

SMALl 0«B(  0CK(l  ) 

LARGLOwBI  OCK(L) 

INBITC»0 

iNBITcB 

GRTSTwP 

CHATNsl 

SUH1N»0 

PRFSPI»2 

PREVPTf? 

NFWFI a? 

OI.ofL«2 


72 


«RlTEn,46»l)  LSN1,13N2 

CALL  TCSTt  (PHtVPT,»»HE3PT,RLAT,RL0NG) 

CALL  PROJEC  (PLAT,«LONG,«» V> 

CALL  TE9T2  (OLOFL  » nLpA  L » *»  V*  »0'JNnx,BOUWDy  ) 

IF(NEPFL,EO,0)  Cn  TO  27S 

INBITil 

TNBITC-1 

N9IJMbN9UM«1 

TAIL  9AVEL9a3N|,L9**2»90ll»rE,N9UM,RniNT) 
TAIL  9AVCXT(X,V, POINT) 

S"**  riM«9U*<T*<4  1 
on  rn  ??5 
r 

c pancfniiiff  FOR  RFMAINIRR  P01mT9  of  line  9FCMCNT 
C 

t29  CALL  TE9T1  (PRf VPT,PRE9PT,RLAT,RL0NC) 

IF  (PRE9PT,E0#F»EVPT)  CO  TO  )5« 
GRT9T«G»T9T^t 
NRITPO.iSS)  RLAT,RLONG 
15P  IF rPRC9PT,C0,0)  GO  TO  229 

CALL  POOJFC  (R| AT,RLONG,X, Y) 

CALL  TE9T2(0LDFL,NF*tFL,X,V,R0lJNnx,R0UNDV) 
IF(NEwFL.FO,»')  CO  TO  2PP 
IF  fOLOFL.EO.n  GO  TO  179 
C 

C PRC9ENT  POINT  IN,  PRCVT0U9  POINT  OUT  (ENTERING) 

r. 

N9IINaN9tlN*] 

CAI.I.  9AVEL9(LSNI,L3N2,90IIRCE,N9iJM, POINT) 

CAI  L9AVEXY  (BOllNOX,HOUNOV, POINT) 
3UnIN»9UMTN*1 

NRITEn,4  7P)  BOIJNOX,  BOUND  Y,X,V 
C 

C*  PRESENT  POINT  IN 

r 

179  CALL  8AV1XV  (X,V, POINT) 

INBIT.1 

INBITC"! 

SIIMlNaSlIMIN^t 
GO  TO  229 
C 

r.  MRE9FNT  POINT  OUT 
C 

2PP  IF(Ol OFl.EO.B)  go  to  229 

r 

C I'RESFNT  POINT  OUT,  PRFWlOU9  POTNT  In  (LEAVING) 

C 

CALL  9AVI XV(BOUNnx,BOONOY,PnTNT) 
9UmIN»91)ntn*| 

CALI  9AVFN0(nU“NYX. POINT) 
r.MAINaCMATNat 

INBITl»P 

XRTTF  (),4RP)  «OUNnX,(iO(INOY,  X,  V 
275  I ARGL A«AM*xi fLARGLA, BLOCK! J) ) 
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SM all A«*min1 (SMALL BLOCK  (J)  ) 

LARGLOaAMAVl  (LARr.LOf  BLOCK  (L)  ) 

SMALLOaAMTNl  (SMAl  LO » BLOCK (D) 

r 

r I sn  nr  TfS^lNG  FOR  tmat  POI^jT 

, 1 ■ rri»JT  JwUE 
7!'/  rpMTTNUF 

f 

r(  NO  OF  TrSTI*Jr.  W ALL  POINTS 
CGTVf  riNAl  SFr.MCNT  SUMMART 

c 

CALL  8AVENO(IH'MMyX, POINT) 

CALL  SAVFNnfOUMMVX, POINT) 

NRITF(A)  XYPTS 

IF  (INBIT. CO. 0)  LSnUT»LS0UT4l 
IF  (INBIT. ro.n  LSINaLSIN  I 
TOTSUMbTOTSUMaSUMIN 
IF(INnlTC.FO.P)  CHAlN»CMAIN-l 
TOTCHAbTOTCMAaCHAIN 
mRITE(3,A3B)  GRTST 
WHITE  (3,440)  CHAIN 
wRITFf3,450)  Sl'MIN 
r 

r"RITF  final  summary 

c 

WRITE  (3,490)  LSIN 
WRITE(3,500)  LSOUT 
WRITF(3,S10)  TOTSUM 
WRTTE(3,520)  TOTCHA 
40M  FORMAT  (213) 

410  F0RMAT(4F6.2) 

420  F0RMAT(4FB.S) 

430  rORMAT(//47H  TOTAL  NUMBER  OF  CROSSINGS  OF  LAT-l.ONG  BOUNDS  a,^) 
440  F0RMAT(29H  number  OF  GFNFRATFO  CHAINS  a. 13) 

4*J0  F0RMAT(3BH  number  of  points  inside  RFCTANGLF  a, 13) 

460  F0Rmat(///?2H  line  Segment  number  a. 13, 14) 

470  F0RMAT(/22H  entering  RECTANGLE  AT,FB,5,3X,F8.5,6M  FROM  ,F8,5,3X, 
rFB.5) 

4"0  F0RmaT(/22h  leaving  RECTANGLF  AT  ,F8.5,3X,F8.B,6H  FROM  ,FB.!),3X, 
CFB.6) 

485  F0Rmat(/31H  crossing  LAT-LOnG  BOUNDARY  AT  ,2E20,B) 

490  F0RHAT(///14H  final  summary. /25M  TOTAL  LS  INSIDE  AT  ALL  a ,14) 

500  F0RMAT(/47H  TOTAL  LS  PROCESSED  BUT  TOTALLY  OUTSIDE  RECT.  a, 14) 

510  F0BMaT(/29m  TOTAL  NUMBER  OF  INSIDE  PTS  a,!S) 

520  F0RNAT(/35H  total  NUMBER  OF  GENERATED  CHAINS  a, 14) 

530  F0RMAT(27H  NUMBER  oF  BLOCKS  SKIPPED  a, 13) 

535  FORMAT(/30H  number  OF  BLOCKS  TO  BE  READ  a, 13) 

540  F0RmaT(/11H  LAT  range  ,F6.2,4M  TO  ,F6.2,13H  LONG  RANGE  ,F6.?,4M 
CTO  ,fh.2) 

050  F0Rmat(/12h  X RANGE  IS  ,FB,5,4H  TO  ,FB.5,12H  Y RANGE  IS  ,FB.5,4M 
CTO  ,rfl,5) 

060  F0RMAT(13H  LAT  HANGE  IS, FB. 3, 4H  TO  ,F8.3) 

0/.1  rnHMAT(14H  LONG  range  TS,F9.3,4H  TD  ,F9.3) 
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A 810  0C8A 
A sti  0C94 
I ,V  0000 


PROGRAMS  : 


8152 

8232 

. P 

83B2 

. MES 

8430 

. U 

9498 

COS 

94BC 

SIN 

959A 

AMOD 

95C8 

AINT 

9662 

$3 

9694 

. COMP 

973E 

97A2 

. RARG 

9770 

$8 

NTRY- 

POINTS  : 

7834 

TESTl 

78FE 

TEST2 

7E96 

S AVEND 

7FA4 

BUMP 

8232 

. P 

82BA 

.0 

8 4 3C 

.U 

8462 

. V 

94bC 

SIN 

9576 
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95C8 

AINT 

9 6 3A 

. 2 
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96BA 
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9770 

$8 

979A 

.5 

OMMOU- 

BLOCKS : 

E6  8E 

BLKl 

FF3E  BLK2 

FFCA 

ULK5 

82BA 

• Q 

836A 

.0 

8 4 5E 

.V 

346C 

01 

9576 
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9588 

A MAXI 

963A 

.2 

964C 

$1 

9 6BA 

. RRARG 

9 6FC 

$6 

9 7 9A 

.5 

979E 
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7BD4 

SAVEXY 

7CE8 

S AVELS 

8012 

PROJEC 
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0J 
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.O 

83B2 

, MES 

846C 

01 

9498 

COS 

9588 
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AMOD 
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$1 
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$3 

96FC 

979E 

$6 

. ZERO 

97  3E 

. RARG 

FF9E 

BLK3 

FFAE 

BLK4 

UNDEFINED: 

NONE 

LOADER 


XOUT 
LOADER 
VA  0400 
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Subroutine  PROJEC 


Purpose 

PROJEC  projects  a point  on  the  globe  into  the  plane 

Procedure  Description 

At  present  the  projection  in  use  Is  a conic  projection  with  two 
standard  parallels,  also  referred  to  as  a secant  conic'  (see  Figure 
14).  The  cone  Intersects  the  globe  along  two  parallels  of  latitude 
which  are  called  standard  as  they  are  projected  at  their  proper  scale. 
Between  the  standard  parallels  the  scale  along  the  parallels  Is  too 
small, resulting  In  a shrinking  of  the  projected  region.  Beyond 
the  standard  parallels  the  parallel  scale  Is  too  large  so  that 
regions  may  appear  stretched  in  the  upper  and  lower  latitudes. 

The  meridians  of  longitude  are  always  projected  at  their  true 
scale. 

The  first  step  in  the  projection  process  Is  to  determine  the 
projected  radius  of  the  upper  standard  parallel  (N'Q*  In  Figure  14). 

If  r»  radius  of  globe 

0°  parallel  « upper  standard  parallel 
0°  parallel  - lower  standard  parallel 
Then  N’Q'  - 2.r  g - coFi) 

and  the  distance  of  an  arbitrary  point  from  the  apex  of  the  cone  is 
given  by 

N*B*  - N'Q*  + 2nr  (^3=^) 

If  B*  is  in  latitude  01°. 

As  the  meridians  of  longitude  are  projected  at  their  true  angles 
a point  in  longitude  0°  (figure  15,  point  B')  will  have  coordinates 
X » N'B'  SIN/J'  where  ~ 0°  ~ N’M"* 

y - N'B'  COS  8-' 
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n’ 


STANDARD  PARALLELS!  UPPER  • PO 

LOWER  > UV 


Figure  14  CROSS-SECTION  OF  GLOBE  PROJECTED  ONTO  A SECANT  CONE 
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Figure  15  SECANT  CONE  DEVELOPED 


The  origin  of  the  map  can  be  adjusted  by  adding  constants  to 
the  X-  and  y-  coordinates.  In  PROJEC  a constant  set  equal  to  the 
distance  of  the  central  parallel^^  from  the  apex  Is  added  to  the  y- 
coordlnate.  Thus,  the  origin  is  the  Intersection  of  the  central 
meridian  with  the  central  parallel. 


Common  Blocks 

Block  Name  Contents  of  Block 

/BLK4/  USTDP 

RUSTDP 


CTRLME 

RGLOBE 

PI 

PIl 

CONST 


Description  of  Contents 

upper  standard  oarallel 

radius  of  upper  standard 

parallel 

central  meridian 

radius  of  globe 

3.14159 

3.14159/180 

radius  of  central  narallel 


In  order  to  Insure  a maximum  scale  error  of  less  than  or  equal 
to  IZ  along  the  parallels  between  the  parallels  40°N  and  60°N  the 
following  constants  are  set  In  the  main  routine: 

upper  standard  parallel  ■ 57°N 
lower  standard  parallel  “ 41°N 

In  addition 

RGLOBE  - 1.0 
CTRLME  - 7.47°E 


Subroutines 

None 


The  central  parallel  Is  the  mid-parallel  from  the  two  standard 
parallels. 
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auBI*OUTl*<t  PBOJEC  («LAT,(*L0NG,XC(»0NT.YC»0NT) 

c 

C THIS  SUSnoUTlNC  PROJECTS  A POINT  AT  LATITUDE  RLAT  (RADIANS),  LONGITUDE  RL0N6 
C(RAOIANS)  onto  the  PLANE  USING  A CONICAL  PROJECTION  MITH  ThO  STANDARD  PARLLELS 
C 

CTARIABLE  OEPINITIONS 

CUSTOPl  UPPER  STANDARD  PARALLEL  (DEGREES) 

CRGLOBEi  RADIUS  OP  GLOBE 
CCTRLHEl  central  MERIDIAN  (DEGREES) 

CRUSTDPi  RADIUS  OP  projected  UPPER  STANDARD  PARALLEL 
CPIll  3.14199/lBe. 

CCONSTI  constant  to  be  added  to  t coordinate  op  projected  point,  it  sets  the 

C ORIGIN  OP  the  graph  at  the  INTERSECTION  OP  THE  CENTRAL  MERIDIAN  hITH 

c the  central  parallel, 

CRPTPROl  RADIUS  OP  projected  POINT  IN  LATITUDE  RLAT 
CXCRONTI  X COORDINATE  OP  PROJECTED  POINT 
CTCRONTi  T coordinate  OP  PROJECTED  POINT 

C 

COMMON  /BLX4/UST0P,RUST0P,CTRLME,RGL0BE,PI,PI1, CONST 
rPTPR0«RU3T0P*2,aRGL08E«PI*((U3T0P-RLAT/PI1)/38»,) 

IP(RL0NG,LT,CTRLME)  go  to  IRR 

c 

C LONGITUDE  OF  POINT  13  GREATER  THAN  OR  EQUAL  TO  CENTRAL  MERIDIAN 
C 

*CRDNT»RPTPR0*3IN(RL0NG-CTRLME) 
rCRONT»-(PPTPRO»rnS(RLONG-CTRLME)) ♦CONST 
CO  TO  2BR 
C 

C longitude  of  point  13  LtSS  than  central  meridian 
C 

(014  XCRONT»-(RPTPRO*3IN(CTRLME-RLONC)) 

TCR0NT»-(RPTPR0»C03(CTRLME-RL0NG))^C0N3T 

200  continue 

RETURN 

END 
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Subroutine  TESTl 


Purpose 

TESTl  makes  an  initial  cut  on  WDBI  based  on  the  geographic 
coordinates  of  the  region  to  be  displayed. 


Procedure  Description 


The  latitude  and  longitude  of  a data  point  passed  as  parameters 
in  the  routine  are  tested  against  the  user-specif led  min  and  max 
geographic  coordinates.  A flag  is  set  to  signal  the  location  of  the 
point  with  respect  to  the  geographic  limits.  A flag  is  also  set 
identifying  the  location  of  the  point  previously  processed.  By 
checking  the  values  of  these  flags  upon  return  to  the  mainline  the 
current  status  of  the  chain  is  known.  The  four  possible  states  are 
(1)  chain  leaves  the  region,  (2)  chain  (re-)  enters  the  region, 

(3)  chain  remains  Inside  the  region,  and  (4)  chain  remains  outside 
the  region. 


Common  Block 

Block  Name  Contents 

/BLK2/  ULAT 

BEAT 

WLONG 

ELONG 


Subroutines 

None 


Description  of  Contents 

Maximum  latitude 
Minimum  latitude 
Western  longitude  limit 
Eastern  longitude  limit 
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subroutine  TeSTI(RRtVfT,RRESfT,RL*T,RLONG) 

C 

C THIS  SUBROUTINE  CHECKS  THE  L*TITUOE  *N0  L0N6ITU0C  Of  * POINT  *NO  SETS  « 
C(PRESPT)»1  IE  THE  POINT  IS  INSIDE  THf  RtMCE  *NO  *8  IP  THE  POINT  IS  OUTSIO 
grange,  it  also  sets  PREVPT»1  IE  THE  PREVIOUS  POINT  HAS  INSIDE  THE  RANGE 
CHS  IE  IT  HAS  OUTSIDE, 

C 

COHMON/BLKJ/ULAT,BL*T,hlONG,ELOnC 
integer  PREVPT.PRFSPT 
PRFVPThPRESPT 

IE((RLAT,LT,BLAT),0R,(RLAT,GT,ULAT))  go  to  IBS 
PRESPTH 

IE((RL0NG,6E.HL0NG)fAN0,(RL0NC,L6.EL0Nr,))  CO  TO  2BP 


IBS  PRFSPTh 

200  RETURN 

END 

I 

BLK? 

0010 

1 

ULAT 

0000 

t 

BLAT 

000A 

t 

HLONC 

000S 

f 

F.LONG 

000C 

K 

TESTt 

0024 

P 

TESTJ 

00BA 

1 

.0 

0000 

1 

,P 

0000 

r 

PRfVPT 

002A 

r 

PRESET 

002C 

f 

RLAT 

002E 

f 

RLONG 

0030 

A 

iriP 

00QE 

A 

?R0 

00AA 

84 


Subroutine  TEST2 


Purpose 

This  routine  tests  the  x-  and  y-  coordinates  of  a data  point 
against  the  user-specified  Cartesian  limits  of  the  map.  If  a chain 
leaves  or  (re-)  enters  the  display  region  and  a boundary  point  needs 
to  be  created,  the  coordinates  of  the  point  are  determined  here. 


Procedure  Description 

As  happens  In  TESTl,  the  Cartesian  coordinates  of  a data  point 
are  tested  again.st  user-specified  uln  and  max  Cartesian  coordinates. 
A flag  Is  set  to  signal  the  location  of  the  point  with  respect  to 
these  Cartesian  limits.  A flag  Is  also  set  Identifying  the  location 
of  the  point  previously  processed. 


These  flags  are  then  examined.  If  their  values  are  not  the 
same  the  chain  has  either  entered  or  left  the  desired  region  and  a 
boundary  point  needs  to  be  generated.  The  algorithm  for  generating 
a boundary  point  can  be  broken  down  Into  three  cases  (referring 
CO  Figure  16)  : 

(1)  only  one  coordinate,  either  X or  Y Is  outside  the 
range  (line  segment  AB) ; 

(2)  both  coordinates  are  outside  the  range  and  a vertical 
boundary  of  the  region  Is  crossed  (line  segment  CD) ; 

(3)  both  coordinates  are  outside  the  range  and  a horizontal 
boundary  of  the  region  Is  crossed  (line  segment  EF) , 


In  all  three  cases  the  first  step  Is  to  find  the  slope  of  the 
line  connecting  the  two  points  (X,  Y)  and  (X*  Y').  The  slope, 

M,  Is  given  by: 


M 


(Y*  - 

(X'  - X) 


Before  a point  can  be  generated,  we  need  to  knew  which  boundary  has 
been  crossed.  This  la  done  by  testing  the  Cartesian  coordinates  of  the 
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A-48,714 


RY;  BOTH  X-ANO  Y - COORDINATES  ARE  OUTSIDE  THEIR 


4 


* 


point  against  those  of  the  rectangular  boundaries . If  only  one  ^ 
coordinate  Is  outside  Its  range,  (case  (1)),  we  generate  the  proper 
value  for  the  other  coordinate  by: 

(1)  when  X'  Is  outside  its  range 

X"  Is  the  X-value  of  the  crossed  horizontal  boundary 
Y"  =.  MX"  - MX  + Y 

- M (X"  - X)  + Y 

■* 

(2)  when  Y'  is  outside  Its  range 

Y"  Is  the  Y-value  of  the  crossed  vertical  boundary 
X"  = (MX  + Y"  - Y)/M 


If  both  coordinates  are  outside  their  respective  ranges,  Y"  Is 
generated  by  equation  (1) . If  the  generated  Y"  is  not  Inside  the  Y- 
range,  Y"  A set  to  the  Y-value  of  the  crossed  vertical  boundary  and 
X"  Is  generated  by  equation  (2) . 

In  Figure  16  the  Y-coordlnate  of  B'  would  be  determined  from 
the  maximum  X-value.  The  Y-coordlnate  of  D"  would  be  determined 
from  the  maximum  X-value.  This  Y-value  exceeds  the  Y-range,  so  Y 
Is  set  to  the  maximum  Y-value.  The  X-coordlnate  of  D'  la  thus 
determined  from  the  maximum  Y.  For  line  segment  EF,  the  Y-coordinate 
of  F'  Is  determined  from  the  minimum  X-value  and  this  value  Is  with- 
in the  Y-range.  The  coordinates  of  F"  are  not  calculated. 


Conmon  B locks 


Block  Name 


Contents  of  Block  Description  of  Contents 


/BLK3/  XMIN 

XMAX 

XMIN 

XMAX 


Subroutines 

None 


Minimum  X value 
Maximum  X value 
Minimum  Y value 
Maximum  Y value 
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9tmBfl.lTt^ir  TrST2(0(.0^L,NfwFl.,X.Y,)ri  ,VH 

C0»M0N/«LK3/X*1IN,X«*X,VMTN,rH*X 
INTtGCH  OLO^LtNCWFL 
ntAL  H 
OLOFL»NtWfL 
C 

CNCN  lst 
C 

IFCNfwFL.OT.n  CO  TO  leo 
OLDXSaX? 

OLoraara 

c 

C TtST  X 
C 

\99  NeaFLal 

X?aX 

IF(X.CT,XM*x)  X2*xmAX 
IFtX.LT.XNlNI  xaaxMlN 
^ IffXa.NE.Xl  MEHFLafI 

C TEST  y 
C 

raay 

IFfV.GT.raAX)  raayMAX 
(y.LT.TNtN)  vaayMlN 
TFEVa.Nt.yj  NfaFLaF 
C 

C IF  AOTm  PBFVimiS  AND  CUHRENT  POIiyTS  At»r  IN  OM  OljT 

COF  IF  TME  POINT  IS  FIP3T  OF  L9,  PetUPN 
C 

IF((OLOFt,rO,NePFL)tO«,(OLOFL.6T.n)  00  TO  3tf 
C 

COCTEPMINE  BOUNOAPy  POINT 

c 

IFIOLOFU.EO.n  CO  TO  2BB 
C 

c SET  xa,ya  to  values  for  previous  point  cause  you  are  cnterins  rectangle. 

XaaOLDXa 

VaaOLOYa 

asR  ifeolox.ne'.x)  oo  to  aas 
X(ax 
Yjava 
GO  TO  see 

aas  Na(OLOY»V)/(OLOX«X) 

IF(xa,EQ.X)  00  TO  aO0 

yiaMa(-OLOXaXa)*OLOY 

Xtaxa 

IF((Yl.CC.YNIN),ANO.(Yl,Lt.YHAXn  GO  TO  SPR 

ase  xia(M#oi.nx-OLDy*ya)/N 
Ylaya 

300  Otoxax 
OLOYay 
return 
end 
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I 

Bl  K3 

BB12 

1 

XHIN 

22(12 

f 

X«** 

2024 

f 

VHIN 

222B 

E 

ypxx 

022C 

K 

TEST2 

2224 

P 

TE9T2 

2294 

1 

.0 

2222 

1 

22(12 

OUOFL 

2224 

NEpFL 

222C 

X 

222E 

r 

2232 

XI 

2232 

VI 

2234 

M 

2242 

inp 

226C 

OLOX? 

2244 

X2 

224E 

OLOY? 

2222 

V2 

22B6 

3B0 

2272 

209 

2129 

OLOX 

222E 

225 

2t2C 

OLOY 

22C2 

292 

2244 
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Subroutine  SAVELS 


Purpose 

SAVELS  creates  a header  for  chain  Identification. 

Procedure  Description 

The  header  for  a chain  has  length  of  eight  words  (see  Appendix 
A,  Table  III).  Because  positions  in  the  output  buffer  are  always 
filled  in  groups  of  four,  the  routine  can  create  the  first  four  words 
of  the  ID,  Increment  the  pointer  and  then  make  a single  check  In  the 
current  location  of  the  pointer  in  the  buffer.  If  the  buffer  is  full, 
it  is  dumped  to  tape  and  the  pointer  is  reset.  The  next  four  words 
are  filled,  the  pointer  is  incremented  and  tested  again.  The  buffer 
is  dumped  and  the  pointer  reset  if  necessary. 


Common  Blocks 
Block  Name 

Contents  of  Block 

Description  of  Contents 

/BLKl/ 

XYPTS 

output  buffer 

/BLK5/ 

IRITE 

signals  that  an  end-of- 

chaln  has  just  been 

written  in  output  buffer 

Subroutines 

BUMP:  copy  output  buffer  to  tape,  reset  pointer 

CALL  BUMP  (IPOINT) 

IPOINT:  pointer  in  output  buffer,  reset  by  BUMP 
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SltflROUTlNF  SAVFLSaS^tlrLSNJ^SOURCCfNLSNrlPOINT) 

DIMENSION  IXVPTSnfeBfn 

COMMON  /BLKI/  XVPTSU6PPI 

COMMON  /BLKS/IRITE 

FOUIVALtNCr  (XYPT8(l)»IXVPTSf  in 

iNTFGfR  SOURCE 

TXVPTS(IPOINn»LSNl 

IXVPT3nPOINT*n«LSN? 

I XYPTSn  POINT*?)  ■SOURCE 
IPOTNT»IPOTNT*it 

ir(TPOlNT,GT,t«0P)  CALI  flllMP(  fPoiNT) 

IXVPTS( IPOTNT)«NLSN 
rXYPTSf  rPOTNT*n»NLSN 
IXYPTS(IP0INT*2)bW 
IXVPTS(1P0INT*3)»P 
1P0INT*IP0!NT*4 

IF  nP0lNT,GT,15PP)  CALL  RUMP  ( IPOINT) 

IRITE-P 

RETURN 

END 


I 

Rl.Kl 

19M0 

c 

XYPTS 

PPM0 

F. 

IXYPTS 

0PK0 

I 

RLK5 

00P4 

F 

IRTTF 

K 

SAVEI  S 

00?4 

P 

SAVFLS 

019A 

1 

.0 

00M0 

1 

.P 

0000 

F 

l.SNl 

00?A 

F 

LSN? 

002C 

r 

SOURCE 

0P?E 

r 

NLSN 

0030 

F 

IPniNT 

0032 

1 

BUMP 

0000 

9} 


Subroutine  SAVEXY 


Purpose 

SAVEXY  copies  the  Cartesian  coordinates  of  a data  point  Into 
the  output  buffer. 

Procedure  Description 

The  routine  writes  the  x-  and  y-  coordinate  In  the  buffer  and 
updates  the  pointer  Into  the  buffer  four  words,  thus  saving  2 words/ 
data  point  deviation  and  rank  to  be  supplied  by  DETAIL 

Cornnon  Blocks 

Block  Name  Contents  of  Block  Description  of  Contents 

/BLKl  XYPTS  Output  buffer 

/BLK5/  IRITE  Signalling  that  an  end- 

of- chain  has  just  been 
written  In  output  buffer. 

Subroutines 

BUMP:  copy  output  buffer  to  tape,  reset  pointer 

CALL  BUMP  (IPOINT) 

IPOINT:  pointer  in  output  buffer,  reset  by  BUMP 
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SUB"OUrtMt  3AVeXV(FlRSr,SCC0N0,IB0INT) 

DIMENSION  IXVRTSdSBS) 
common  /BLKl/XYfTS(16Sf) 

COMMON  /BLM9/1MXTE 
COUIVAUENCriXVMTSiD.KYRTSttn 
C 

C ON  ENTRY  ME  ASSUME  THAT  IROINT* I , 8. 9, , . , 1 B*7  AND  THAT  IROINT  R01NT8  TO  T 
C NEXT  available  slot,  MUST  INCREMENT  IROINT  AfTfR  RXLLIN6  AND  MUST  DUMR 
C BUFFER  IF  FULL. 

C 

XYRTS(IR0INT)»F1R8T 
X YRTS  ( IROINT*  IXSECONO 
XYRTS(IR0INT*2)«B,B 
IXYRTS(IR0INT*3)»R 
IR0INT»IR0InT*4 

IF(IROINT,6T,16O0)  CALL  RUMR(IROINT) 

IRITE-S 

RETURN 

END 


I 

BLXt 

190B 

t 

XYPTS 

«B0B 

1 

IXYPTS 

000B 

I 

BLKS 

0RH4 

E 

IRITE 

0000 

M 

SAVEXr 

0024 

P 

SAVEXV 

00FL 

L 

.0 

0000 

1 

.r 

0000 

F 

first 

002A 

r 

SECOND 

002C 

r 

IPOINT 

002E 

i 

BUMP 

0000 
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Subroutine  SAVEND 


Purpose 

SAVEND  places  an  end-of-chaln  mark  In  the  output  buffer. 

Program  Description 

SAVEND  copies  the  end-of-chain  mark  (1000.0)  and  updates  the 
buffer  pointer  by  four  words. 

Comnon  Blocks 

Block  Name  Contents  of  Block  Description  of  Contents 

/BLKl/  XYPTS  Output  buffer 

/BLK5/  IRITE  Signals  that  an  end-of- 

chaln  has  just  been  written 
in  output  buffer 

Subroutines 

BUMP:  copy  output  buffer  to  tape,  reset  pointer 

CALL  BUMP  (IPOINT) 

IPOINT:  pointer  in  output  buffer,  reset  by  BUMP 
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subroutine  SAVENOCFXRST, IPOINT) 
DIMENSION  IXYPTSMSaa) 

COMMON  /SLKl/  XYPTSU6PP) 

COMMON  /HLKS/IPITC 

F0UIVAUENCF(XYPTSn),IXYPT8(in 

XYPT3(lP0TNn*FIRST 

IXYPTSdPOlNT^n  ■ 0 

IXYPTS(IPniNT*?)«P 

IXYPTS(IPOINT*3)«0 

IP01NT«IP01NT+4 

IF(1POINT,GT,1600)  CAU.  BUMP(IPOINT) 

IRTTE»1 

»<F  TURN 

END 

I BLKl  19M0 
C XYPTS 
f.  IXYPTS  MP00 
I 8tK5  MP04 
[.  TRITE  00M0 
K SAVEND  «»P24 
P SAVFNO  «0Ff> 

I ,Q  0PM0 

( ,P  00WP 
r FIRST  00?4 
r IPOINT  002C 
I.  BUMP 
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SUBPOUTIME  BUMPfIPOiNT) 
COMMON  /BLKl/XYPTSMfidP) 
WPTTEf4)  XYPTS 
IPOINT«1 

PFTURN 

ENO 

I BUKI  19(1(3 

E XYPTS  PCIM0 

K RUMP  eiP24 

P BUMP  (1P);S2 

I .0  PPfelP 

I .P  0PM0 

r IPOINT 
I PJ  001/10 
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DETAIL 

Purpose 

DETAIL  is  a point-ranking  mechanism.  By  assigning  a priority 
level  to  each  point  in  the  data  base  the  user  controls  the  amount 
of  data  to  be  displayed  by  specifying  a minimum  rank  which  a point 
must  have  to  be  retained  In  the  map. 

Procedure  Description 

Assignment  of  Rank.  The  projected  and  reformatted  data  pro- 
duced by  SUBSET  Is  organized  In  chains.  A chain  Is  defined  as  a 
series  of  adjacent  points  forming  a coastline  or  boundary  line  seg- 
ment. 

Data  points  are  assigned  a rank  or  detail  level  on  a chaln-by- 
[ chain  basis.  The  first  step  In  the  process  Is  to  assign  to  the  end- 

points of  the  chain  the  highest  rank.  A "trend  line"  Is  constructed 
between  these  points  and  the  deviation  of  each  data  point  from  the 
trend  line  is  computed  (see  Figure  17  (a)).  The  point  having  the 
greatest  deviation,  point  C,  Is  assigned  a detail  level  and  a new 
trend  line  Is  constructed  between  points  A and  C.  The  deviations  of 
all  points  between  these  endpoints  are  determined  (see  Figure  17  (b)). 
This  time  point  D Is  assigned  a detail  level  and  another  trend  line 
Is  constructed  to  Include  D.  This  process  continues  until  all  points 
between  two  endpoints  have  detail  levels  assigned.  When  this  occurs 
new  endpoints  are  chosen  such  that  the  first  endpoint  Is  the  last 
successive  data  point  with  a detail  level  already  assigned  and  the 
other  endpoint  Is  the  first  ranked  data  point  following  tha  newly 
assigned  endpoint.  For  example.  In  Figure  17  (c).  If  all  points 
between  A and  D had  detail  level  assigned,  and  If  both  of  the  points 
between  D and  C had  not  been  given  detail  levels,  D would  remain  an 
^ endpoint  and  C would  become  one.  If  the  maximum  deviation  of  points 

from  the  trend  line  Is  held  by  more  than  one  point,  all  points  with 
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C 


ENDPOINTS  A,  B HAVE  BEEN  ASSIGNED  THE  HIGHEST  RANK. 
C HAS  THE  GREATEST  DEVIATIAN 


FIRST  ITERATION 
D HAS  THE  GREATEST  DEVIATION 


SECOND  ITERATION 
E HAS  THE  GREATEST  DEVIATION 


Figurt  17  RANK  ASSIGNMENT 
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that  deviation  are  ranked  the  same. 


Rolntlonshlp  Between  Rank  euid  Deviation.  Each  rank  corresponds 
to  a band  whose  width  has  been  determined  by  the  user.^^  A point  Is 
assigned  the  rank  corresponding  to  the  largest  band  which  la 
exceeded  by  either  the  deviation  of  the  point  or  the  product  of  Its 
deviation  and  trend  line  length. 

It  Is  possible  that  the  deviation  of  a particular  point  In  the 
nth  Iteration  of  the  deviation  calculations  will  be  greater  than  that 
calculated  In  an  earlier  Iteration.  In  fact,  the  deviation  of  a 
point  In  the  n^^  Iteration  could  be  greater  than  or  equal  to  the 
deviation  calculated  for  the  endpoints  of  that  section  of  trend  line. 
For  example,  referring  again  to  Figure  17,  the  deviation  of  point  D 
in  (b)  Is  greater  than  that  calculated  for  point  C or  D In  (a). 

Thus  a restriction  on  detail  level  assignment  Is  necessary.  If  the 
deviation  of  a point  Is  so  great  so  as  to  assign  that  point  a higher 
rank  than  the  minimum  assigned  to  the  endpoints  of  Its  associated 
line,  the  point  receives  the  detail  level  equal  to  that  minimun . 

Connon  Blocks 

None 

Subroutines 

CDLEV:  calculates  the  deviations  of  points  from  a trend  line  and 

determines  the  point  with  the  greatest  deviation 

CALL  CDLEV  (X,  Y,  I,  J,  D,  P,  MD,  ILL) 


for  Information  on  band  specification  see  Section  III,  User’s  Guide 
to  DETAIL. 

13 

The  choice  of  using  the  deviation  of  the  point  or  the  product  of 
deviation  and  trend  line  length  for  rank  assignment  Is  made  by  the 
user.  For  Information  on  metric  choice  see  Section  II,  Data  Re- 
duction. 
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X:  array  containing  the  x-coordlnates  of  all  pointa 

In  the  chain  being  proceaaed 
Y:  array  containing  the  y-coordlnatea  of  all  pointa 

In  the  chain  being  proceaaed 
1:  pointer  to  the  first  endpoint  of  the  trend  line 

J:  pointer  to  the  final  endpoint  of  the  trend  line 
D:  array  containing  >1  for  each  element  yet  to  be 

assigned  a detail  level  and  the  rank  of  those 
elements  which  have  been  assigned  one 
P:  pointer  to  the  element  In  chain  with  current  maxiimi 

deviation  (l.e.,  to  the  element  to  be  assigned  a 
detail  level  upon  return  to  DETAIL) 

MD:  deviation  of  point  pointed  to  by  F 

TLL:  current  trend  line  length 
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I 
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Figurt  18  DETAIL  FLOWCHART 
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tNTPGCW  “i.MtPI 

C . nh,  nr  LfVFLS,  M|.Pt  ■ m|.*1 

TNTEfiCW  OFnijG 

INTEGER  I,J,P 

integer  R 

integer  L 

C 

integer  li.Jl 

INTEGER  12,13 

t 

NEAL  NTNAM 

C . MINIMUM  aANONlOTH  eactop 

REAL  Hi»(2n,BME(2n 

c . .SANOWIOTMS  ANf>  BANOWinTH  FACTHRS 

REAL  EACTON 

r . INGEPMETRir,  THE  HULTIPLIER 

r.  - IN  linear,  the  increment 

real  mcw 

r •>  MINIMUM  CPORPTNATE  range 

integer  mcppE 

C - metric  cope 

integer  RTP 

integer  TTV,TAPEI,TAPE0 

c - logical  units 

INTEGER  MCS 

r . NO  PE  chains  ikpiit 

integer  nps 

c - NO,  or  points  input 

integer  npi.c 

c . NO,  PE  POINTS  IN  IONGTST  chain 

Integer  npsc 

r - NO,  PE  points  in  shortest  chain 

integer  h 

c - HA*  CHAIN  length 

integer  n 

c - ACTUAI  chain  length 

REAl  »(B«PT,T(RPP) 

C - CHAIN 

REAL  OEVE8P0) 

c - PCWTATIONS 

integer  O(SRR) 

C - PE  I 4 II  levels 

INTEGER  NPfPl) 

c - NO,  nr  POINTS  AT  I EVEI  T 

integer  nc(?n 

r - longest  Chain  at  i evel  T and  ahpvi 

integer  CSCfRl) 

c - CHORTEST  chain  at  lEVEl  I AND  ABOVE 

integer  CNPr?n 

r . NUMsfp  nr  PIS  AT  LEVEL  I ANP  ABOVE 

integer  TCCI  (21) 

c - NO,  or  Pis  AT  LEVEI  I AND  ABOVE  IN  A CHAIN 

integer  TNPIPI) 

r . NO,  or  Pis  AT  levii  i in  a chain 
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iNTECr" 

c 'MAX  OLfv  tULnwcn  for  a point  in  a tmain 

REAL  IBl.lfrRf  1«P0),OBliFFP(1«HP) 
iNTcr.rR  iBuFFi(i60P)»ORurri(iAP0i 

COl'IVALFNCr  (iBUFFPd  ),TBUFFirin,(nmiFFim,OBUFFPm) 
c - INPUT  ANO  output  PUfFtBS 

integer  IRP,OBP»IBB»OBa 
C - BUFFER  pointers  *N0  SIZES 

INTFGFR  «c 

c . INPUT  B| OCX  COUNT 

INTFGFR  R,0 

C - R IS  fxctss  chain  length  on  input 

c - 0 IS  FL*6  for  -too  long  chain-  conottton 

HfFL  EOC 

C - END  OF  CHAIN  CODE 

REAL  MO 

c - MAX  deviation  or  PT  from  trano  line 

RFAL  ill 

C - TREND  I IMF  I EMGTH 

integer  LSNt.LCN?,aOURCF,8SN,NL8N 
C - HFAOER  data 

C 

C INITIAl  17ATION 
C 

PTH  • .1 
tty  • 5 
TAPEI  ■ 1 
TAPrr  • ft 
C 

C GET  input  PaHAMETFRS 
C 

mR1TE(TTV,9P0) 

900  FORMATOAHOFTAIL  OFUNFR) 
mRITF (TTV,B99) 

PR9  FORMAT (2BHENTER  1 FOR  OFBUG,  ELSF  P) 

RCA0CTTV,B9B)  OEBUG 
B90  FOHMAT(TI) 

C 

mRITE(TTV,9.M) 

931  format (piHFNTER  COOE  FOR  METRIC! 
b NRITF f TTY,93?) 

93?  FORhaTHTMENTER  1 FOR  DEVIATION  FROM  TREND  LINFT 
MRTTF (TTV,933) 

933  rORMATCSAHF NTER  9 FOR  PRODUCT  Of  DEVIATION  AND  TRFMO  LINE  1 EHCTHT 
HFA0ITTY,93A)  MCOOF 

934  FORmaTOO 

TF  ((MCOOF. NF  , IT. ANl>, (MCOOF  .NF,?!)  GO  TO  B 

r 

10  mHTTF(TTY,901T 

9E1  format  (AIHENTFR  no,  of  OtTAIl  LEVELS  (NN,  UP  TO  2P)T 
OF A0(TTy,902)  ml 
90?  FP»mat{I2) 

IF^hl.Gy.PO)  go  to  10 

TF  (Ml  .1  T.  1)  GO  TO  tP 
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MLPl  • ml*1 
C 

wRTTFfTTY,9^3) 

9M3  rOR*«AT(;»7HFNTEl»  RAmOwIOTM  SPEC  MFTHOO) 

2B  WRTTE fTTV»9R4) 

904  F0RMAT(21Mf NTtR  1 FOR  RFOMCTRIC) 

KR!TF(TTY,9P5) 

905  F0RM4T(1«HCNTER  2 FOR  LINEAR) 

WRTTF(TTY,O06) 

906  F0RHAT(21MFNTER  3 FOR  ARBITRARY) 

REAP(TTY,9R7)  I 

907  FORMATCII) 

IFd.EO.l)  GO  TO  30 
TFd.FO.Z)  GO  TO  4P 
IFfl.FO.T)  GO  TO  60 
GO  TO  20 

C 

30  wRTTE(TTY,908) 

90S  format  (29MENTtR  MIN  RANOiiIOTH  (EM, 7 )) 

RE  AO(T7Y,909)  minr-* 

909  rnRMAT(E14.7) 

RmF(2)  ■ MTNRw 
IFfMl.  ,E0,1)G0  TO  70 
60  MRITF(TTY,910) 

91»  F0RMATf39Hf nTFR  RANO  MULTIPLIER  fN,NNN,OVFR  1,0)) 

REA0(TTY,01  DFACTOR 
Oil  F0RmaT(F6,3) 

IFfFACTOR,LC,l,0)GO  TO  AP 
00  80  I»3,mlP1 

80  BMFfT)  ■ rACTOR*HwM  T-n 

on  TO  70 

c 

40  wRTTF (TTY,90B) 

RFAOf TTY,909)MlNnw 
Bt»F(2)  ■ mjnOw 
IF  fM|  .EQ.1 ) GO  TO  74 
10(1  MRITF  (TTY.91?) 

01?  FORMAT  (?BMf  NTER  INr.MfMFNT  (N.NNNNN)) 

RE *O(TTy,90O)  FACTOK 

IF  TF  ACTOR.lt, 0.0)  C’O  TO  100 

on  90  I«3,MIP1 

OF  Bi>F(I)  ■ F ACTOR^flUF  ( I«l) 

GO  TO  70 

r 

6F  MAT FF f tty. 01 3)  “L 

91J  FORMAT fSMF^'TER, 13, ?1H  BAMOmIOTMS  (n.nnNNM)) 

00  110  I«?,M|.P1 
J ■ T-1 

1?0  -MTE(TTY,914)  J 

014  f ORMATI IBHFNTER  BAMOmIOTH, I O 

t’F  AOf  TTY,9F9)BMF(t) 

TF  (I, to, 2)  GO  TO  114 
M fMMrf T).GT,BkFTT-l))  GO  TO  IIP 
GO  TO  1?0 
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tl0  CONTINUE 
C 

7»  WRlTf (TTV,915) 

915  rORMAT(29HrNTI.R  HIM  COORO  RANGE  fFH.Bl) 
REAO(TTY,909)  hCW 
00  13M 

130  RwfT)  ■ RWE(1)*HCW 
BUM)  u 0.0 

HwFM)  • 0.0 

C 

C TNITIA1.TZATT0N 
C 

NCS  a 0 
H»S  a 0 
NPLC  a 0 
NPSC  a 32767 


€••*****•**•***•*•« 

C********THIS  ShOULO  be  CHANGfO  NNfw  CMANGING  OINFNSTONS 
C**»*****FOR  *,V,0,PFV 

r************ 


Mb550 


NRITF(PTR.916)  •* 

916  FnRHAT(19MMA*  CHAIM  LENGTH  IS.T5I 
OP  14M  laJ.MLPl 
NPfl)  a 0 

n r(n  a 0 

M0  cscm  a 3P767 
r. 

MEAI»(TAPrn  IBIIFFR 

TRP  a -3 

PHP  ■ -3 

IMS  ■ 1600 

ORB  ■ 1600 

HC  • (I 

IOC  a 1000,0 

r 

C Gl  T THE  NEXT  r.HAlN 

C 

C 

C IBP  ->  LAST  4-HORP  mOC< 

C 

r Cf  T THF  I.SNl  4-WpRO  BLOCK 

r 

]60  IBP  a IBP«4 

TE  (THP.IT.IBSI  go  To  161 
r>FAO  fTAPEIl  TBUEFR 
IBP  a 1 

161  II  (IBMriRdBPI.FO.tnCl  GO  TO  201 
IBNI  a IBUFMdOP) 

IBN2  a IBiirriCTBPai ) 

50"wcE  a iniiFEKiRP*?) 

T1  a I"f'rrT(IRP*31 
r GET  TM»  I BN?  4.aOpp  BLOCK 
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r. 

T«f»  ■ TBP+4 

If  CTBP.LT, IBS)  GO  TO  1*«? 

PE*n  (T*PF1)  IBUFFB 
IPP  ■ 1 

lb2  SS^i  ■ TBUFFICIBP) 

*IL3N  • TBUFFICIBP*)) 

T2  ■ lailFFTCiBP*?) 

1.1  • IBtiFFTCIBP*3) 

r 

f PASS  TT  TO  TMF.  OUTPUT 

C 

OBP  • OBP^A 

IF  COHP.LT.IBBP)  go  to  )63 
HRITF  CTAPFO)  OPUFFB 
OHP  ■ 1 

163  OBi.irr  j (fiBP)  ■ l3*J1 

OBIIFFI  (OBP*n  ■ LSm2 
ORUn  ICn8P*2)  ■ SOUt^CF 
0RUFFI(0BP*3)  • T) 

OBP  ■ ORP*A 

IF  CPPP.LT'.OBS)  r.O  To  1(S4 
w«ITF  (TAPFO)  OBUFFM 
OBP  • ) 

Ifid  nnuFFTCORP)  ■ SS'J 

OBMFF I CORP*) ) ■ Ml  SN 
OMilFf  T C0BP*2)  ■ T2 
0RijFrnoBP*3)  • 13 
p ■ P 
N ■ P 

n « P 

c 

l7p  IBP  ■ TBP*d 

TFCTBP.l T.IHS)  go  TO  IBM 
rtFAOf TAPFT)  IBUFFP 
IMP  • 1 
HC  » BF.*1 

IFCROOCHC,  IP). TO  )BP 
(.PTTFCPTM.O)  7)  Rf 
017  FnPMATCRHHLOCK.ld) 
r 

IBP  trCIBUFF»CTRP),FO,FOC)  GO  TO  2PP 
TFfK.CO.M)  GO  TO 
M 

)R5  VC‘I)  ■ IBHfrP(tBP) 

VC^i)  • TBDFFBCTRP*)  ) 

GO  TO  17P 

)0P  0 ■ 1 

t>  • B*! 

GO  TO  )RB 

r 

2<)p  rrco.^'f  ,n  go  to  ??d 

0 • Nc;S*l 

-PTTF  CPTP.01R)  Or.’J.B 
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018  FORNAT(5HHLnCK, I4,«h  CM*Im,IS,7h  FXCIS8»I8) 

GO  TO  2?8 

r 

r. 

C ».Nf>  OF  INPUT  file 

r 

201  08P  ■ nnp*4 

IF (onP,LT,0R8)G0  TO  238 
HRITF(T»PF01  OBUFFR 
OBP  a 1 

230  0BUFFR(0BP1  a FOC 
nBUFFKOHP*!)  a P 
OB'IFFI  (0BP*2)  a P 
nRUFFI(08P*3)  a p 
aPlTF(TAPFr)',  OBUFFR 
CALL  COF(TAPCn) 

C 

CNPfHtPl)  a NP(MLPt) 

1 a ML 

2^0  IFd.rO.PIRO  TO  25P 

CNP(I)  a NPrn«CNP<l«n 
I a T-1 
r.n  TO  240 

r 

260  MPrTFFTTT.OlO) 

010  f 0R“AT(2H0L»  l**3HBi«iF,5lf,?HRM,ftX,?MNP,4K,3HCNP,3X,3HCl  C»3X,lHCSn 
no  26M  lal.MLPl 
J a T-1 

?f>U  wRITF(TTV,O20)  J»nMF(n,RPfT),NP(T).CNPn>»CLC(n»C8Cm 

020  F0RMAT(I2,  iX,F7,«.,1>‘»E^.S»»*»I^»f  *.I5»1*»IS,1X,IRT 
MRlTFrTTY,071)NCS 

021  FORMATdSflXflftHCMAXNS  PROCfSSEP) 
i«iRITE(TTV,022)NPS 

022  FnRMATfIB,1X»lftHP0TNT3  PMOCESSFO) 
wRTTFf TTV,0231NPlC 

023  F0PMAT(tS,lX»?3HP0TNT8  IM  LONGFST  CHAIN) 
wRTTKTTV,024)nPSC 

024  FnRMATfIB,l*»2AHP0TNT8  IN  ShORTFST  CHAIN) 
taRJTF fTTY,025) 

026  F0PMAT(4H0nnf) 

CAIL  PEafTAPCn 

CAI  I.  OfNfTAPFO) 

STOP 

r 

r HAVL  A Nf.M  chain  in  X ANO  Y 
r 3CT  r.l  ORAL  8TATI8TTL8 
r, 

220  NCO  a Ncfl^t 
NPO  a NP!)*M 

IF /N.OT.NPI C)  NPI C ■ N 
IF  (n.I.T  ,MP5C)  NPSCaN 
r 

f INITTAI I/f  FOR  CHAIN  PRUtFBSINR 

r 
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00  27^*  t«l 

?/(l  Of!)  ■ -i 

om  • ►*(. 

orw)  ■ “I 

OfV(l)  • 

OEV(N^  ■ 

00  2«t‘  T«1,MLP1 
?«H  TNP(n  • (» 

TNPfMi  pn  ■ 2 

1 • 1 
P ■ N 

c 

r riNp  TMT  MFVT  T J 

r 

?9P  .1  ■ P 

ir  (,i.r.T,fi  + i ))  on  to  jop 
31d  1 » T*1 

irn.EO.M)  r.o  TO  5.i»’ 

IF (Of T) .NF.-l ) GO  TO  JiP 
T « I-t 
J ■ W2 

32P  IF  (0(3)  ,MF  .-n  GO  TO  JPP 

J ■ 

r.n  TP  32P 
t 

C PFOFOPM  TPEnp  Lt'^f  C*LCO(  ATTPM 
( 

3PH  CPNTIMIIF 

IF(OFPljr,,FO.P)  GO  TO  3P1 
'»PTTF(PTP,«)Oft)T,J,P,M 

fl<l6  rpOMAT  (2HIi,  I*S,3H  t»,l«>,3H  P«,I5,3w  *<»,I^) 
wRITF(PT«,«Q7)(t<,)((«),V(K),nfK),K»1,w) 
n<>7  rnMMAT(I*>,1X,?H*»,F/,b,3H  V«,F7,!),3H  0«,IM 

1H1  COMTIMUF 

CAI.L  FOlfV  (X,V,  I,  )»n,P,Mr,,TlL) 

TF  (OF  mJG,F'’.P)  GO  TO  302 
tt».ITF(PTR,P<}5)T,J,P»^< 

no?  FPHmat (2WI», I?,3M  3«,I?,3H  P«,IS,3H  N«,I?) 

wRTTKPTP,nQ4)fK,xr«)»V(K),0(tO,X»1,M) 

H9<1  FPHmaT  (T?,  i X,?mx»,F  7 3m  y»,F7.?,3w  D«,I?) 
J/2  COFTI‘i|if 

C********** 

IF  (‘•CPOF  ,F0,  n GP  TO  330 
IF  (TLl.Mf.P.P)  GO  TP  33« 
wRTTF  (PTR,93b1 

935  FPPmat  (i7MZFMP  (FwGTM  chaTm) 
on  TP  339 
3JH  “0««0«TLI 
339  CPSTTMjf 

r********** 

T 1 • 1 

34iA  lrm,GT,“|Pn  GP  TO  35t« 

IF  f "P.l  F ,P»  f t n ) r.P  TO  3?P 

T 1 ■ M *i 
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GP  TO  34^ 

11  ■ HLPl 
II  ■ 11-1 

J1  • H 

L • HfP) 

IF  fpF0UG,fO,P)GOTO3/l 
wPITE(PTR,P03)L,P 
FORMAT(2Ml  ■,  T*),3M  P«,I5) 

CPMTINUf 
MAyML«0f n 

ir  fnm  .NE.otJi ) MAXML»MTNMfn(p,ncjn 
IFfll.GT.HAXML)  I1»fi*XHL 
0(P)  ■ 11 
r>EV(Pl  ■ Mn 
J1  ■ .Il*1 

IFO  .EG.-DGO  to  3«k) 

P » I 
GO  TO  37P 

TNPf  ii*niTF)Pm*n*.ii 
IF  (Jl.CO.n  GO  TO 
WRTTF(PTR,«>26)  Nr3,Jl 

^OWMAT(^«CWAIN,  1*5,  AHf  q PTS.HI 
GO  TP  ?qp 

r 1 NO  nr  cnain 
r 

33'>1  wPTTF(PTH,93P)NCS,LSNl,LSN?,Smi»CF  ,SSN,Nl  SN,N, 

1 (TNP(n,i»i,MLPi) 

9J«  rn9MAT(4Mr.HN«,I3,nf»T3,I4,T?,TR,TO,JH  NB,I3,?|Ih) 
no  3P'<i  i«i,N 
n«P  ■ ORPfA 

IF fOHP.LT.ORS)  GO  TO  4HM 
WRITE  f TaPEOT  OmiFFR 
ORP  ■ 1 

nHijrrpiPRP)  ■ »fTi 
OBitFFR(nBP4lT  a Y(I) 

PHliFrRfORP^?)  a OFVd) 

OP'IFTI  (ORP»3)  a 0(1) 

39P  CONTINUE 

OPP  a 0RP*4 

If  (OOP.l.T.nRSI  GO  TO  41P 
wRTTF(TAPrn)  OnUFFP 
npp  a 1 

4)a  OR'iFf  R(OnPT  a For 
nHMFFT(PRP*n  a P 

OHiiFF  T (0BP*2)  a p 
OMiirr  T (pqp»3l  a o 

r 

on  42P  Tal.NlPI 
4?P  NP(n  » ‘IPm  ♦ THJP(T) 

Trn  f“lP1  1 a TNPfMI  Pj  ) 

I a »'l 


3*oP 

3F>P 

t 

370 


BP3 

371 


C 

3R0 


926 

r. 
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433 

irn.rQ.-n  go  to  443 

TCCUn  • TNP(n*TCCL(T*l) 
T • I-t 
on  TO  4351 

44;) 

no  45M 

irfTcr,L(n.GT.cLCfT)) 

rtcfn-TccLd) 

TrfTca(T).LT.CSC{T)> 

card)  ■ Tfci  d) 

4b»* 

coMTiMur 

GO  TO  160 
ENO 

1 

.0 

0030 

E 

MU 

1026 

1 

MLPI 

1024 

C 

nrouG 

1B2E 

r 

I 

1032 

f 

j 

1036 

r 

p 

1034 

r 

K 

103E 

c 

1. 

1042 

f 

n 

1 046 

r 

11 

104A 

r. 

12 

104E 

f 

13 

1052 

{ 

;olMpw 

1056 

1 

Hw 

1054 

r 

pwF 

lOAf 

f 

f 4CT0R 

lC<-'2 

f 

MCW 

1CH6 

f 

McnoF 

1 CM4 

c 

OTR 

1C3E 

r. 

f TV 

tCi  2 

r 

TAPE  I 

1 r 1 6 

i 

TAPf  0 

ici  4 

t 

MCS 

IClE 

f 

MPS 

1C22 

f 

^(’1  r. 

1C26 

r 

►JPSC 

1C24 

f 

H 

1U2C 

[ 

>) 

ir.12 

L 

X 

1C  36 

r 

V 

2006 

L 

",f  V 

3536 

f 

P 

41»6 

f. 

NP 

4F  36 

r 

r.  LC 

4FR4 

1 

f.sc 

4FPE 

r 

C'.p 

4F  32 

f 

Tcn. 

4T  R6 

f 

TnP 

41  PA 

r 

M4»M| 

502E 

f 

IhHT  r P 

5P  12 

r 

’pnrr I 

12 

» 

iHHFf  P 

0032 

r 

nniirr  I 

50  12 

110 


( 

• J 

0000 

len 

088* 

161 

08C* 

201 

0C*2 

162 

098* 

163 

0*2C 

164 

MACE 

170 

0848 

160 

0PnC 

L 

MOO 

0000 

917 

0RtC 

200 

0C7* 

190 

0C82 

185 

MCI  8 

220 

1006 

918 

0CC4 

230 

0032 

1 

roF 

00M0 

240 

0OOf 

250 

0E36 

919 

0E4* 

260 

0E9C 

920 

0F72 

921 

0*HE 

922 

07  r * 

923 

1036 

924 

1078 

92S 

10H4 

1 

HI  w 

0000 

1 

.s 

0000 

?70 

1 1 22 

280 

1194 

290 

1 10E 

300 

1286 

310 

1 1*C 

330 

176* 

320 

125* 

301 

13*8 

896 

12CC 

897 

137f. 

1 

COI  1 V 

00M0 

102 

1470 

895 

1404 

;j94 

1 486 

339 

1542 

138 

15  36 

935 

1518 

340 

154* 

l*>0 

|58* 

360 

1592 

370 

15*8 

171 

1 5«6 

893 

\51  0 

1 

'MH0 

08'^M 

Ill 


IBP 

8?  3? 

OBP 

8236 

IR3 

823* 

OHS 

B23F 

er. 

82*2 

8246 

Q 

824* 

r»c 

824E 

HO 

8252 

TLL 

8256 

1 RMl 

825* 

1.5N2 

825E 

SOUPCF 

8262 

3SN 

8266 

NLSM 

826* 

90P 

PP38 

•I 

M8i4P 

899 

PP64 

898 

PPA2 

9J1 

P8HE 

•> 

P«riC 

9i? 

PPF8 

933 

PI  32 

934 

P18C 

tit 

P1H6 

9(11 

Pl(.* 

9(1? 

P218 

9(13 

P26  4 

?ll 

(1288 

904 

(129C 

908 

P2C  F 

9(16 

22Fr 

9(17 

P33B 

3(1 

432* 

.1(1 

H4HF 

5(1 

P5P8 

9(18 

H38F 

9(19 

P398 

/(I 

P6(12 

(ill 

P3F  * 

Olfl 

P44F 

91  1 

P45* 

H(1 

P4/f 

1 (IP 

pB;iC 

{'8?P 

9(1 

P57B 

9|3 

M*)04 

1 IP 

P6l  P 

12’» 

f»81  2 

914 

(16  9F 

916 

461  6 

1 3P 

IP 

918 

•i/(  2 

1 4'» 

.'Ml  2 

r 

X 

007A 

r 

V 

002C 

r 

t 

O02f 

f 

J 

onje 

r 

1. 

0032 

r 

M 

0034 

r 

>40 

0036 

f 

It 

0038 

i 

»<Drv 

03CC 

r 

XBAff 

0300 

r 

VPAR 

0304 

f 

TO 

0300 

f 

03OC 

L 

XNJ 

03E0 

E 

VNJ 

03E4 

L 

TfMP 

03E0 

K 

03f.t 

K2 

03F0 

?I)P 

0344 

N 

040C 

?5 

0220 

0274 

1 

«es 

0000 

i&p 

031  A 

50 

0200 

iUP 

0302 

750 

0300 

1 

.* 

0BitO 

^HOCrtAHU: 


U6Ct 

07At 

.P 

J336 

.'2 

D3E6 

.0 

D*<2£ 

. MLo 

D9AC 

. U 

09DA 

.V 

U9E8 

CIA14 

i'iOO 

EA4A 

A«S 

EA6  6 

. COM? 

EA3C 

• 1 1 ARC 

UAC8 

• RARtj 

£AFA 

EOF 

EUl  A 

»E\; 

EIJ3A 

MI'.IO 

t B4C 

.1 

EB5C 

51 

EU7  2 

52 

EDDO 

. S 

LB£4 

. A 

EC96 

ALOG 

C09E 

EXP 

EEA4 

AIXV 

CF16 

LF8A 

S6 

£F5d 

S3 

CF8  2 

.5 

CF8  6 

. ZERO 

IT  <Y- 
L>2D<> 

POIUTJ : 
CDLEV 

OoCE 

<*J 

D7AE 

.P 

0836 

.0 

D3E6 

. 0 

D92E 

• MES 

09U8 

.U 

D9DE 

• V 

U9EJ 

EAl  4 

MJO 

EA4A 

ABB 

CA6  6 

. COMP 

1.A8C 

.IIARO 

EACH 

. RARti 

EAFA 

EOF 

EDIA 

REW 

EB3A 

Mi:<o 

EU4C 

. 1 

EB5C 

SI 

EB72 

S2 

LUBU 

• 9 

EHE4 

• A 

EC9  6 

A LOG 

ED9E 

EXP 

LEA4 

AIMT 

EK16 

$6 

EF58 

S8 

EFS2 

.5 

CK86  .ZIjRU 


>•  Jii  b 


V.  ,i  DCK  I I4LI) : 


113 


Subroutine  CDLEV 


Purpose 

CDLEV  calculates  the  deviation  of  points  from  a trend  line  and 
returns  pointer  (a)  to  the  point  (s)  with  maximum  deviation. 

Procedure  Description 

The  deviation  of  a point  from  its  associated  trend  line  is 
calculated  as  follows  (refer  to  Figure  19  (a)). 

let  XBAR  - X(J)  - X(I) 

YEAR  - Y(J)  - Y(I) 

, .r  » year 

slope  of  line  % - m^- 

equation  of  llnei:Y  - Y(I)  • mg (X.  - X(I)) 

0 - mj,X  - m£(X(I)  - Y + Y(I)) 

- YEAR  X - YEAR  X(I)  - XEAR  Y + XBAR  Y(I) 

- YEAR  X - XBAR  Y - Y(J)  X(I)  + X^J)  Y(I) 
let  Z - -X(I)  Y(J)  + X(J)  Y(I) 

and  I is  given  by 

YEAR  X - XEAR  Y + Z - 0 

The  distance  d,  of  a point  (X(N),  Y(N))  from  line  i is  given  by 
d - YEAR  (X(N))  - XEAR  (Y(N))  Z 
V YEAr2  + XBAR^ 

The  deviations  of  data  points  from  trend  lines  are  all  cal- 
culated in  this  manner  with  the  exception  of  those  cases  in  which  the 
endpoints  of  the  trend  line  are  coincident  (l.e..  Island  chains).  In 
such  cases  the  deviation  is  measured  as  the  distance  of  the  point  from 
the  endpoints  of  the  chain  referring  to  Figure  19  (b) . 
if  CkX  - X(N)  - X(J) 

^Y  Y(N)  - Y(J) 

d ^X^  + 
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iwbrowtliw 


p 


SlIBWOUTINf  cnLEV(X,Y,I,J,L,P,HD,TU 

RFAi  Hn,TL»MOCv,¥n),vm 

REAL  XRAR»VBAR»TD»Z 

HEAL  *NJ,VNJ 

INTEGER  P,L(n,TFMP»I,J 

INTFGr.R  N,K2 

XHARbXM)-X(J) 

VRAR  «Vn)-Y(J) 
TL"XBAR*XRAR«YBAR*vBAR 
ZBX(I)*Yf J)-X( J)*Yf I) 

K«l4l 

K2«J-1 

MOFVbM.P 

Pb-1 

no  2PC<1  HaK,K2 

r 

r t«CCK  IF  chain  forms  an  island 
c 

If (TL.Nr.R.P)  GO  TO  25 
XNJsXrN).X(J) 

YNJ«Y(N)-Yf J) 

TO«XN.J*XNJ»YNJ*YNJ 
GO  TO  .5P 

2i  TD«AB5(XfH)*YBAR-YfH)*XBAR*ZI 

I'T  irfTD.LT.MnEV)  GO  TO  2BP 

TF  (TO.FO.mdCV)  go  To  IRP 
C 

c max  deviation 

c 

50  TFTP.EQ.-l'  GO  TO  100 

TFMP«L(P) 

I (P)»-l 
PbTEMP 
GO  TO  50 
100  MDFVbTD 

PbN 

GO  TO  200 
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c 

C TDbMOEV 

r 

15h  L(N)»P 
PaN 

2»P  COMTINUF 

ir  (TL.nc’.P.P)  go  to  28P 
MD«MnEV**fl,5 

return 

25U  T|.«(TI.)**P,S 

HO«Hnr;v/TL 

RETURN 

FNn 


K 

CPI.FV 

Pf1?4 

P 

( OLEIV 

M3C4 

1 

t 

.P 

0000 

A 

iftri 

A 

026 

I74A 

A 

93P 

1/FC 

A 

JO0 

tflF8 

A 

^(10 

1864 

A 

4i0 

1948 

A 

42P 

19H2 

A 

4J0 

lAlC 

A 

44P 

lA7fl 

A 

4b0 

IBPC 

1 

.V 

00(40 
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EXCLUD 


Purpose 

EXCLUD  deletes  entire  chains  from  the  data  base.  It  is  used 
primarily  for  removing  Islands  which  are  of  Insignificant  size  for 
the  scale  of  the  display. 

Procedure  Description 

After  an  array  of  numbers  Identifying  those  chains  to  be  deleted 
In  order  of  their  appearance  In  the  data  base  has  been  read  In,  the 
search  for  the  chains  begins.  The  data  Is  read  Into  core  In  6400 
byte  blocks  and  Is  examined  In  4-word  sections.  Each  section  Is 
Identified  as  a header,  a data  point  entry  or  end-of-chain.  If  It 
Is  a header  the  chain  nuaber  is  checked  against  the  first  chain 
nimber  In  the  array  which  has  not  been  found  In  the  data  base  up  to 
that  time.  If  a match  Is  not  made  the  chain  Is  copied  to  the  output 
buffer.  The  output  buffer  Is  the  same  size  as  the  Input  buffer  and 
so  the  data  remains  blocked  at  1600  words. 

Conznon  Blocks 

None 

Subroutines 

None 


119 


iA«4mt 


Flgurt  20  EXCLUO  FLOWCHART 


120 


OIMfNSlnN  Iff  (16910),  lH(2),LINK(l(*0),nB(  1009)  » 100(21  »TSAV(i) 

tEO'iTVALFNCF  (H (11,10(1)1,  (00(1),  100(1)1 

<t«TTF(5,9Q01 

0FA0(6,9«21  IN 

Wi9TTF(5,997) 

0FAn(5,9P'71  TOUT 

WRTTF(S,QQ91 

L»1 

1 MF4n(S,9Cll  1 I INK(L) 

u a iNK(i,  1 .ro.o)  CO  TO  10 
L»L*1 

Cn  TO  1 
JP  I.*! 

NL«LTNK(n 

ThFAObI 
2?  RFADdN)  0 

0(>  100  M«l,1600,4 
r,n  TO  (201, ?02, 203, 2041, IMtAO 
pffi  rF(TO(Mr,GT,|000)  GO  ro  50« 

AD»? 

TSAVf J )»T0(M) 

ISAV(21*I0(M*ll 

TSAV(.J1ii?R(M*2) 

TSAV(/O*T0(M*31 
GO  TO  100 
?02  IMFAObA 

IF(H|  .n?,T0(N411)  GO  TP  205 
IMF  Ar)*3 

TORdlBTSiVd) 
in0(r+1)«I8AV(21 
IPR( J*2)«ISAV(31 
IPP(T*31»T8AV(41 
T«T44 

TFd.UT, 15901  GO  TO  101 
I»1 

wRITEdnilTI  OR 

101  IO0(T1«T0(N1 
IO0(I*l)«TR(M*n 

102  tO0d*?)«T0(N»21 
(0Hd*31«in(M^3) 

T«I*4 

tFd,l.T,1590)  CO  TO  100 
I«1 

wHTTEdOIIT)  OR 
CO  TP  100 
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?{>.!  riRfT)«R(M) 

IFfBfM)  ,ro,lCH00,)  THEARbI 

no  TO 
?0*S  lBl4l 

Ml  bL  INXfU 
GO  TO  IBB 

?04  TF(R(M),m.1BBB.)  THFAObI 
IfIB  rONTTMlJF 

no  rn  ?B 

bf'B  ORfnBRfM) 

WHITEftnilTl  OB 
STOP 

001  KIWMATfTT) 

on?  roooATn?) 

907  FOHmaT  f3BHf  niter  OUTPUT  device  KiUM0FR(NN)) 

000  rORN!AT(29MFNTF9  T^PUT  DtVICF  N"IM0ER(NN)) 

000  r ORHA  1 ribWFWTtR  CHAlMS  TO  0F.  EL  INIINATE  D (NNN)  , / , 
1 24HEMTER  71R0  *S  L*8T  CHAIM) 
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ISLAND 


Purpose 

ISLAND  performs  two  calculations  on  the  chains  In  a data  base 
which  have  coincident  endpoints.  It  calculates  the  horizontal  and 
vertical  extent  of  the  chains  and  determines  their  separation  from 
nearby  ones. 

Procedure  Description 

The  calculations  are  performed  In  two  passes  of  the  data  base. 

The  data  Is  handled  In  the  same  manner  for  both  passes.  Data  Is  read 
in  6400-byte  blocks  and  4-word  sections  are  examined  and  Identified 
as  header,  data  point  or  end  of  chain  entry.  In  the  first  pass  an  array 
Is  created  containing  the  minimum  and  maximum  x and  y coordinates  of 
each  chain  with  coincident  endpoints  (l.e.,  4 pairs  of  coordinates/ 
chain).  From  these  figures  the  horizontal  and  vertical  extent  of  the 
Island  are  computed.  In  the  second  pass  the  coordinates  of  each 
point  In  the  data  base  are  subjected  to  the  following  tests  to  de- 
termine chain  separation  (refer  to  Figure  21): 

(1)  Vertical  Test 

The  y-coordlnate  Is  tested  for  Its  proximity  to  each  Island; 
that  Is,  Is  YMIN>y>YMAX? 

For  each  Island  In  the  y-range  the  minimum  x-coordlnate  and 
Its  y-value  are  compared  to  the  coordinates  of  the  point.  If  the 
x-coordlnate  Is  within  a certain  distance  of  the  x-mlnlmum  of  the 
Island  and  similarly  for  the  y-coordlnate,  a message  Is  sent  to  the 
printer. 

The  mlnlmi^m  separation  In  the  horizontal  direction  Is  4/256 
and  in  the  vertical  direction,  4/240.  These  values  were  empirically 
chosen  with  respect  to  the  resolution  of  the  display  medium  as  the 
minimum  distance  required  for  distinct  representation  of  two  lines. 
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Figure  21.  Island  Analysis 
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Figure  21  ISLAND  ANALYSIS 
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Figurt  22:  ISLAND  FLOWCHART 


Thus  If 

|XMIN  - x|  < 4/256  and 

I^XMIN  “ y| ^ 4/240  a oeasage  is  printed 

The  aaae  test  is  performed  with  the  island  z-maximun  and  its 
y- value . 

That  is 

IXMAX  - x|  < 4/256  and 
fYXMAX'-  y|  i 4/240 

(2)  Horizontal  Test 

The  x-coordlnate  is  also  tested  for  its  proximity  to  each 
Island,  that  is,  is  XMIN^x^XMAX? 

If  so,  the  minimum  y-coordlnate  and  its  corresponding  x-value 
are  compared  to  the  coordinates  of  the  point. 

If 

|YM1N  - y|s  4/240  and 

l*YMIN  ■ *1^  4/256  a message  is  sent  to  the  printer. 

The  maximum  y and  its  x-value  are  also  checked; 

|YMAX  - y|  < 4/240  and 

IXYMAX  “ *1  ^ 4/256 

For  each  point  fulfilling  either  set  of  conditions,  the  chain 
numbers  of  the  Island  and  the  chain  to  which  the  point  belongs,  the 
coordinates  of  the  island  and  nearby  chain  points,  and  vertical  and/ 
or  horizontal  separatlon(s)  are  listed. 

Common  Blocks 
None 

Subroutines 

None 
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OIMf»l3lO*J  an6»W),lB(160») 

OI**f*<SIOw  TFST(8?0)/ITEST(820) 
f8UIVALCNCCMB(l),8(l>) 
fOllIVALENCr  (TEST(»)»ITE8Tm) 

WHITC<3,998) 
iMf AO«t 
L*8??4 
K«1 

lOir  n(AO(6i  B 

PO  1 I»l,1688»4 
GO  TO  (18,?8,3«*r48)»lHEAD 

If  LW1«I8(I) 

IF(lNi,GT,lfff)  CO  To  2«f 
LM2»1B(I*1) 

IORlC«lBf 1*2) 

IMEAP«2 
GO  TO  1 

2f  1SE0«IB(I*1) 

IDSEO«IB(I) 

IHEAP»3 
GO  TO  I 

3f  X9TART«B(n 
VSTART«8(I*n 
XMAXbXSTART 
XMJNsXSTART 
X0VMAX"X»TART 
XpyMl^BXSTART 
V***X«VST  ART 
YMTNbYSTART 
Y0X‘<AX«Y9TART 
YOXMJi^svBTART 
IHF  AOaA 
CO  TO  t 

4f  TFfBfn.EQ.lfff.)  GO  TO  SR 
XENO«B(T1 
YE*<0«B(l*n 

XNAXaAMAXl  (XMAX.xr.NO) 

IF (XM*X,E O.XENO)  YnXMAXaYrwn 
YMAXbAMAXI (YMAX, YEND) 

TFTYMaX.F Q.YENP)  XOTMAXaXENn 
XMTN«AM|M1 (XMIN,XFM0) 
tr f xMiN.rn.xcNO)  YOXMiKaVENO 
YMTN«A*^IN1  (YMIN,YC^0) 

IF(YMl»*,ro,YE^O)  XOYmIbbxTMO 
GO  TO  t 

5R  Tr((XSTART.EQ,XFMn).ANO,(Y8TART,CQ.YtiYO))  L«lfF84 

»<RTTF(3,999)  L^t ,Ln2, I0RIC,IO8f0, T8rO,XSTA»T»Y8TART»XtMO,YFNO,U 

TF <L.CQ,B224)  GO  TO  Bf 

t»RT.tF(3,997T  XmIm,yOxMIm 

wRTTFrJ,Q96)  X*<AX,YOxHax 

w»TTE(3,99hT  Y“IN,xUyOIn 

WBITF(3,994)  Y“AX,xOYM4X 

OFl TAX«XM4X-XMTN 

OFl  TAY«Y“4X-Y**TN 
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WRITF(3,993)  OCLTiX»OEl TAY 
1*8224 

ITE8TtK)«lStO 

TrST(K*n*  XRAX 

Tt8T(K*2)«Y0*M4X 

TE8T(X*3)*  XMIN 

TFSTCK*4)*Y0XM1N 

TE8T(K*5)«  YMAX 

TE8T(K«6)*X0YHAX 

TE8T(K*F1«YMIN 

TE8T(K*8)«X0YMIN 

K*K^9 

68  IHEAD*! 

1 CONTINUE 
CO  TO  10R 
C 

C AT  THIS  POINT  WE  H«Vf  «N  ARRAYr  TEST  (ITEST).  CONTAINING 

C the  min  and  max  X AND  Y COORDINATES  OF  ALL  ISLANDS  IN  THE  DATA  RASE 

C 

218  CALL.  REWfAT 
IHEADsl 
C0N8T*4,/2A8, 

C0NSTX>4./2S6, 

WRITE(3,990) 

380  HEA0(6)  R 

DO  2 Iai»ie08»4 

GO  TO  (31P.320,330*A40)r  IHEAD 
310  IFfIB(I),GT,18P0)  GO  TO  400 
IHFAO  •? 

GO  TO  2 

320  ISEQ*IRn*l) 

IHEAP*3 
GO  TO  2 
330  IHFAD*4 

GO  TO  340 

340  IFfBm.NE. 1000.0)  UO  TO  34? 
iHEAOaf 
r.O  TO  2 
34?  XaR(I) 

Y«RIl*n 

C 

C TE3TX 
C 

DO  347  K«6. 630*9 
irCITCSTfW-RJ.FO.lSLO)  CO  TO  347 
TrtCTCSTeK),LT.Y).0R,{TE8T(XA2),6T,Y))  60  TO  344 
TF((rTF8T(»<-3)aCONST),LT,Y),OB,r(TF8T(K-3)*CON8T),6T,Y)) 
cr.o  TO  343 

XDlFFaARS(TESTrX-4)-X) 

IF  (XOIFr,6T,0,01)  60  TO  343 

WRITE  (3,99l)ITFST(R-5),TE8TfX-4),TE8T(X-3),18F0,X,Y,X0IFF 
343  lF(((TEST(K-l)*CONSl),LT.Y).0R,(fTE8T(K-n-C0NST),GT,Y)) 
cr.o  TO  344 

XDTFraAnS(TF8T(X.2)*X) 
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IF  (XO1FF,GT,8I.01»  eo  TO  344 

WRITE (3, 901 )ITE8T(K-S), TEST (K.2), test (K-l),I8t0»X,r,K01FF 
C 

C TEST  y 
C 

344  IFt(TCSTfK-4),LT.X).0R,(TCST(K-2J,CT,X))G0  TO  347 

IFf  ((TE8T(X*n*C0NSTX),LT,X)  .OR.  f (TEST  (K«l  )»CONSTX)  ,6T.  X)  ) 

IGO  TO  345 

y0tFF**BS(TE3T(K)*Y) 

ir(yOtFF',GT.0,*n  go  to  34S 

wRITF(3,992)ITEST(K-8),TFST(K*n,TpST(K),ISEO,X,y,yoiFF 

345  lF(((TE8T(K*3)*C0NSTX),LT,X).0R,((TrST(K*3)*C0MSTX),GT,X)) 

ICO  TO  347 

yOIFF«AISS(TEST(K*2)-ty) 

IFfyniFF.GT.B.Bn  GO  TO  347 

WRITE (3, 992) 1TF8T(K-5), TEST (K*3J,TF8T(K42),18E0»X,V,V1)IFF 
347  CONTINUE 
2 CONTINUE 
GO  TO  J0R 
400  STOH 

•90  F0RNAT(6HI8,1.8N,2X,8H1S. (X,y),8X,7MM,L,L8N,2X,9HM,L,CX,Y),3X, 

C7M0ELTA  X,3X,7H0ELTA  Y) 

991  FORMAT (2X,I3»2F 8.4, 3X,I3>3X,3Fn. 4) 

092  FnRMAT(2Y,I3,2F8,4,3X,I3,3X,2F8.4,9X,E»,4) 

993  F0RMAT(6X,  1 |H  OEITA  X ■ ,F8,4,HH  OEtTA  Y • ,F8,4//) 

994  F0RMATC6X, 1RM  MAX  Y IS  ,F8,4,8M  AT  X • ,F8,4) 

995  F0RMAT(6X, IRH  mjn  Y IS  ,Ffl.4,8H  AT  X ■ ,F8,4) 

996  F0RMAT(6X,SRH  MAX  X IS  ,FB,4,«M  AT  Y ■ ,F8,4) 

997  rOPMATtSX, IBM  MIN  X IS  ,FB,4,8H  AT  Y • ,FB,4) 

99B  F0RMAT(3X,BH0L0  in,5X,6HSnURCE,3X, 1 IMORIGIN  SE<3.,3X, 

1 IMMTOTAt  St0,,5X,6HXSTART,5X,6MYSTART,6X,4HXEND,7X,4MYEN0) 
999  F0RMAT(3X,T3,I4,3X,I5,6X,I5,8X,I«,5X,  4 ( 3X,F8.4 ) , 2X,  A 1 ) 
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lORlG 

3306 
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33RA 
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996 

0C96 
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0C6A 

994 

HC3E 

t>f  LTAX 

3306 
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330A 

993 

OCOC 

RCN 
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const 

33C6 

CONSTX 

33F2 
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prograhsi 

6416 

• J 

84F6  ,U 

6924  .V 

6932 

• I 

95SE 

AMAXl 

98FB  .2 

9962  61 

9996 

63 

9SC* 

ABS 

9SE6  ,CO«R 

966C  .RARG 

963E 

REK 

96SE 

,s 

9692  ,MFS 

9719  AMim 

9722 

lntry- 

6416 

POINTS! 

• J 

8902  |il' 

8926 

.V 

8932 

PI 

999E 

AMAXl 

9970  ,2 

9982 

91 

9998 

S3 

99CA 

APS 

99E6  .CONP 

960C 

.RARG 

963C 

REW 

969E 

.8 

9692  ,M!3 

9710 

AMINI 

COMMON-BIOCKSI 
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ENDPT 


Purpose 

ENDPT  lists  the  chain  numbers  of  the  chains  in  a data  base  and 
the  coordinates  of  the  first  and  last  point  In  each.  It  Is  the  most 
general  tool  In  the  package  for  checking  data  Integrity. 

Procedure  Description 

Data  Is  read  Into  core  In  6400—byte  blocks.  Four-^ord  segments 
are  Identified  as  header,  data  point  or  end  of  chain  entries.  If  a 
header  is  being  processed,  the  WDBl  line  segment  number,  chain  number 
and  type  code  are  saved.  While  examining  the  data  points  the  coor- 
dinates of  the  first  and  last  are  saved.  When  an  end  of  chain  Is 
detected  the  stored  Information  about  the  chain  Is  printed. 

Common  Blocks 

None 

Subroutines 

None 
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Figur#  23  ENOPT  FLOWCHART 
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OIMfNSlON 

FOUlVALCNCt 

hPITF(3,9Q0) 

IHEAO*l 

La6??4 

Pal 

100  RF*D(6)  8 

DO  1 Ial,t600»4 

CO  TO  n0,20*3et40)f  IHEAt) 

10  lA'lalBd) 

TF(L*<1  .GT.1000)  GO  TO  200 
LN2aIB(l*l) 

I0RIGalB(I«2) 

IHEA0a2 
GO  TO  1 

20  ISEOaIBn«n 
lOSEQalBd) 
lHE*Oa3 
GO  TO  1 

30  XSTABTaRd) 

VSTARTaHd^l) 

THfA0a4 
GO  TO  I 

40  TrTRdl.FG.l*****.)  '■,0  Tn  00 
XFSOaBd) 

TF*<Oaad*n 
CO  TO  t 

t>0  ir ((X3T4RT.FO,xrMn).*ND,Cy3TART,F9,YENO))  Lal0704 

WRTTE(3,999)  1^1 ,Ln2, lORIG, inSFO, I8C9» XSTAPT » VST ART, XtNO» YFMO,L 

IFfL,CO,0?24)  60  TO  00 

La8224 

60  iHFAOal 
1 C0*JTIMIJF 
on  TO  100 

200  -RITCf5,201) 

201  F0HMAT(5M.  STOP) 

990  FORMATflX,6H0Ln  TO, bX,6MSnUHCF . 3X, 1 1H0RI6IM  SE0,,3X, 

1 10MTPTAL  SC0.,9X,6HXSTART,5X,6HVSTART,6X,4HXFND,7X,4MVENn) 

999  rnHMAT(3X,T3,l4,3X,l9,6X,19,0V,in,5X,  4 (3X,F0.4)  ,2X,  A1  ) 
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PRINTM 


Purpose 

PRINTM  produces  a map  of  a data  base  at  a user-specified  level 
of  detail  on  the  printer.  It  Is  a crude  tool  for  checking  the  con- 
tents of  the  data  base. 

Procedure  Description 

Through  a subroutine  call  to  INITMP  a 60  x 100  matrix  Is  Initial 
Ized  with  blanks.  Data  Is  processed  In  6400-byte  blocks  and  4-word 
sections  are  Identified  as  header,  data  point  or  end-of-chaln  entries 
If  a data  point  of  rank  greater  than  or  equal  to  the  user-specified 
minimum  Is  being  examined,  the  coordinates  of  the  point  undergo  a 
linear  transformation  (see  Figure  24)  producing  a point  in  the  range 
of  the  matrix  indices.  A zeal  to  Integer  conversion  is  performed  and 
an  'X'  Is  placed  In  the  cell  of  the  matrix  corresponding  to  the 
transformed  and  converted  point  coordinates.  After  all  the  data  has 
been  processed  the  matrix  Is  printed. 

Common  Blocks 

None 

Subroutines 

INITMP  - Initializes  a 60  x 100  matrix 


CALL  INITMP 

(MAP,  LX,  LY) 

MAP; 

60  X 100 

LX: 

defined  as  'X' 

In 

the 

subroutine 

LY: 

defined  as 

in 

the 

subroutine 
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DIMENSION  I8UFF(160tf) 

FQUIVAUFNCr.  (M'lrFfDrlRIJFFd)) 
integer  D,DD 

CALL  INITHP(MAP,LX»LV) 

WRITF (5,9QP) 

REA0(5»9<»8)  00 
9R  N«1 

m f?FAD(6)  BUFF 

DO  IBM  M«ld6BBrA 
IF(N,CQ,0)  60  TO  2B 
IFfN.EO,?)  GO  TO  3B 
IF(IBurF(MI,6T,l0Bfl)  GO  TO  2B0 

N>2 

GO  TO  30B 
3B  N«M 

GO  TO  3BB 
20  X»RUFF(M) 

Y«RUFF(M*n 

0«IBUFF(M4S) 

IF(X,f Q.1B0B#)  GO  TO  100 
IFCO.LT.DO)  GO  TO  300 
IX«IFTXr(X*,4)«147>.b) 
IV«60-IFIXr(T*,27)*00,*,5) 
TF(IX,6T,100)  IX-100 
IFflX.LT.n  IX«1 
IF(IY,GT.«0)  IY«60 
IFflY.LT.n  IY«1 
MAP(IX»IY)«LX 
GO  70  300 

100  IFfN.EO, n 60  TO  200 
N»1 

300  CONTINUE 
GO  TO  10 

200  DO  101  J«l,60 

NRITEf3,900)  (MAPfK»J),KBl,100) 

101  CONTINUE 
PAUSE  99 
LX«LY 

GO  TO  98 

900  FORMAT(100A1) 

950  FORMATflH  ,2110) 

990  r0RMATfI2) 

999  rORMAT(22HFNTER  DETAIL  LFVEL(NN)) 
END 
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UNDCFINCOI 

MONf 
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Subroutine  INITMP 


Purpoee 

INITMP  is  an  asseably  routine  that  Initializes  the  nap  matrix 
and  sets  variables  LX  to  'X*  and  LY  to  . These  variables  are 
used  by  PRINTM  to  Identify  those  cells  of  the  natrlx  corresponding 
to  the  coordinates  of  nap  data  points. 
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0O4AR  FND 

* INITMP  ndflSR 
LOOP  P01AR 

SAVR  P00PR 


0(500R 

entry 

INITMP 

0O0PR 

SAVR 

03 

8 

0O0RR 

O0C0 

0000R 

INITMp 

STM 

12. SAVR 

0O0CR 

OlEF 

0002 

LN 

14,2(16) 

0O1PR 

CRP0 

1770 

LHI 

13,6000 

001 4R 

CSC0 

2020 

LrtI 

12,C«  * 

oniRR 

40CC 

0000 

LOOP 

sth 

12,0(14) 

00  ICR 

CAE0 

0004 
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14,4 

0O20R 

CBO0 

0001 

SHI 

13,1 

0O24R 

4220 

001RR 

BP 

LOOP 
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LHI 
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0O2CR 

40CF 

0000 

STh'  ' 

12,0(16) 

0O30R 
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0000R 

LH 

12, SAVR 

0O34R 

C8C0 

RF20 

LHI 

12,C»0  1 

0O38R 

48EF 

0006 

LH 

14,6(15) 

0O3CR 

40CE 

0000 

sth 

12,0(14) 

0O40R 

O1C0 

0000R 
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12, SAVR 

0O44R 

4AFF 
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AH 

15,0(15) 

0O4RR 

0O4AR 

030F 

BR 

FNO 

15 
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CHAP 


Purpose 

CHAP  creates  an  Image  In  a drum  file  of  all  the  map  points  with 
tank  equal  to  or  greater  than  a user-specified  level  through  sub- 
routine calls  to  PALLET. 

Procedure  Description 

Data  Is  read  Into  core  In  6400-byte  blocks.  The  rank  of  each 
data  point  Is  checked  and  If  It  Is  above  the  user— specif led  mxnlmum, 
the  X-  and  y-  coordinates  of  the  point  are  stored  In  X and  Y arrays. 
Upon  completion  of  the  processing  of  a chain,  the  contents  of  these 
arrays  are  written  on  the  drum  file  as  an  Item  (l.e.,  line)  of  the 
map  Image.  If  a chain  has  more  than  200  elements  of  sufficiently 
high  rank  to  be  kept  In  the  Image,  more  than  one  line  Is  created  for 
that  chain. 

Common  Blocks 

Block  Name  Contents  Description  of  Contents 

RG  BUFF  Input  buffer 

Subroutines 

GET:  fills  Input  buffer 
no  parameters 
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dimension  X(8«R),V(S08),NAME(2) 

DIMENSION  IBUFE(16«0)#BUEFCt6t0) 

EQUIVALENCE  (IBUFE(t)»Burrri)) 

EXTEKNAL  WRKSB»SAVEA 
INTECER  SAVEA»EBBOB 
INTE6E0  D,DD 

C K«BFGINNIN6  or  COMMENT  STMING 
c J*LAST  entmy  of  CUMMENT  STMINC 
C M»  CUMMENT  BUFF  ENTMV 
C NUMa  numbem  or  subimaqeQ 
MMITE(S,1BI) 

MPA0(9»1BB)  NAME 
NMITE(9,99S) 

REA0(9»999)  D 

MaEBOB 

NUMaB 

IMEAOal 

lal 

CALL  SETSAV(WMKSP,9AVEA) 

CAIL  INTOMm(E9,1E8,Q>0,«,EMMOM) 

IF  (EMMOM.NE.B)  STOP  IB 
call  1NTFMP(29,1,SAVEA) 

CALL  DEFILFtMMK8P,B«MEA0FlLf »9AVEA,|) 

CALL  OPENfHMKSPf BHRLAOPILf  f 9AVEA) 

10  NAMEIPIaNAMt (2)«t 
NUMaNUMal 

CALL  0MENIfNAME,-,4,-,27,.2B,.41,B> 

Na4B 

16  MaH44 

IF(M,r.T,159B)  CALL  6lT 
IFdHCAO.FQ.B)  GO  TO  17 
IF  (TMC*0,tO,2)  GO  To  10 
IF(BUFF(H),EO,100B.)  60  TO  2B0 
lO  IHEADa2 
GO  TO  16 
IHEADaB 
GO  TO  16 

17  xmaBUFF(M) 
r(I)aBUFF(M*n 
Tm»,41B«Y(I) 

IF(X(I).E0.1BBP.)  r.O  TO  30 
00aIBUFr(M*3) 

IF(0,6T,0D)  GO  TO  16 
NaN«B 

IP(N,6T,17BB)  go  to  21 

lalal 

GO  TO  10 

21  IFd.GT.nCALL  LINF9(0mlINE  »X(1),Y(1)»I,7) 
CALL  CLOSEKNAmEI 

xmaxm 

vm«v(i) 

Ia2 

CO  TO  10 
30  Ia!-1 
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200 


201 


100 

101 

998 

999 


CALL  LtNtSrSHLTNC  . » ( 1 ) , V 1 1 ) , I , T ) 

IMCADal 

NsN*12 


!■! 

GO  TO  18 

CALL  CLOSFKNAMn 

NAMF(8)«igAMf  (2)-num*1 


f0..0..51t.,479.,0) 


CALL  OFFNlieHNAP 
DO  201  1«1,NUH 
CALL  INCLUD(NAME,-:i4fl,,«332,,6T8.,1100.,0,0,, 
NAME(2)«N4MC(21«1 

> 

*0.*0.*911.,479.,0) 
»0*»0«»011»»AT9,»0f0,»0 
) 


CALL  CLO9FKOHM40 
CALL  OPEMUSMXORLO 
CALL  INCLU0(0NM*p 
CALL  CLOSEKSHWOPLO 
STOP 

F0RM4T(2A4) 

FORMATdlM  enter  N*ME) 

F0RMAT(23N  CNTFR  OF  TAIL  LCVfLFNN) 
F0RMATFI2) 

END 
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1 
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17 
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01  Af. 
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.^R.rRHHAP  ) 
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RETURN 

END 
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8EBE 
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^TGHB 

ClCC 
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.8 

C24E 

•J 

C32E 

iP 

esaa 

.8 
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.0 

C4*f 

.HES 

C92C 

css* 

PI 

0966 

COS 

099* 

SIN 

0644 

.* 

0BP8 

^LOG 

OPPE 

EXP 

0904 

HIN0 

0916 

MAX0 

0928 

.1 

D9S6 

AMOO 

0966 

AINT 

090B 

ai 

D9EC 

• 2 

0A2C 

FLOAT 

0A6E 

IPIX 

o*ap 

.T 

0812 

.** 

OBP0 

.enpp 

0996 

.RRAPG 

OBDB 

96 

oct* 

.PAPG 

0C4C 

•a 

0CP6 

.9 

OCP* 

.ZERO 

OCPf 

iTAUg 

ooep 

atan 

OOP* 

CNTRVaPOINTSi 
69BA  GET 

69PE 

getcor 

6602 

GETLNC 

6604 

CORPTR 

6606 

GIVCOR 

660A 

GIVLNC 

660C 

8XVPTR 

660E 

PLT911 

6612 

PLT4P9 

6616 

FSPll 

661* 

P4P79 

661E 

OTOP 

6622 

N63 

6624 

N90 

66&B 

HESS 

6678 

LENGTH 

66PA 

SUPON 

66PF 

supopf 

6682 

ENHUX 

6684 

IHPTAB 

660e 

PINO 

6790 

driver 

6894 

PICTUR 

PIEE 

AMPCAP 

P238 

018 

736C 

OOTBPR 

7374 

UP0I9 

7904 

naming 

P662 

lOPPN 

768* 

OPfNl 

79P2 

CL08EI 

7744 

INCLUD 

P7BC 

BINO 

7920 

BLOCK 

78*0 

CHAR 

7948 

POINTS 

79BC 

L1NC8 

7C9* 

BUPPRP 

7C9C 

ATTRPL 

7C7* 

P0HI9 

7014 

REPL*C 

7089 

rgepp 

7092 

PHFAO 

7P9* 

9TPACT 

7E88 

8ET9AV 

7ECf 

chkBav 

7EEC 

HMEAD 

7P9B 

HPITEB 

a92A 

PFAOB 

8990 

SAVE* 

80*9 

NRK9P 

8138 

POLL 

6249 

PEC 

9494 

XHlT 

8936 

AMR 

8934 

TPAKUP 

9E78 

0I9PUT 

9020 

EP*8E 

9082 

NEH*HP 

914A 

UPHI8 

919C 

EPM9G 

9ica 

FIXREC 

99PC 

CONVT 

9688 

A0080N 

970E 

OPPILE 

9874 

oetoat 

9922 

pndson 

9A1C 

INTOPM 

9B8E 

INTFhr 

9C69 

AOK 

9C6* 

Onf 

9C6C 

8AVE1 

9C9C 

9AVE2 

9CAC 

saves 

9CCC 

SAVE4 

9CPC 

PMP09C 

90BC 

OPEN 

9002 

PUTOAT 

9ECE 

PEMINO 

9PE2 

ADBRO 

9F9A 

AONOO 

*164 

fchain 

A2F8 

9128 

A368 

gthain 

A3DF 

PMP0P3 

«3E2 

FMPGP2 

A3F6 

PMPOPl 

ASP* 

PMPGP8 

A418 

HASH 

*486 

NAHGEN 

*HH* 

NOOOK 

AB2C 

8AV0K 

*949 

PHRJRO 

*990 

GMINO 

*960 

NXTSON 

*62* 

PA080N 

A694 

PMPRO 

*P16 

savPag 

AP82 

ERROR 

A79E 

SEOSON 

*BP6 

PHPTRC 

*868 

FTOPF 

*89* 

FTON 

*9*8 

lOENT 

A926 

I0FMT2 

*980 

RfClVt 

*9A6 

LISTEN 

A9CC 

ONCIVE 

A9P2 

UNLI9N 

**CC 

AL0U2 

AC30 

RSHTCM 

AC90 

M9MTCM 

Acaa 

IMPINT 

*03* 

STHP70 

*E82 

BUPSIZ 

AF86 

BUPNUH 

AERA 

MAIT 

*P9C 

HAIT2 

0190 

PATCN4 

9198 

INTBIT 

BI9A 

P*TCN7 

9isr 

»P 

H178 

“PMATN 

930r 

SMITCM 

J 
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B3i* 

FRTRtG 
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TQhBnT 

B336 
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,B6IN0 

BD4A 

.BCIN 
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BE72 
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BE74 
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.exitn 
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BFFt 
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BFAF 
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BFCB 

TOFF 

BFD2 

RTRAC 

CR2E 

ORFF. 

CPF4 

TOUTl 

CllE 

T0UT2 

C198 

STQHB 

C1B6 

nSFCT 
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SIZBIIF 

C1D8 
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ClOA 

flUFlNO 
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MSGMAX 

Cl  or 

Busro 

cice 
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JOIN 


Purpose 

JOIN  connects  chains  which  have  coincident  endpoints.  Its  pri- 
mary use  is  to  connect  a sequence  of  short  chains  so  that  minimum 
chain  overhead  is  produced  by  CMAP  when  creating  the  map  Image  on 
drum. 

Procedure  Description 

A list  of  chains  to  be  connected  is  read  into  an  array  in  the 
order  in  which  they  are  to  be  connected.  The  data  is  then  processed 
in  6400-byte  blocks.  A block  is  searched  for  chain  headers  and  when 
one  is  found  it  is  compared  to  the  chain  number  array.  If  a match 
is  found  the  entire  chain  is  copied  onto  a temporary  file,  usually  a 
drum  file.  If  the  chain  is  not  to  be  connected  it  is  written  to  the 
output  buffer.  When  all  chains  have  been  examined,  the  connection 
process  begins.  The  temporary  file  is  repeatedly  rewound  and  searched 
for  each  of  the  chains  in  the  array.  As  a chain  is  located  its  points 
are  copied  to  the  output  buffer.  If  a chain  is  to  be  connected  in  the 
opposite  order  to  its  current  one  it  is  Inverted  first  and  then  con- 
nected. 

Conmon  Blocks 

None 

Subroutines 

None 
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DIhfwSTOW  INH(l6SI(»Wl0llTRf  16Jl^),T8ri6(K»ll,INT(2mfW) 

DIMENSION  tSAV(-0 

dimension  lfIN(2),KnOT(2),*S(2)»LlNK(20) 

FOllIVALENCMlNB(l),*INn)),(IPUTS(l),XOUTm),(ISm,XSm) 
EouivAiENcr  f iNB( n»iNT(m 
1 wBlTFf5,Q99> 

RFAO(5»930)  in 
MPm  (5,990) 

READ(5,90«)  IOUT 
mRTTE(5,997) 
nFAD(5,90O)  I3K 
I»1 

lMEAn«t 

TI»1 

WRITE (5,99ft) 

L*1 

3 nfAD(5,99)l)  LINK(t) 

if(link(l).eq,r)  go  to  1H 

L»L*1 

GO  TO  3 

)!» 

IFd.LF.Pt)  GO  TO  5RH 
i r?EA0(IN)  inB 
on  M,i,i6ffj»,4 

GO  TO  ( 1B1  , 1B2,  1B3, 1B4,  ISIS,  IBB)  , IHFAD 
101  IF(INB(M),GT,1B0^)  go  to  4B(1 

ISAVd  )«INH(M) 

TSAV(?)«IN0(M41  ) 

I8AV(3)«I*gR(M4?) 

ISAV(A)»INRfM*l) 

IHEAD-2 
CO  TO  ?0 
10?  DO  ASI  J«1,L 

IF(INH(M*n,EO.LlNK(J))  GO  TO  |r4 
IF(INP(HAn,FO.-LINft(J))  GO  TO  104 
40  CONTINUE 

iniiTBf  i)«rsAvM ) 
imiTHI 1*1 )«ISAV(?) 

I00TB(I*2)«I3AV(3) 

Tni)TB(U3)«ISAV(4) 

I ■ T«4 

IFd.l  T.1598)  go  to  11 
I»1 

WRITFdOUT)  inilTH 
31  lniiTH(I)«lNO(Ml 
IMFA0»3 

I0llTB(l4l  )»INH(M4l  ) 

21  IOUTH(Ia?)»INH(m*?i 
TOiitP(T*3)«INB(M«.3) 

I»l44 

IF(I,LT,159B)  go  to  ?11 
I • 1 

wHTTfdO'in  lOUTR 
r,o  TO  ?o 
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103  *0'lT(I)«XTNfM) 
X0UT(I*t)«XIN(M*|) 

IF(*lN(Ml,ro,l»U»^,)  THFAOat 
GO  TO  21 

104  ia(II)"ISAV(l) 
IS(IU1)«ISAV(2) 
1S(II*?)«1SAV<3) 
IS(II«3)«ISAV(4) 

rF(M,l  T.1S98)  GO  TO  107 
II-l 

W0ITEM3K)  IS 

105  CONTINUI 

107  1S(I1)«I»JB(*») 

IHEAObO 

I9(IT*l)»INB(M*i) 

100  IS(1U2)»INB(0*2) 
IS(1I*3)«INB(M*3) 

!I»II*4 

IF(I1,LT,159S)  GO  TO  20 
1I»1 

W»1TE(19K)  IS 
GO  TO  20 

106  X8(l  n««TNfM) 
X8(II*ll«XTNfM*n 
IF(XIN(mi .tO.lOOO.)  IHEAOal 
GO  TO  lOB 

20  CONTINUE 
GO  TO  2 

400  X8(m«1000, 

NBITFflSK)  19 
TOUTBdlaO 
I0UTB(T*1 )»0 
inilTB(U?)«4 
lOilTBf  1*3)»0 
I«I*4 

lFn,LT,l89B)  GO  TO  414 
1»1 

nBITF(TOUT)  TOuTB 
414  iniiTBm»i 
j»i  iN«(n 
urJ.iT.o)  j«-j 

TniiTM(i4j)«j 

I0UTB(I*?)»0 

T0UTB(U31«O 

I«I*4 

IFd.l  T.189B)  1.0  TO  404 
T»1 

NBtTFdOiiTl  TOIITO 
404  00  401  .I"1»L 
CALL  BCwdSM 
t7»0 

NL»LTn4c,T) 

TFfNL.GT.O)  GO  TO  443 
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IZ"1 
NL«>NL 
4R3  !HEAO«t 
II»1 

42P  PFAOdSK)  IS 

DO  492  Mat,l«C>F>»4 

r.o  rn  (6(ii,6fl2,6iA3,6H4,en5).iMEAn 

601  IME*0»2 

IFnS(H),GT,10»'»')  *»AUSE  100 

60  TO  40? 

602  IHEA0«3 

IF{I3(“An.NE,NU  r,0  TO  402 
IMFAD»4 

IFdZ.EO.l)  IHEAObS 
r.o  TO  40? 

603  !F(X8(»«).EO,10O0,)  iHEAOal  , 

GO  TO  40? 

604  TF(X3(O),EO,10O0.)  Go  to  401 
XOUTin^XSt*^) 
xouTd*n«xs(M*n 

I0UTBd4?)»I3(H^2) 

TnuTBd*3)»I3(N*3) 

1«I*4 

IFd.l  T.159S)  GO  TO  402 
I»1 

wRTTFdnuT)  TOUTS 
GO  TO  402 

605  irfXSfM),t0,>0OO,)  Go  TO  606 

XT^dn«XSfM) 

xiNdMn«xs(s*n 

INTdT*?T«ISfM4?) 

TNTd|431iIS(M43) 

IT»IT*4 
40?  rONTiNUF 
GO  TO  420 

606  TI»II-4 
TFdl.Lr.O)  GO  TO  401 

xouTd)«xiNf  in 
xoiiTinnaxiNdi^n 
I0UTHd4?l«IMTdI*?) 
Tni)TBd43)«INTdl43) 

T»T*4 

TFd,LT,l!>90)  GO  TO  606 
T«1 

wMTTFdOUT)  TOiiTB 
GO  TO  606 
401  CO'^TimiF 

XOUTdIalOOO, 

I«T64 

IFd,LT,1606)  GO  TA  601 
T»1 

*PTTF(tnun  TOUTS 
tot  XOliTT  daMOO, 

»P  r Tf  dOllT  1 TOi'TR 
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50(1  STOP 

900  F0RMAT(I2J 

901  F0PMAT(I4) 

996  FORMATOIMFNTCR  CHAINS  TO  BC  LINKFO(NNNN) , 
t 34HNCGATIVE  NUNBFR  FOR  OPPOSITE  ORDCR*/* 

1 40HFNTER  ZERO  AS  EHO  OF  CHAINS  TO  RE  LINKIO*/, 
1 18HNULL  CHAIN  TO  EXIT) 

997  rORMATOlHCNTER  SCRATCH  DEVICE  NUMRER(NN)) 

998  F0RHAT(3PHENTER  OUTPUT  DEVICE  NUHRFR(NN)) 

999  F0RHAT(29HENTFR  INPUT  DEVICE  NUMBER(NN)) 

FNO 


1 

.u 

00MH 

r 

INR 
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E 

INT 

0E5E 

E 

XIN 

0E5E 

E 

lOtITR 
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C 

KOUT 

2P9F 

t 

13 

469E 

E 

XS 

469F 

ISAV 

5F9E 
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5FAE 
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• 1 

0000 

900 

onjf 

IN 

5FFE 
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OEOC 

lOilT 

6002 
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ODE  4 

ISK 

6006 

1 

600A 

IHFAD 

6012 

II 

6016 

996 

0D4E 

L 

601  A 

3 

0OCe 

901 

0046 

10 

0)26 

500 

0038 

? 

0144 

1 

• J 

0000 

?0 

06(4 

M 

6022 

101 

0162 
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021F 

103 

043F 

104 

0460 

105 

0562 

106 

0652 

400 

060A 

4 0 

026A 

J 

602E 

31 

034f 

?1 

03A2 
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1B7 

0562 

ms 

0506 

414 

07A4 

404 

0872 

401 

0CA4 

L 

RCw 

0000 

IZ 

604F 

ML 

6052 

403 

08C4 

420 

0004 

402 

0080 

601 

0910 

602 

0940 

603 

098C 

604 

0904 

605 

0A0A 

1 

0000 

606 

0O9F 

501 

0O0A 
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0000 
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.V 
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EDITDB 


Purpose 

EDITB  Is  a point  editor.  The  routine  allows  addition  and  de- 
letion of  points  and  changes  In  the  coordinates  and/or  rank  of  points. 

Procedure  Description 

A list  of  chain  numbers  In  which  the  points  to  be  edited  are 
located  Is  read  Into  an  array.  Data  Is  processed  In  6400-byte  blocks. 
Each  block  Is  searched  for  chain  headers  and  when  one  Is  found  It  Is 
compared  to  the  chain  number  array.  If  a match  Is  found,  the  entire 
chain  Is  copied  onto  a temporary  file.  Those  chains  for  which  a 
match  Is  not  found  are  copied  to  the  output  buffer. 

Upon  completion  of  the  chain  search  the  actual  editing  begins. 
Editing  Is  performed  on  a chaln-by-chaln  basis.  For  a given  chain 
the  user  specifies  the  position  and  coordinates  and/or  rank,  when 
applicable,  of  the  point.  The  new  data  Is  sent  to  the  output 
buffer  and  the  next  position  to  be  edited  In  that  chain  Is 
specified  If  there  Is  one.  If  not,  the  remainder  of  the  chain  Is 
copied  to  the  output  buffer.  The  process  Is  repeated  for  each  chain 
In  the  temporary  file. 

Comnon  Blocks 

Block  Name  Contents  Description  of  Contents 

/BLKl/  lOUT  output  device  number 

/BLK4/  USTDP,  RUSTDP,  CTRLME  constants  for  conic  pro  - 

RGCOBE,  PI,  PIl,  COPST  Jection 

Subroutines 

MTBUFR:  copies  output  buffer  to  tape  and  resets  the  pointer  to 

the  top  of  the  buffer. 

CALL  HTBUFR  (lOUTB,  I) 

lOUTB:  output  buffer 

I:  pointer  In  the  buffer 
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PROJEC:  Projection  of  point  from  the  globe  onto  the  plane  via  a 
secant  cone. 

CALL  PROJEC  (RLAT,  RLONG,  X,  Y) 

RLAT:  latitude  of  point  In  radians 
RLONG:  longitude  of  point  In  radians 

X:  projected  x-coordlnate,  returned  by  routine 
Y:  projected  Y-coordlnate,  returned  by  routine 

See  SUBSET  documentation  In  this  section. 
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Figurt  28:  EOiTDB  FLOW  CHART 
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OI>«ENSION  I^Bn60fn»  I0liTRM6M^1, 

OImCNSTON  iafcV(4> 

DIMENSION  XIN(?),XOUT(?),KS{?),LIN»<(?«) 

COHMOM  /RLKl/  lOUT 

C0MM0»(/PLX4/UST0P,RU8T0P,CTRI MC,R6L0BE»PI,PI1,C0N8T 
EOU!V*LFNCCnNR(n,*IN(n).(IOlJTR(l  ),XOUT(D),  (IS(n,X5(D) 
rouTVALE^icr  (inn(n»T*<T(n) 

INTtGfR  TASK, TVPE, blank, GCTPns 
integer  MHrRUP,PCI9rTN,RANK 
r 

RGLOBC-J 
PI«3, 14159 
PIl-PI/lBP, 

IIST0P»57>J1 
BST0P»41,PP 
CTRl Mr»7,47 
CTRL*<t«CTRLMF*Pll 

RU3TPP«?,*PI*R6LOBF*((USTOP-B8T9P>/3«0,)*((CO8(U9TnP*PTn)/ 

c(ro5(B9Tnp*pii)-cos(iiSTOP*Pin)) 

CNTRLP«Ei»8TOP^RSTt>P)/2, 

CONST»R|iflT9P*?,*PI*NGLOnE*(IISTf5P.CNTRLP)/3flP, 

r 

1 WRITF(5,Q991 
RfA0(5,9«P)  IN 
wPTTF (5,99R) 

READ{5,9PP>  lOUT 
«RITEf5,997) 

RCA0(5,9PP1  T8K 
C 
r 

BL ANKsB 
r.E  TPns»i 
T»1 

IMFAf>«l 

n«i 

PRITFf5,994) 

L»» 

3 PFA0(5,9Pn  LINKfU 
IF rLlNK(l ) .FQ.P)  GG  to 
I •!  ♦! 

r 

GO  TO  3 

r ?t4pt  prccfssingi 

C f.PFATE  TF-P  FK  t FOR  TO  BF  COITFO 

IP  l.»)  -1 

TF(l ,LC,P1  GO  TO  5B0 
? RFAD(TN)  TNR 

on  ?p 

GO  TO  (fP1,lP?.lPJ,lP4,lP5,U«6l,lMF4n 

r 

UM  TF flNPfMi^GT.IfnP)  uo  TO  4PP 

ISAVf t)«TNR(M) 

TS*v(?1»Tnr(M*i) 

TSAvrn«TNO(H*7) 
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I3*V(4)«TNB(M*1) 

IHe*0a2 
r.O  TO  20 

102  DO  40  Jal,L 

TF(INB(Man,EO,LINt<(J))  GO  TO  104 
40  CONTIMJC 

IOliTB(I)«I»AV(l) 

I0UTB(!*1)«1SAV(2) 

imiTB(T*2)»ISAV(3) 

I0UTB(U3)aT8AV(4) 

I*I^4 

TFfI,GE,159B)  CALL  0T8UFB ( lODTB, n 
10I)TB(T)«1NB(H) 

IHEAD*3 

IOUTP(!*n«INB(M*n 
21  I0UTB(I*2)«INB(M42) 
inilTP(I*3)aINflfH*3) 

T»I*4 

1T(I.GE.1S98)  CALL  «TBUFB ( lOUTB, T ) 
GO  TO  20 

103  xniiTniixTNTB) 

xnuT(i^n«xiN(M^i ) 

ir  (XTMTMl  ,rO, 11*00, 1 IHEAOal 
on  TO  21 

104  1ST  in»I8AV(l) 

1ST II*l)«T8AVr2) 

IS(TI*2)»ISAV(3) 

TSni*3)»T8AV(4) 

IIaIT*4 

TF  f IT.I.T.1BQS)  GO  TO  107 
TI«1 

«BTTF(I8K)  is 

105  CP*iTT*UF 

107  ISdDaTMBf*-) 

IMF A0»6 

isfii^i)«iNB(M>n 
104  IS(II*2)»TMB(M*21 
!S(II*31aTNB(H*3i 
T !•! 1*4 

tf  fTT.LT,lB9B)  GO  TO 
rial 

ttPTTfnSF)  IS 

r,p  TP  ?'* 

lOA  XS(II)aKIW(M) 

*5(II*1  IallT^T>'*n 

TF f *IM(K) 10PP.1  THFAOal 

r.P  TP  i«B 

20  rpi'TiMiif 
GP  TP  7 
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4W0  *sf I n ■! 

KRTTPnSK)  IS 

c 

C St*9T  rOTTlWR  nw  CW4TW  Br  CW4JW  BASIS 

c 

WRITE  f5,995) 

WRITE  f5,97B) 

OP  4B1  J«t,L 
CALI  RFwnSK) 
r,FTPOS«1 
WHfRIJRX* 

^L-LI^iKf  J) 

WRITE  (S,9BPI)N|. 

A(13  IhF40»1 
C 

4?0  RFADdSKI  IS 
pn  4P? 

r.p  TP  (6B1 ,6P2,6«3,6B4>  , IRCAP 
r 

601  IMEAO*? 

IF  f IS(W)  ,GT,  ►’^HSr  IBM 

LS»Jt«I3(M) 

LSN2.IS(H41 ) 

TYPE ■I3(M4?i 
r.n  TP  AM? 

r 

to?  TWEAn«.l 

IT  fI5(w*n.NF,^iLT  r.O  TO  4B2 
THEAPaA 
inUTRf t)al SN1 
70UTB(I*l)tLSN? 

TOUTBf I*?)bTVPE 
IPHTB(T*3)aRl AMK 
laT^A 

rF(r.Gt,i59fli  CALI  'ttbuerttputb, n 
iPi)TPn)»is(><) 

Tni'TB(i*n»is(H*n 
imiTB( I*?)bIS(w*2) 
inilTBf  U1)aIS(“*,3) 

TaUA 

im,Gt.1S9M  CALL  *’TBlirR(ir)IITn,TT 
r,n  TO  AB? 
t 

603  IFfXS(W).EQ,1000.)  iHFAOal 
r.p  TP  ab? 
r 

6BA  whERIJHawHE  RURa  1 

IF  roETPos.rQ.Bi  r.n  to  rb*. 

WRITF (5,oq?) 
t>FAn(R,6Qn  POSITM 
Ad*'  wRI  TE  (6,R9a)  POSTTN 
REAnfR.Oqil  TASK 

Tf  f f P.r.T  , T 4SK)  .DR,  f 6.1,  T.TiStn  1 GO  TO  4BR 
Tf  Oi.l»nna  1 
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CFTPOS*0 

60?  ir((WHtRUB.EO,POSITN),*NO,(IFOlJNO.FO,l))  CO  TO  607 
XOliT(I)»yS(H) 

xouT(i*n»xs(M^n 

lOlJTBf  I^?)«IS(M42) 

I0UTR(I*3)»I8(M*J) 

I*I«4 

IF(I,GE.199B)  CALL  MTSUFR ( TOUT0, T ) 
!FfX9(H),EQ.1000,«)  GO  TO  4H6 
GO  TO  401? 

607  IF((Ta9K,E0,2).0R,(T*SX,E0.R))  CO  TO  6(99 

w»ITF(5,9901 
PEA0(5,9891  RL*T,RL0NG 
CALI  PP0JEC(PL4T,RL0nG,X,YT 
XOUT(I)»X 
xnuT(i4ij»v 

IF({TASK.E0,3).nR,fTA3K,FQ,4))  GO  TO  689 
TnilTB(l4?)«I3(H*?) 

IOUTn(t*3)«lS(M*31 

T»I*4 

IFf  I,GE,159B)  call  OTBUFPCTO'iTB.n 
GO  TO  613 
C 

609  TF (TASK.FO.R)  GO  TO  613 
WRITF  f5,9fl«) 

PFAO(R»9081  RANK 
inilTP(T43)«R*NX 

T»r*4 

tF(t.GE.lR9B1  CALL  mTOHFR(  lOIITP,  I) 

IF fTASK.Nf.A)  GO  TO  613 
C 

M«M_4 

r. 

613  IF { TASK ,F Q.s)  wHFRMMbKMEPDP-I 
IFOUND*!* 

»“RITF  (5,966) 

PEA0(5,993)  IPACK 
IF r IHACK.EO, 1 ) GFTPOS«l 
4 0?  CONTTnijF 
r.n  TO  4?0 
r 

406  T*"  riFnijNO.Fn.0)  GO  rn  401 

KRITF (5,965) 

4(11  CONTiNiJf 

xnnT(T)»l880, 

KPT  TE  ( iniiT)  IOIJT6 
500  STOP 

999  format (?9HFNTE»  INPUT  OCVICF  numrer(NN)) 

900  FORMATd?) 

996  FORMATOPHrNTFR  OUTPUT  OFVICF  NI|M«EP  ( NN ) ) 

99/  rORMAT  (31MFNTF.R  SCPATCw  OIVICF  niimRFR(NN)  ) 

996  1 ORmaT(31MFnTLR  CHaInS  TO  PE  EO I TFO ( NNNN ) , / , 

1 40HFMF9  7FRO  AS  F Nf)  np  CHAINS  TO  Pf  FOITED) 

901  rr)PMAT(T4) 
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99?  ro»M*Tf4?M  fiKTFB  POSITION  OF  POINT  TO  BE'  FD I Tt D (NNN)  ) 

991  FOPN*Tn3) 

98tt  FOPNATMBH  CHAIN  NUHREP  ■ ,13) 

994  F0«haT(19h  POSITION  NUMBER  • ,13,??M  ENTER  T*SR  HUMBER  (N)) 
993  FORMAT(II) 

995  FORmAT(13h  fOIT  COOES  l,/,21H  1-CHANGE  C00R0IN4TFS,/, 

12?H  ?-CmanGE  detail  level, /,7H  3-H0TH,/,I?M  4.AD0  POINT) 

986  FORMAT(31H  FURTHER  EDITING  IN  SAME  CHAIN?,/, 

118H  FnTER  1»YE8,  B»N0) 

986  F0RHATM9H  POSITION  NOT  FOUND) 

988  FOBmAT(?4H  enter  OETaII.  LEVEL  fNN)) 

989  rORH4T(E14.e,2X,E  H.B) 

99b  rORHATT47H  ENTER  L*T  AND  LONG  IN  RADIANS  (F 1 4 , 8, 2X,F t 4 . 8 ) ) 
9/B  rORMATflSH  5-DELETF  POINT) 

F NO 


1 

BL41 

BBB4 
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lOtIT 

BBFP 
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PLX4 
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PUSTOP 

BBB4 
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ctrlhf 
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PCI  DBF 
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PI 
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PP14 

r 
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?P1  fl 
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.u 
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I NR 

1 ;tc 

r 
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t 
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r 
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XOUT 

311/iC 

r 

T3 

4 AUC 
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r 

ISAV 
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636C 

r 
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BLANK 
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63/P 
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[ 
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90P 
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( 
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A 
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I 3- 

6jr( 
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f 
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63F4 

c 

IMCAD 

63F0 

F 

11 

63FC 

A 
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0F8? 

c 
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3 

OlPfl 
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ntbufr 
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MERGE 


Purpose 

MERGE  combines  dats  flies  (e.g.,  coastline  and  boundary  data) 
sequentially. 

Procedure  Description 

A data  base  Is  copied  to  the  output  tape  In  6400-byte  blocks . 
Each  block  Is  checked  for  the  end-of-flle  mark,  and  the  block  con- 
taining this  mark  remains  In  core.  The  next  file  Is  copied  beginning 
just  before  this  mark.  The  chain  numbers  of  the  second  data  base  are 
Increased  to  reflect  the  new  positions  of  these  chains  with  respect 
to  the  previously  copied  data. 

Common  Blocks 

None 

Subroutines 

None 
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APPENDIX  A 


DATA  BASE  FORMATS 


Table  I 


World  Data  Bank  I 

Block  Size:  1600  bytes 

Record  Size:  80  characters 

File  1 and  File  2:  Coastline  and  Boundary  data  files 

Field  No.  Field  Length  Data  Description  Format 

(Bytes) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

File  3: 


Map  Area 


7 

Line  Segment  Number 

17 

20 

Latitude  In  Radians 

E20.8 

20 

Longitude  In  Radians 

E20.8 

3 

Blank 

3X 

2 

Latitude  - Degree  part 

12 

2 

Latitude  - Minute  part 

12 

2 

Latitude  - Second  part 

12 

1 

Dlrectlon-N (North) ,S(South) 

A1 

3 

Longitude  - Degree  part 

13 

2 

Longitude  - Minute  part 

12 

2 

Longitude  - Second  part 

12 

1 

Dlrectlon-E(East) ,W(West) 

A1 

3 

Blank 

3X 

2 

Rank 

A2 

1 

Blank 

IX 

9 

Record  Sequence  Number 

19 

The  third 

file  of  WDBI  Is  actually  a merge  of 

two  files — 

the  Map  Area  Code  Index  and  World  Data  Bank  I Index.  The 
format  of  each  Is  given  below. 

Code  Index 
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FIlMi^D 


Field  No 


Format 


(Bytes) 


1 

11 

Map  Area  Code 

All 

2 

1 

Blank 

IX 

3 

24 

Map  Area  Description 

A24 

4 

World 

44 

Data  Bank  I Index 

Blank 

44X 

Field  No 

. Field  LenKth 

(Bytes) 

Data  Description 

Format 

1 

11 

Map  Area  Code 

All 

2 

9 

Blank 

9X 

3 

1 

Map  Feature 
B “ Boundary 
I - Island 
C - Coastline 
L “ Lake 

A1 

4 

4 

Blank 

4X 

5 

1 

Rank 

1 ” appears  on  all  aaps 

2 • single  point  Islands 

A1 

6 

4 

Blank 

4X 

7 

7 

Line  Segment  Number 

17 

8 

43 

Blank 

43X 
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Table  II 


Reduced  World  Data  Bank  I 
Block  Size:  1600  bytes 

Record  Size:  16  byte  binary 

File  1 and  File  2 : Coastline  and  Boundary  data  files 
Field  No.  Field  Length (Bytes)  Data  Description 

1 4 First  three  digits  of  Line 

Segment  Number 

2 4 Last  four  digits  of  Line  Seg- 

ment Number 

3 4 Latitude  In  radians 

4 4 Longitude  In  radians 


177 


Table  III 


Projected  Data  Base 

Block  Size:  1600  bytes 

For  each  chain: 

Field  No.  Field  Length 

(Bytes) 


1 

2 

3 

Chain  ID 
Field  No. 

1 

2 

3 


4 

5 

6 
7 


32 

. Variable 
16 

Field  Length 


4 


4 

4 


4 

4 

4 

8 


Data  Description 

Chain  ID  (q.v.) 

Data  points  (q.v.) 

End  of  chain  (q.v.) 

Data  Description 

First  3 digits  of  original  line 
segment  number 

Last  4 digits  of  original  line 
segment  number 

Type  of  chain 

1 - coastline 

2 - boundary 

3 “ box 

4 “ other 

Unused  integer  field 
Original  source  chain  number 
Cumulative  chain  number 
Unused  integer  field 


Data  Points 
Field  No. 

1 

2 

3 


Field  Length 
(Bytes) 


4 

4 

4 


Data  Description 

X - coordinate 
Y - coordinate 

Deviation  of  point  from  its  trend 
line* 

Rank  of  point* 


this  value  is  determined  by  DETAIL 
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End  of 

chain 

Field  No. 

Field  Length 
(Bytes) 

Data  Description 

1 

4 

1000,0 

2 

12 

Unused  Integer  field 
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Table  IV 


WDBI  Index 

by  Line  Segment 

Number 

Field  No. 

Field  Length 
(Bytes) 

Data  Descrintion 

Format 

1 

3 

First  part  of  Line  Segment 
Number 

13 

2 

4 

Second  part  of  Line  Segment 
Number 

14 

3 

2 

Blank 

2X 

4 

1 

Map  Feature 

A1 

5 

2 

Blank 

2X 

6 

24 

Map  Area  Description 

6A4 

7 

2 

Blank 

2X 

8 

40 

Map  Area  Code 

10A4 
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APPENDIX  B 


COMPARISON  OF  MAP  SCALE  DISTORTION  BY  PROJECTION  TYPE 
REGION  DESCRIPTION 

Assume  the  region  to  be  mapped  has  a 40°  latitudinal  span  and  a 
60°  longitudinal  span  (a  somewhat  larger  east-west  extent  than  the 
European  region  we  map) . The  location  of  the  region  will  vary  among 
projections  In  order  to  simplify  calculations  and  to  minimize  distor- 
tion. 

For  the  zenithal  projections,  the  region  Is  assumed  to  be  cen- 
tered at  the  North  Pole.  The  scale  error  of  a point  whose  position 
is  measured  relative  to  the  point  of  tangency  to  the  plane  is  unef- 
fected by  the  absolute  location  of  either  point.  Thus,  for  example, 
a point  20°  south  of  the  North  Pole  with  the  projection  plane  tangent 
at  the  pole  will  undergo  the  same  scale  increase/decrease  along  a 
radial  line  from  the  point  of  tangency  as  a point  in  latitude  30°  N if 
the  projection  plane  were  tangent  along  the  parallel  50  N (in  the 
figures  below  NPl'  - TP2').  In  the  former  this  radial  line  Is  a me- 
ridian. 


N PI  ' 

1 

r P2' 

^ — 

N 

/pi\ 

.A 

/ \ 

/ \ 

50^/^^ 

1 

V \ 70  \ 

( \ 

\ 

E / 

o 

o 

\ 

\ 

V 

1 

e\^ 

V 

N = North  Pole 

T is  the  latitude  50°  N 

EQ  = Equator 

P2  is  in  latitude  30°  N 

PI  Is  in 

latitude  70°  N 
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In  the  cylindrical  group,  the  region  is  centered  along  the 
equator  for  both  the  Mercator  and  Simple  Perspective  Projections  so 
that  the  normal  position  of  the  cylinder  is  used  and  calculations 
are  minimized.  For  Gall's  projection,  the  area  is  centered  near  the 
45°  N parallel  as  the  cylinder  intersects  the  globe  at  45°  N and  45° 
S parallels  and  thus  the  minimum  distortion  of  the  region  is  calcu- 
lated. 

In  both  cases  of  the  conical  group  for  which  scale  distortion 
is  given,  the  region  is  centered  at  the  49°  N parallel.  Thus,  the 
tangent  cone  is  tangent  along  the  49°  N parallel  and  the  secant  cone 

O O 

intersects  the  globe  at  41  N and  57  N parallels. 


COMPARISON  OF  SCALE  DISTORTION 


For  each  projection  two  types  of  scale  distortion  have  been 


determined;  meridianal  and  parallel.  The  percentage  of  parallel 
scale  distort  on  for  projections  in  the  three  classes  has  been  de- 

O 

termined  by  calculatin),  the  ratio  of  the  length  of  the  parallel  N 
away  from  the  'center'  of  the  map  to  its  corresponding  length  on  a 


globe  of  unit  radius.  The  'center'  of 

for  all  zenithal  projections  - 

for  tangent  cones  and 
cylinders 

for  secant  cones 
for  Gall's  projection 


the  map  is  defined  as  follows: 

point  of  tangency  to  the  plane 

intersection  of  the  central 
meridian  with  the  standard 
parallel 

Intersection  of  the  central 
meridian  with  the  central 
parallel  (mid-parallei  between 
the  two  standards) . 

Intersection  of  the  central 
meridian  with  the  45  N par- 
allel 


The  percentage  of  merldinal  scale  distortion  is  a projection/ 

globe  distance  ratio  and  is  determined  as  follows: 

zenithal  projections  - meridlanal  distance  from  the 

center  of  projection  to  a point 
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in  parallel  N°  away  from  the 
center 

cylindrical  projections  - meridianal  distance  from  the 

equator  to  a point  in  parallel 
N away  from  the  center 

conical  projections  - meridianal  distance  from  the 

equator  to  a point  in  parallel 
N°  away  from  the  center 

Tables  V a and  b list  the  scale  distortion  for  the  stereographic, 
orthographic  and  gnomonlc  cases  of  the  zenithal  ^roup;  Tables  VI  a 
and  b lists  the  same  for  three  cases  of  the  cylindrical  group  and  Ta- 
bles VII  a and  b for  two  conical  projections. 
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Table  V 


ZENITHAL  PROJECTIONS 


(a)  % MERIDIANAL  SCALE  DISTORTION 

Projection 

Stereographic 

Orthographic 

Gnomonic 

(b)  % PARALLEL  SCALE  DISTORTION 

Projection 

Stereographic 

Orthographic 

Gnomonic 


Distance 

from  Point  of 
(Degrees) 

Tangency 

20“ 

30° 

1 % 

2.3  % 

2 % 

4.5  % 

• 

4.2  % 

10.3  % 

Distance 

from  Point  of 
(Degrees) 

Tangency 

10° 

20° 

0.8  Z 

3.1  % 

0 % 

0 z 

1.5  % 

6.0  % 
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Table  VI 


CYLINDRICAL  PROJECTIONS 
(a)  % MERIDIANAL  SCALE  DISTORTION 


Projection 

Distance 

from  Standard 
(Degrees) 

20° 

o 

o 

Mercator 

2.1  X 

4.9  % 

Simple  Perspective  (one 
Standard  Parallel) 

4.2  X 

10.2  X 

Gall's  (two  Standard 
parallels) 

13.3  X 

14.2  % 

(b)  Z PARALLEL  SCALE  DISTORTION 

Projection  Distance  from  Standard  Parallel 

(Degrees) 


10° 

20° 

Mercator 

1.5  X 

6.4  X 

Simple  Perspective  (one 
Standard  Parallel) 

1.5  X 

6.4  % 

Gall's  (two  Standard 
parallels) 

15.8  X 

21.9  % 
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Table  VII 


CONICAL  PROJECTIONS 
(a)  % MERIDIANAL  SCALE  DISTORTION 


Projection 

Distance 

Parallel 

from  Standard/Central 
(Degrees) 

20° 

30° 

One 

Standard  Parallel 

0 % 

0 % 

Two 

Standard  Parallels 

0 % 

0 % 

PARALLEL  SCALE  DISTORTION 

Projection 

Distance 

Parallel 

from  Standard/ Central 
(Degrees)* 

9° 

19° 

One 

Standard  Parallel 

1.2  % 5.7  % 

Two 

Standard  Parallels 

0.2  / 

i 3.8  % 

^ o o o o 

9 and  19  are  used  Instead  of  10  and  20  because  distortion  at 
this  distance  had  been  previously  determined  for  the  conic  projection 
with  central  parallel  at  49°  for  the  European  map  displays. 
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